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INTRODUCTION 


The problem of motion in gravitational theory ■was first solved 
in a paper by Einstein, Inf eld and Hoffmann in 1938. The calcula¬ 
tions were so troublesome that we had to leave on reference at the 
Institute for Advanced Study in Princeton a whole manuscript 
of calculations for others to use. After that, Einstein and I made 
some progress together on this problem. Twenty-two years have 
elapsed since the first paper was published and I have again worked 
on this problem with my students in Warsaw during the last few 
years. This book presents the final results of all our work. Parallel, 
independently of us and a little later, W. Fock and his school in 
Leningrad tackled and solved the problem of motion in relativity 
theory, too. His results arc also piesented in his book entitled 
“Theory of space, time and gravitation” Pergamon Press, London, 
Hew York 1960. Though our approach is different from Pock’s 
and more in the spirit of Einstein, this work is not intended to 
be polemic. 

I have written this book with Dr Jerzy Plebahski. We discussed 
the contents carefully over the four years it took to write it. Un¬ 
fortunately we finished only the first chapter and appendix when 
Dr Plebahski received a Rockefeller Fellowship to go to the 
United States. Before leaving, Dr Plebadski prepared a sketch 
in Polish of the rest of the manuscript with the exception of the 
last chapter. This presentation was later very much changed by 
me for which I take the full responsibility. 

This book presupposes only a knowledge of the general principles 
of relativity theory. Readers of greater mathematical inclination 
are advised to read the appendix first and not merely the short 
chapter on notation which is a summary of it. 



8 


INTRODUCTION 


In. -writing the book, w© were greatly helped by Dr Andrzej 
Trautman -who made many critical remarks, checked the formulas 
and prepared the bibliography. Our thanks are also due to Dr 
W. Tulczyje-w who helped me greatly in preparing the last sections 
of Chapters TV and V. 


Warsaw 1960 


Leopold Infold 



NOTATION 


A. NOTATION OF GENERAL RELATIVITY THEORY 

We shall use throughout tho tensor analysis of General Relativ¬ 
ity Theory (G. R. T. for short). 

Wo shall denote by 

:c°, a? 1 , a»*, as* (0.1) 

the time and space coordinates of a Riemannian manifold. If we 
assume Special Relativity and a Cartesian coordinate system, then * # 
corresponds to the time t from 

x° = at, 

c being the velocity of light; for k —. 1 , 2 , 3. the x k (or x) denote the 
space coordinates. 

All Greek indices run from 0 to 3, Latin indices from l to 3. 
Repetition of indices implies summation. 

Tho geometry of tho Riemannian space-time continuum is char¬ 
acterized by a symmetrical metric tensor 

(«- 2 ) 

To distinguish between time and space in all possible coordinate 
systems we must assume that the metric tensor always satisfies 
the condition: 


iloo •''' b, 



0 

000 / 


(0.3) 


for arbitrary y a ^ 0. 
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Ins tead of these, we may assume the equivalent Hilbert condi¬ 
tions restricting the arbitrariness of space-timo transformations: 


0oo > 0, 





9 oo ? 

ffon 

002 

9 oo y 

001 

<0, 






010 > 

9m 

018 


011 








P 20 ? 

9 211 

0a 2 


(0.3a) 

0, g = detail < 0. 


The metric tensor g aP is a generalization of the Minkowski 
metric tensor ^ of Special Belativity Theory, defined by 


Voo — Voa — 0 , Vab — ^06 — 


for 

for 


a — b, 

a b. 


(0.4) 


To the covariant metric tensor g^ there corresponds a contra- 
variant metric tensor g afi defined by 


irg* = n = 


for 

for 


a = /5, 

a ¥* P- 


(0.5) 


We shall denote the determinant of g^ by g and all quantities 
that transform like 

X tensor 

we shall call tensor densities and denote by 

s S!::a - M 

The ordinary derivative will usually be denoted by a stroke: 

dS 


S ' a ~ dx a ' 


(0.7) 


The Ohristoffel symbols, which do not have tensor character, 
are: 


[tf> %{9a Y \p + 9fly\* 9afi\y) 7 

(0.8) 

II 

NJ 

/o 

1—l 

(0.9) 
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Those symbols allow us to differentiate tensors in a covariant 
way. We shall denote such a covariant differentiation by a semi¬ 
colon: 


jf.a. 




..a.;/? 



( 0 . 10 ) 

( 0 . 11 ) 


The indices written after the semicolon have tensor character and 
can be shifted up or down according to the ordinary rules. 

From the Ohristoffel symbols we form the full Biemannian 
tensor: 


From it, by putting p = a we form the contracted Biemannian 
tensor (Eicci tensor): 

.+M W -- WW- (013) 

and the curvature scalar: 

It = (0.14) 

The Einstein tensor is: 


d afi = B afi -\g ae B. (0.15) 

All the quantities g^, | a ^|, etc. are of a geometrical 

character. With at least some of them one associates a definite 
physical meaning. In a certain sense the Ohristoffel symbols repre¬ 
sent the intensity of the gravitational field and the metric tensor 
its potentials. We shall have much more to say about this later. 
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B. THE 8 FUNCTIONS 

We shall introduce “good” 8 functions — that is, 8 functions 
which, besides having the properties of ordinary Dirac <5 functions, 
also satisfy the condition 

= 0 for p = 1,2, ...,1c. (O.lti) 

All calculations in which <5 functions appear and no other assumption 
is made, are performed with “good” 8 functions. A more extensive 
theory of these functions is given in Appendix 1. 

C. THE FIELD VALUES ON THE W011LD-L1NES 

A 

Let us call £ a (s A ) the world-line of the A’tli singularity. Usually 
we take as the parameter on which the world-line depends, not the 
“eigentime” s A , but the time *° = el. Then the motion is eharaetor- 

A 

ized by | TO (®°) and f* = ar. Let us assume 

Then the “twoedling process” belonging to the A’th singularity 
and denoted by ~ written over the “tweedlod” expression means 

A 

two things: first, the singular part at x = f of / is ignored; second, 

it introduces the expression instead of x k into the regular part 
of f. Or, using our “good” <5 fmictions wo define: 

£ A 4 

/ = J dxd(x—i (x 0 ))/(®«, «*, h. C0.17) 

Prom this definition follows (if wo omit tho A’s): 


Vyt — c P\$ J r c P\a£ > \<i — V|a £ tt |o J 


(0.1.8) 
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We shall assume throughout that the functions with which we deal 
obey the rule for “twoedling” the products, which is: 

<py — (prp . (0 19) 

For more about this process, see Appendix 2. 


P. THE COVARIANT CHARACTER OP THE <5*8. 

TENSORS ON WORLD-LINES 

The four dimensional Dirac d (4) fnnotion is a scalar density. 
This follows from the relation: 

/ A w (x)da> = 1. (0.20) 

By “world-lino tensors” we shall understand tensors defined 
only on the world-lines £ a . To such tensors we may apply the rules 
of tensor algebra but not those of tensor analysis. To apply the 
latter, we must have tensor fields. We may change a tensor defined 
along a world-line, at least symbolically, into a tensor density field, 
in the following way: 

J^V)(*-^(A))%::^(A), (0.21) 

—00 

The transformation properties of <5 (3) follow from: 

OO 

<5(3) = f S {i) d£° (0.22) 

— OO 

which, together with (0.20), gives: 

fd (g) dx = l. (0 23) 

This means that S (3) dx can be treated as an invariant. If the (5 (3) ’s 
are our “good” 5’s, we may regard (0.22) as the definition of the 
“good” S (4) . 

For more about this, see Appendix 3. 



CHAPTER X 


GRAVITATIONAL INTERACTION AND THE GENERAL 
THEORY OE MOTION 


1. GENERAL STATEMENT OP THE PHYSICAL PROBLEM. 

PARTICLES AND THE GRAVITATIONAL FIELD 

We shall now state our problem in general terms: it is that of 
particles, and of the mutual gravitational interaction between 
them. 

The particles move in a Riemannian continuum. Its geometry 
depends on the positions and motions of the particles; this is 
determined by the distribution of matter. Such is the picture that 
we usually have before us in G. R. T. Is it the right one? In this 
. section we shall try to sketch an answer to this question. 

First, a word about the concept of point particles. Obviously 
this concept is an invention of the human mind. With its help we 
describe the much more complicated reality. Its use implies that 
in a given case we may neglect the dimensions of the bodies and 
their structure, if we wish to describe only the essential character¬ 
istics of motion. 

For this simplification we must pay a price. The picture of 
matter is simple, but that of the field becomes complicated. The 
field becomes singular on the world-lines! Yet the origin of this 
difficulty is of a non-physical character, therefore we are justified 
in treating it non-physically — that is purely formally. In doing 
so we shall use the mathematical tools developed in the Appendix, 
and sketched briefly in “Notation”. 

Later we shall have much more to say about the difficulties 
connected with the singularities of the field. They will turn out 
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to bo important in what we are doing from the formal point of view, 
as our aim is to find the essential characteristic features of the 
motion of point particles. One should, however, not overestimate 
these difficulties from the physical point of view. They are not 
essential and they are not connected with the logical structure 
of G. R. T. 

If we accept any coordinate system with respect to which we 
wish to describe the motion of point particles, then its description 
consists in providing a a 

1°: The world-lines of the particles |° = f°(A), 

2°: The metric tensor g aP that characterizes the geometric 
properties of the continuum. 

Thus, to be able to find the world-lines and the metric tensor, 
we must know the laws that govern the motion and the laws that 
govern the metric field which in turn appears to be formed by the 
motion. 

This is a difficult problem which we shall have to approach 
gradually. To characterize the mechanism that connects field and 
motion in G. R. T., we shall start with some remarks concerning 
the field concept generally. At the end of this section we shall see 
which of these concepts can be applied to G. R. T. and which 
cannot. 

In Newtonian physics, the concept of a field of force, and more 
specifically that of the gravitational field, is not of primary im¬ 
portance. The Newtonian gravitational theory speaks mainly about 
the force acting upon particles, depending only on their distances 
apart and their masses. This force is uniquely determined by these 
distances and masses, and the Newtonian equations of motion allow 
us, at least in principle, to find the coordinates of the point part¬ 
icles as functions of time. The fact that the force between particles 
depends essentially on their positions is usually understood in the 
following sense: the action of one particle upon another propagates 
with infinite velocity. Or, as it is often put, Newtonian Physics 
introduces action at a distance. 

According to Relativity Theory, no action can be propagated 
with a velocity greater than c — the velocity of light. Thus the 
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field propagating from one point to another, largely because of 
Relativity Theory, becomes a physical reality, as real and as 
material as corpuscular matter. In all descriptions of nature to-day 
•we have two principal concepts: field and matter. Besides corpus¬ 
cular matter there exists the field with its laws of changes in space 
and time, which are of a deterministic character. These laws arc 
as important as, if not more important than, those governing corpus¬ 
cular matter. 

The presumption of the existence of such real fields with a finite 
velocity of propagation, implies some laws concerning their structure. 
The fields must be described by vectors or tensors in a four-dimen¬ 
sional manifold as functions of x a , satisfying partial differential 
equations of the hyperbolic type, so that Cauchy’s problem can 
be correctly set. This is the mathematical expression of the physical 
statement that the field is governed by deterministic laws. Also we 
demand of these equations that they express the physical fact that 
the field is caused by moving particles. Thus, in our mathematical 
description, on the right-hand sides of those equations some ex¬ 
pressions characterizing the sources of the field usually appear. 
It is obvious that Dirac’s <3 functions are a convenient mathematical 
tool for the charecterization of the sources and their motion if point 
particles are assumod. 

This is the dualistic view which wo usually accept in physical 
theories, for example Maxwell’s theory. That is, we assume the 
existence of both the field and its sources. But what about the 
motion of these sources? Let us state: no linear field theory can 
determine the motion of its sources except in trivial cases. This can 
be seen in the following way. Let us assume that the motion of 
a particle A is determined by the field equations. We now take 
another solution in which the motion of another particle B is also 
determined by the field equations. Then, because the equations 
are linear, there exists a solution in which both particles A and B 
move exactly as before — that is, as in the solutions in which only A, 
or only B, moved. Thus we arrive at an absurd conclusion: the 
motion of A is determined by the field equations but it is not changed 
by the presence, or motion, of B. Therefore, from the linear theories 



GENERAL STATEMENT Off THE PHYSICAL PROBLEM 


17 


1.1 


motion can. only be deduced in the trivial case of non-interacting 
bodies. 

Thus in every linear field theory of interacting particles we 
must add to the field equations the equations of motion of the 
sources of the field. In the case of point sources, these are ordinary 
differential equations of the second order, in which the acceleration 
of a particle is determined by the value of the field at the point 
where the particle is. 

Thus we say: the laws of propagation of a field are usually 
described by linear partial differential equations of hyperbolic 
type; the motions of its sources, that is, of the particles or charges, 
are described by ordinary differential equations in which the 
accelerations are determined by field functions at the points where 
the particles are at a given moment. It is this situation which we 
have in mind when speaking about close-range action. 

It would seem that all linear and non-linear relativistic theories 
should bo based on close-range action in the sense just described. 
Yet this is not quite true. In particular we shall have to change 
slightly the meaning of close-range action so that it will fit Einstein’s 
theory of gravitation. 

Einstein’s gravitational equations are hyperbolic differential 
equations for the metric field g^. The components of the metric 
field are also understood as the potentials of the gravitational field. 
The sources of the field form the tensor of energy-momentum; 
they appear on the right-hand side of the field equations and char¬ 
acterize the distribution and motion of corpuscular matter. There¬ 
fore, such a situation would seem typical, fitting perfectly the scheme 
described before, and very similar to the case of Maxwell’s equations. 
But there is one very essential difference. Whereas in the situation 
described up to now, the equations of motion have to be added to 
the field theory (e. g. the Lorentz equations of motion) this is not 
true for G. R. T. It is a non-linear theory in which the motions of 
the sources are determined by the equations of the field. Our aim 
in this section is to investigate provisionally just this peculiar 
situation. 


Motion and relativity 


2 
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For the sake of greater clarity, we shall try to say more precisely 
what we understand hy the words “equations of motion” in grav¬ 
itational theory or, for that matter, in any field theory. It will he 
convenient for onr purpose to distinguish among equations of motion 
of three kinds. 

1°. The equations of motion of the first kind arc the equations 
of motion for a test particle; we shall understand hy them the 
equations of motion of a particle on which a given external field 
acts, while the field does not depend on the motion of the particle 
itself. 

By a test particle we mean, therefore, a particle that can he 
ignored as a source of the field; this is so when the coupling of the 
particle with the field is infinitely small, so that we may neglect 
the modification of the field caused hy the presence of this particle. 
In the case of a gravitational field we mean hy a test particle one 
with infinitely small mass. 

According to the principle of close-range action, such an external 
field acts upon the test -particle; the second derivatives of its 
coordinates with respect to the proper time are determined by the 
external gravitational field. Since the particle itself does not part¬ 
icipate in the formation of the field, the field is regular along the 
world line of the particle. It suffices to introduce into the external 
field the coordinates of the test particle. 

In the early days of G. B. T., these equations of motion for 
a test particle were postulated as the geodesic equations derived 
from a variational principle. Thus in G. B. T. (in the case of a test 
particle), as in any linear dualistic theory, we had two independent 
sets of equations:.those of the field and those for the motion of this 
test particle. Both the field equations and the equations of motion 
can generally (in a gravitational theory as well as —■ for example — 
in Lorentz’s theory) he deduced from one variational principle: 

»{W,+Wx) = 0, (1.1.1) 

where Wj depends only on the field and Wj on both the field and the 
world lines of the particles. In the case of a test particle we may, 
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however, split this equation into two independent variational 
principles 

dWf+dWi = 0, <5Wj = 0. (1.1.1a) 

In the first we vary the field, in the second — the world-line. Essen¬ 
tially all that was known in the early days of G-. R. T., as far as 
motion was concerned, were the equations of motion for a test particle 
which were independently assumed. 

This situation lasted until 1927, when a paper by Einstein and 
Grommer appeared. In this paper, it was claimed that the equations 
of motion for a test particle do not need to be added to the field 
equations, that they are a consequence of the field equations. One 
of the reasons that it is possible to deduce the equations of motion 
for a test particle from the field equations alone is that the field 
equations are non-linear. Thus in the theory of G. R. T., we do not 
need to postulate any equations of motion for a test particle. The 
field equations are sufficient! 

2°. By equations of motion of the second kind we shall mean 
equations of motion for finite masses whose accelerations are 
functions of the field on the world-lines of the particles. Thus the 
field at the particle may depend on the mass of the particle itself 
and on its motion. But in the case of a point particle the self-field 
is infinite on the world-line of the particle. Therefore, we must find 
a definite meaning for the concept of the field at the particle. It is 
here that the tweedling process described in “Notation” and Appen¬ 
dix 2, will play an essential role. 

On the other hand, we know that the field remains undeterm¬ 
ined until the motion is determined. Therefore, the equations of 
motion of the second kind will have only a symbolic meaning, since 
they will contain an unknown quantity: the field at the point of 
the particle. To make these equations of motion definite, we 
obviously do not need to know all the field components everywhere. 
It is sufficient to know only the combination in whioh they appear 
in the equations of motion and to know them only on the world¬ 
line of the particle. 
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These general equations of the second kind can he found in 
practice in any field theory from a variational principle: 

<5 W x = 0 (1.1.2) 

in which Wj is essentially the same as before in (1.1.1a) with the 
sole difference that the field at the points on the world-line is not 
assumed to be given but to be deduced from the field equations. 

We shall show later that in the case of Q-. B. T. we can obtain 
the same equations of motion that we deduce from a variational 
principle just as well from the field equations. The practical value 
of these equations is, however, small since the field which enters, 
into them remains unknown; it must be deduced from the field 
equations. 

3°. The equations of motion of the third kind are found from the 
equations of motion of the second kind by introducing into them 
the actual values of the field. 

How is it to be done? The answer to this question is not simple. 
We shall try now only to indicate it ip general terms and shall do 
it in detail later. 

Let us go back to the case of a linear dualistic field theory. Let 
us assume that we have the field equations for arbitrary motion. 
In other words, we assume the case in which the field equations 
do not determine motion and in which the equations of motion 
are independently added to the field equations. It is in this case 
that we can solve the field equations for arbitrary motion. The 
field tensors — solutions of the field equations denoted symbolic¬ 
ally by / — depend on position and functionally on the world-lines . 
Thus we can write symbolically: 

/ = /(*“)(K (1.1.3) 

A. 

Then, putting into / in place of ®“, say, f°(®°) we find the field at 
the JL’th particle. Introducing these values for the field into the 
added equations of motion, we obtain a set of differential-functional 
equations which we shall call the equations of motion of the third 
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kind — for the world-lines of the N particles. Thus in these equa¬ 
tions the field is already eliminated and what remains are equations 
formed only of the |’s. Their integration would give us the motion 
of 2F particles. 

Generally speaking, such is the situation in every linear theory 
in which the motion of the suces is arbitrary. But we are mostly 
interested in' the equations of motion in G. B. T. In Einstein’s 
gravitational theory the situation is in some respects more compli¬ 
cated, in others simpler, than that sketched above. 

In what respect is the situation more complicated? 

In G. B. T. we do not have field equations with arbitrary motion 
and then, additionally, the equations of motion. The same set of 
equations that determines the field, also determines the motion. 
How, then, can we find the solutions of such equations, solutions 
giving at the same time the field and the motion? The field is deter¬ 
mined by the motion and the motion is determined by the field. 
What is the way out of this difficulty? We shall see later that there 
are two ways out. We shall sketch them only very briefly and 
simplify their description (for the moment) to such an extent as to 
make them not quite truthful. 

The first way is: change Einstein’s gravitational equations into 
an artificial system of different equations so chosen that it permits 
arbitrary motion. (This is done by the help of a non-physical field, • 
which we shall call the dipole field). Then find the field corresponding 
to these 'non Einsteinian equations. However, in the end we wish 
to return to Einstein’s equations. This can be done only by removing 
the non-physical field, and so restricting the motion to that pre¬ 
scribed by Einstein’s equations. Thus the procedure consists in 
finding the field for the non-Einsteinian equations and then, at 
the same time, rsstricting both the field and the motion so that 
both satisfy the equations of G. B. T. 

The second way is based on the approximation procedure. 
Boughly speaking this can be described in the following way. We 
start with the Newtonian field which, in the first approximation, 
satisfies the field equations of G. B. T. Then wo find the motion 
belonging to this field and determined by the field equations. It is 
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the Newtonian motion. Then, having done this, we can find the 
field in the next approximation. To this field there belongs a motion 
calculated (with the help of this field) to the post-Newtonian approx¬ 
imation. At each approximation level the existence of the equa¬ 
tions of motion forms a condition that makes it possible to carry 
the approximation procedure for the field one step further. In 
this way, step by step, we can determine field and motion with 
increasing accuracy. 

These are the ideas that we shall employ in our study of the 
problem of motion in G-. E. T. They make the case of the gravita¬ 
tional field more complicated than, say, that of the electromagnetic 
field. However, as we mentioned before, in one respect the situation 
is simpler than that for an electromagnetic or for some other linear 
field. 

In what respect is the situation simpler? 

We have said before that the equations of motion of the third 
kind should be (generally speaking) differential-functional equations. 
Or, in other words, differential equations of an infinite order, since 
to describe the f’s at an arbitrary moment, all derivatives of the f’s 
at that moment are necessary. Tet the gravitational equations of 
motion, at least up to the post-Newtonian approximation, are 
a set of 3 N ordinary differential equations of the order 2. Therefore, 
# their structure is much simpler than we could have predicted from 
the general case. 

How does this situation in G. R. T. fit our concept of close-range 
action? This concept was formed for linear theories in which the 
motion was not determined by the field equations. Its mathematical 
expression consisted of the hyperbolic field equations and the 
differential-functional equations for the motion of the third 
fcind. 

In G. E. T. we also have hyperbolic equations for the field, 
as in the linear theories. However, after using the approximation 
method, the equations of motion that follow from the field equations 
turn out to be ordinary differential equations of the second order; 
therefore they have a form similar to the Newtonian theory of 
motion, that is, as in the case of action at a distance. In this respect 
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the concept of close-range action formed on the basis of linear field 
theories does not fit G. E. T. 

Can these equations of the third kind be gained from a varia¬ 
tional principle in -which only the world-lines and not the field 
appear? We shall show that this can, indeed, he done for the equa¬ 
tions of motion up to the pcst-Newtonian order. 

We shall call such variational principles leading to the equations 
of motion of the third kind those of Fokker’s typo. 

Summarizing the essential point briefly: It is the field which is 
of primary importance in G. E. T. The field equations determine 
not only the field; they determine the motion as well. 


2. NEWTON’S THEORY OP GRAVITATION 

Newton’s theory of gravitation is a theory based on the concept 
of action at a distance. In it the gravitational field appears only 
as an auxiliary mathematical tool. Because of its simplicity, this 
theory is well suited to be an example illustrating the three kinds 
of equations of motion and the simple variation principles that 
lead to them. 

A 

Let |(t) be the world lines of JVpoint particles (A = 1 , 2 ,..., N), 

A 

m their masses and <p(x, t) the Newtonian gravitational field. We 
shall assume, as the foundation of Newton’s theory, the variational 
principle: 

SW[<p, I, I, 0, (1.2.1) 

where the action W is a functional of the field <p(x, t) and the world- 

A 

lines § (t). The variation has to be performed with respect to both <p 
and the |’s. The action W consists of two parts: W f belonging to the 
field and Wj belonging both to the world-lines of the particles and 
to the field, that is, to the field-matter interaction: 


W = Wf+Wi. 


(1.2.2) 
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We define W f and TTj: 

F,= 



Tc = gravitational constant) 


(1.2.3) 


N tz A A N a A A 

W x = y J dt^mf — y J dtjdxmd(x—§{t))<p{x,t). (1.2.4) 

jfcl A-l «! 

Lot ns start "with the field equations. They are obtained by 
varying W = W/+W r -with respect to <p. We have 
■v <2 A A 

dW = — J dt JdxmS(x— §)d<p(x, t) + 

aZ i 

+ / dtjdx^aadcp, (1.2.5) 

*1 

and from it: 

N A A 

<P\aa = J£md(x —|). (1.2.6) 

A-l 


IE we admit only solutions that vanish at infinity, then the solution 
of the above equation is unique: 


N 


<p = 


-'S' 


A 

m 


A-l X- 


A ' 

'll 


(1.2.7) 


Thus we can solve the Newtonian field equations with arbitrary 
motion, that is with arbitrary functions # [t). The field does not depend 
functionally on the world-lines. The field q> is simply a function 
of the position coordinates of the particles, all taken at the 
same moment of time. 

Now, as to the Newtonian equations of motion: 

The equations of motion of the first kind, or, the equations of 
motion for a test particle in a given gravitational field, are obtained 

1 A 

by putting in W x to = Am and to = 0 for all the other particles 
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(that is, for A =2,3,..., N). Assuming that <p is regular along 
the -world-line §(t) of the test particle, we take into account only 
the following part of W x : 

h 

W'j = j dt[\Am??-Amq>{£,t)). (1.2.8) 

h 

Varying W'i with respect to f, we have: 

SW' X = - J dtAmd£ a (?+y l8 ), (?,. =|pr) (1-2.9) 

from which we obtain, as the equations of motion for a test particle 
in a given gravitational field <p: 

Am£ a = — Amcp\ ai (1.2.10) 

and, dividing by Am, we have, finally: 

? - -9V- ( 1 - 2 . 11 ) 

The equations of motion of the second kind we obtain by varying 
the full expression (1.2.4) for Wj with respect to f(t): 

N h A a - A 

dWx = - Jg f dtm-d£*(£*- f dx3(x-£) ls <p(x,t)). (1.2.12) 

A -1 ti 

Integrating by parts under the integral, we obtain, because of 
d[dx = —d/di, the equations of motion of the second kind: 

* A A 

= ~9\b — — J dxd{x—§)<p\ 8 . (1.2.13) 

In these equations there appear the field intensities along the 
world-line. But <p may be singular along a world-line. We get rid 
of this singularity by using our good S functions. In the case of a test 
particle it was immaterial whether we used our good functions or 
Dirac’s 3 functions, since there their only purpose was to introduce 
the § in place of the x. In the case of a test particle it was also 

immaterial whether we used 5 >i* or <p la , since there q> did not depend 
on £. This is not immaterial when varying (1.2.12). In passing from 
(1.2.12) to (1.2.13) we must assume that 9 0 , at least in its regular 
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part, does not depend on f explicitly, or, if it does, that we vary 
only the <5’s and not the (p’s with respect to the f’s. 

The equations of the third kind are gained by substituting into 
the equations of the second kind the field variables as functions of 
the world-lines of the particles; this means that we have to introduce 
into equation (1.2.13) for <p the value (1.2.7). We have: 


<Pi 


7 -?“(*) 
= -5 s 

R-l X — f(t) 


(1.2.14) 


so that 


<P\C 




N „ A li 
B _ ga 


R-l 


B.S 
X—I 


m- 


B—1 
B,tA 


A B 3 

l-l 


(1.2.15) 


A A 

This equation is true, because d(x— f) replaces x by f (f) in the 
regular part of the function. We may remark, in passing, that this 
equation would have been true had we used, not our good 8 func¬ 
tions, but ordinary Dirac 8 functions. This would have been so 
because the integral of the singular part of <p 1(t would have vanished 
on account of the spherical symmetry of the ordinary 8 functions. 
Therefore, to obtain a regular expression for the field on the world- 
lines in Newton’s theory, it is sufficient to appeal to the spherical 
symmetry of the ordinary 8 function. We do not need to introduce 
our good 8 function. Thus, introducing (1.2.15) into (1.2.13), wo 
obtain equations of the third kind: 


a w a^ b 

r—kSi'-' 

B+A 


A B 3 * 

l~l 


(1.3.16) 


or multiplying by m, the more conventional form: 


a xr 

A - V-J AB 

m f = — Tc > mm 

Brn 1 

B+A 


A 

! a - 


B 

■! a 


A B 3 

l-l 


(1.3.17) 
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Now. we can easily give Fokker’s variational principle leading 
to (1.2.17). This is the variational principle, well-known from 
theoretical mechanics: 


U h , N a a N AB 

' * (2 * AfV+* y fe -^-)] = o, (1.2.18) 


where the variations have to be performed with respect to the 
world-lines. 

What is the connection between this variational principle and 
that for the equations of the second kind? The answer is: If into 
W = Wf+Wj we introduce as the field <p the concrete solution 


N 




A 

m 


A - 1 x— # 


(1.2.19) 


of the field equations, then we obtain the right expression for W F . 

To prove this theorem we start with W f . We have the following 
equality which we introduce only for the sake of variation with 
respect to q>: 

1 h 1 <2 

W 1 = dx / dt^<P\a<Pa\ = f dtj dx<p laa <p. (1.2.20) 

<i 

Substituting for <p from (1.2.19) and making use of the field equa¬ 
tions (1.2.6), we have: 


N N 


w, = 


— $ y y *hmm J'dtj'dxd(x—§) |x—f | -1 

A-l B—l <i 

N ti A B N it 

>. 2 J**-^-* 2 J**'f* 

§ — § 


Aqitjimml <2 


A ■“ 1 


— * ( 1 . 2 . 21 ) 
x—# 
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The last integral vanishes, because our 3’s are the good 3’s. 
Now substituting <p from (1.2.19) into Wj, we have: 


$2 N ^ ^ <2 N ^ 

W 2 = f + f <2t 2: Zc-m J dx8[x~ |) 

-4-1 <1 A,R-1 


R 

m 


R 

•# 


( 1 . 2 . 22 ) 


ti N A A h N 

= fdt Y im?? + f dt 2J 

<1 A-l U A t B~l 

* * A » 


A R 



f 2 J <Mm 2 J dx-y*- |y } . 


Again the last integral vanishes for the same reasons as before. 
Adding the last two equations together we find: 


W f +W x = 


it . N A A N AB < 

i \iA jg* I' 


(1.2.23) 


which is identical with Fokker’s Wj- which we quoted in (1.2.18). 
But this complete agreement is due to the use of our good 3’s. 
Using ordinary 3’s we should have found in (1.2.23) an additional 
expression of the form: 

li N 

= i f dt Y Jem 2 f dx-^ (1.2.24) 

HaZi |X| 

which is infinite in the case of Dirac’s functions, zero in the case 
of our good functions, and generally could take any pro-assigned 
value, since we can find models of the <3 functions for which 


/ 


I* 


= <u (1) 


where co (1) is an arbitrary pre-assigned number. In any case this 
additional expression would be constant and would have no influence 
upon the variation. Thus, as far as the variation of W& is concerned, 
our statement is independent of the choice of the 3’s. 
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The Newtonian theory of gravitation as summarized in this 
section uses the idealized concept of point particles. With its help 
we represent the essential characteristic feature of planetary motion 
on which the dimensions of planets and their rotation have a small 
influence. 

G. B. T. is essentially a generalization and a broadening of 
the simple Newtonian theory. There too, the introduction of point 
particles and the use of the 6 functions will allow us to describe 
the important features of their motion. 

The Newtonian theory has been developed in this section in 
such a way as to allow us later to stress the differences and, of course, 
also the similarities between it and G. B. T. The most essential 
difference consists in the fact that in G. B. T. the equations of 
motion follow from the field equations. Yet it is the Newtonian 
theory which will serve us as a guide to G. B. T. Indeed, we shall 
look for solutions which for e -> oo go over into solutions of New¬ 
tonian theory. The principle of correspondence is just as important 
in the relation between G. B. T. and Newtonian theory, as it is 
between quantum theory and classical mechanics. 


3. INTERACTION IN G. R. T. THE EQUATIONS OR MOTION 
OP THE FIRST AND SECOND KIND 


Wo have N point particles and their world-lines 1° = f°(A). 

A 

We denote their masses by Formally these may be regarded 
as the coupling constants of the particles with the gravitational 
field. 

For the interaction Wj between the particles and the gravita¬ 
tional field, we shall assume the following expression: 

W x - ^m {0 )0 f (gZpdPdP) 112 . (1.3.1) 

A*ml Ol 


We shall interpret g^ according to what was said in “Notation”, 
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or, more fully in Appendix 2, that is: put in f* instead of oc 7 *, ignoring 
the singular part of g a p. We can write the covariant definition: 

= f dx 3 (4) (oo—£W)9afi(.®) (1.3.2) 

■where 3, 4) is defined by (0.22), with the help of our good 3’s. 

If is regular in ®°, as is always the case, we may, again making 
use of (0.22) write the above equation in the form 

t>=S &«(*-W*!?)> ^- 3 * 3 > 


where the 3 is our good 3 function. We shall therefore assume that 
all integrals of the kind 



3(x) 

i*r 


for an arbitrary given p vanish. Finally, going back to the notation 
(1.3.1), we may remark that a y and <r a denote two arbitrary hyper¬ 
surfaces which are space-like, that is, such that the normals to 
them n a satisfy the condition > 0. The integrals must be 

taken along the world-lines between these two surfaces. 

The expression Wj defined in such a way (if we ignore for the 
moment subtleties arising from the definition of field components 
along the world-lines) is the natural generalization of the correspond¬ 
ing one in Special Relativity Theory. Indeed, in the case of Min¬ 
kowski space-time, that is, in the case according to (0.4) 

our equation (1.3.1) for Wj takes the form: 


Wj = - 21 m (0) e / (( dt 0 ) 2 - dt&ty *, (1.3.4) 

A-l 

that is, the well-known expression for inertial action in Special 
Relativity Theory. We see that our original expression (1.3.1) 
for Wj is the only sensible generalization of the above expression 
to the case of a curved world. 
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Thus, having Wj, we can easily obtain the differential equations 
of motion of the first and second kind. 

We start with the equations of motion of the first kind, that is, 
the equations of motion for a test particle. To do this, we have 
to vary in (1.3.1) only the expression with mass m, for which we 
assume a small value Am. We assume, too, that Am does' not influ¬ 
ence the given field, which remains regular along the world-line of 
the particle. Thus the <5 (4) in (1.3.2) does only one thing: it intro¬ 
duces the world-line, that is, |° in place of x a . Therefore W x for 
a test particle is, as far as its variation is concerned, 


Wj = -Amo J (g afi (S)d^d^) l l 2 . (1.3.5) 

o i 


We obtain the equations of motion in the usual way by varying W x 
with respect to Writing 


ds = (g^apd?) 112 

we have: 





(1.3.6) 


If <5|“ vanishes at o t and cr 2 , we obtain from the condition 
SWi — 0, the equations for a test particle 


<1 l a?\ , d? d r A 

ds \ 9t * ds ) ^ 9ma ds ds 


(1.3.7) 


where, of course, g a? and its derivatives have to be taken at the 

d I d^\ 

point £°. Performing the differentiation —I&,«-=—I in the last 

ds \ ds I 

equation, and raising the index a by multiplying this equation by 
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3 *®, we obtain: 


a*? a? dr 

ds 2 + W ds tin 


(1.3.8) 


where the Christoffel symbol has to be taken at the point I*. Using 
the symbol D for covariant differentiation, we have: 


DT a = dT°+{^}TW. 


We can, therefore, write (1.3.8) simply as: 

T 


D dr 

ds ds 


= 0. 


(1.3.9) 


(1.3.10) 


These equations in either form, together with the definition of the 
line-element ds, fully define the world-line T as a function of s. 
These are the equations of a geodesic. 

We can also treat Wi as a function of the upper limit of ct 2 . This 
means: we integrate along the geodesic world-line, from a fixed 
lower limit to an arbitrary upper limit which is the intersection of 
the world-line with an arbitrary spaoe-like hypersurface. Then Wj 
becomes identical with the function satisfying the Hamilton-Jacobi 
equation. In this case we have: 

8Wt 

m = - -Amcg^d? 

so that: 

dWr di fi 

W i\* = -Jp = (1.3.12) 


, (1.3.11) 

a 


_ dr dr 

Because g af — —- = 
p ds ds 

tion for a test particle: 


1, we obtain the Hamilton-Jacobi oqua- 


= (cAmf. 


( 1 . 3 . 13 ) 
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It may sometimes "be convenient to find the motion of a particle 
by finding the total integral of this equation. Differentiating it 'with 
respect to the parameters, we may find all the integrals of the 
equations of motion. In this way we may avoid the sometimes 
troublesome calculation of the Ohxistoffel symbols. 

We shall discuss now the equations of motion of the second type. 
We treat the field as unknown; later we shall see how to find 
it from the field equations. The g^ will, in general, be singular on 
the world-lines. We must therefore be careful with the variation 
procedure, so as to get rid of the singularities. 

We shall vary g a0 with respect to £°. We have: 

= -fdx <3 (4) (0 - f) w g a0 4 e 

= /d®S (4) (®-|)^ 1# 4 c =^ 4* . (1.3.14) 


where g^ has to be treated as a given function not to be varied with 
respect to the f’s. 

We introduce 

ds A = (^dpdp) 1 ' 2 . 


(It is convenient to write ds A instead of di.) We then proceed to the 
variation of W z : 


<5W r 


N 

-1 
A~ 1 


itsP 

ds A 




<n 


+ 


N o 2 


+f 

A—l ai 


ds A dP 



A dp 


, A. dP df 1 , 

ds A da A ( 


1.3.15) 


Motion and relativity 


3 
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A 

Prom the vanishing of the <5£ on a x and <r a we deduce the equa¬ 
tions of motion of the second kind: 


d 3 dp f ^ dp dp 

ds A 9a * ds A * gMa ds A 1*2 ~ °* 


(1.3.16) 


We may change the form of these equations, remembering that 




d 
ds „ 


£)g a ? 


A _ 

r a dP 

1^2 J £)\a9ap— 9a&\»’ (1.3.17) 


Introducing this into (1.3.16) we have: 

_ A 

-+l> 

»A 


A - d*p 


9ap 1?~ 


A A 

df_ d£ 
ds A d8 A 


= 0 . 


(1.3.18) 


We shall now make an assumption that is practically always 
fulfilled, that is: 




7 


p] = g^ [>v, /S] = |® v J- (1.3.19) 

That is, we shall assume the law of the tweedling of products 
quoted in “Notation” for the g 's and their derivatives. Under these 
conditions we can write (1.3.16) in the form: 


A 

Q a = 


d*P fTj dp dj 
d*\ + W ds A ds A 


= 0. 


(1.3.20) 
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A 

As in (1.3.9) we can introduce covariant differentiation D for 
tensors on the JL’th world-line: 


A 

JDT 0 = dT°+ 



T^df. 


(1.3.21) 


Using this notation we have 



4 A 

D d£_ 
ds A ds A 


= 0 . 


(1.3.22) 


Thus we have obtained the equations in a form very similar to 
those of a geodesic line. But it must he remembered that the last 
form was derived on the assumption that fpy> = <j>y is valid for 
the g’a and their derivatives. 


4. THE EINSTEIN PIELD EQUATIONS 

We shall recall here the famous Einstein equations for the 
gravitational field. They can be derived from a variational principle 
which we postulate in the form dW — 0 where W = W f +Wi and 

(1 -°> 

Here Te is the gravitational constant, the same that appears in the 
Newtonian equations: 

Te = 6.67-10~®cm 3 g -1 sec -2 . 

The integral has to be taken between space-like hypersurfaces. 
The function under the integral is the simplest scalar density that 
can be formed from the g^. Thus it is also the simplest Lagrangian 
that can be postulated if we ignore the possibility of a constant times 
V —g which leads to the additional “cosmological term”: a constant 
times g# in the field equations. 
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Indeed, under certain natural assumptions, V —gB is the only 
possible Lagrangian. This follows from the fact that the Lagrangian 
must be a scalar density and that we wish the gravitational field 
equations to contain derivatives of only up to the second order. 
This we assume, because the Newtonian equations for the gravi¬ 
tational field are of the second order. 

An invariant Lagrangian can be built only from the full Ric- 
mannian curvature tensor B^. Therefore, the only scalars to be 
taken into account besides B = g^B^, are those like B^Bf, 
B mia B ma , etc. The Lagrange function must be a scalar function 
of these invariants multiplied by a scalar density. Or, it could be 
a scalar density in the form det||jR oj5 ||. Generally, all the Lagrangians 
with the single exception of V —gB, lead to field equations of the 
fourth order. The constant in front of the integral in (1.4.1) is so 
chosen as to give the transition to Newton’s equations for e -*■ oo. 

If there is matter represented by point singularities, then the 
entire action should be a sum of W f from (1.4.1) and Wj from (1.3.1): 

W = TP/H-Wj. (1.4.2) 

A 

The constants m (0) in Wj play the role of coupling constants between 
the point particles and the metric field. 

We obtain Einstein’s equations for the gravitational field 
interacting with point particles by varying W with respect to &/)(*) 
and putting dW — 0. 

Let us calculate the variation of W, starting with W f i 

iw '“i 

<*I 
02 

= ^ / dxiV^B^dg^+Bdf^+^gg^dB.^) 

3 a 2 _ 

= f dx[V — giBtf— to af B)d!f f, +V-g<f l, dB al> l (1.4.3) 
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In the above equation we made use of the fact that 
&9 = 99°* ^9ap = ~99apSg ap . 


Because 


we have: 


where 




*W f = jdxV-g(R aP -^g aP JR)dg all +, 


J dxV—i 


q afi 8R„.. 


(1.4.4) 

(1.4.5) 


(1.4.6) 


(1.4.7) 


To find SS, let us remark that although j^J is not a tensor, 

its variation dj^j is. Indeed j“ v j A a datP is the change of the 
vector A p in a parallel displacement from x e to cP+daP. Therefore 
(^{“ v })JL a da^ is the difference between two vectors at the same 

point aP+doP. Both of them are obtained by a parallel displacement; 
one of them in the old field g^, the other in the varied field dg^. 
The difference between two vectors at the same point is a vector. 

Since A, and doP are arbitrary vectors, <5 j^j is a tensor. 

Using (0.13) we have: 

Using a locally Galilean coordinate system, that is, one such 
that at a given point 


9n,\ t = « ^ 

we have in place of (1.4.8) 


, = 0 and therefore j^j = 0, 
of (1.4.8): 


(1.4.0) 
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Because of the tensorial character of the <3 j^J we have in every 
coordinate system: 


where 


= * Q a ;a =.~L:(V-gdQ% 
y-g 

«r- fc} 


(1.4.10) 


(1.4.11) 


is a vector. Therefore, going hack to (1.4.7) we find: 

#2 

**=f <i-^> 


Ass umin g that dQ a vanishes on the space-like surfaces a y and cr a 
we have: 


<3 S = 0. 


(1.4.12) 


Therefore, using the definition (0.15) of the Einstein tensor, we can 
write (1.4.6) in the form:. 


dW f = 


J 


(1.4.13) 


Now let us vary the action TFj with respect to g^i 

L ^ITi ffj 




A—l ai 


1 d£_ 

2 


(1.4.14) 


where 


^9ap — J*a (4) (a> — {) dg^. 


( 1 . 4 . 15 ) 
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Substituting this into (1.4.14) we find: 


a d£ a A r A 

bWx = — Jjv m (0) c J - - dg 0 j dx <5 (4) (x — $) dg a ^(x) 


A -1 
<r 2 


i °* ^ dr 

= “a / dxd9a *>Zj m <°> c J (a? “^"SJ'• (1 * 4 - 16) 

A —l -oo 


Thus we have performed the variation of W. We demand that 
dW should vanish, if the dg^ are arbitrary, subject only to the 
condition dQ a = 0 on ff x and <r a for every choice of the space-like 
surfaces o 1 and <x 2 . Therefore, putting together the results expressed 
by (1.4.13) and (1.4.16), we obtain Einstein’s equations 


©<* = 


OTZK ~ 

~ z m <°> 

A-l 


J < 5 ( 4 ) 1 ®- 


—OO 


as A <*8 A 


^ (1.4.17) 

These equations have the form of field equations representing 
close-range action. Their left-hand sides consist of non-linear 
differential operations performed on the metrical field g^. Their 
right-hand sides consist of sources of the field formed by the moving 

particles. The masses have the character of coupling constants. 

Introducing the metrical energy-momentum tensor density for 
the moving point particles, we have: 


N oo A 

Z* 0 = JTto (0) c 2 J <5 {4) (a— £) 


A 

ar 


dP 


ds A ds A 


ds , 


(1.4.18) 


With its help we can write the Einstein, equations in the simpler 
and more general form: 


©<* = 


8 nk 

~T 


z*. 


(1.4.19) 
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We assume x° to be the parameter for all world-lines: 


A A A 

= !*(»»), £ 0 =®°. 


Then if we introduce 

c 2 ^ 

we can write (1.4.21) in the simple form: 


8 nJc 
d 


A** 1 


Therefore Einstein’s equations (1.4.19) become: 

n . 


a a. 


(1.4.22) 


2^°) d(x-§)i a lo P ]0 . (1.4.23) 


©“* = -8* ^(a®) 5 (x-1) |“ 10 £ V (1.4.24) 


These are the equations which will be especially convenient for 
finding the equations of motion of the third kind. 

One more remark:, if we put a® = at, then, as follows from 
(1.4.22), we have: 


A 

ti(i) = 



(1.4.25) 


Therefore, the development of fi in power series of o~ l starts 
with the order two with a constant that has a definite physical 

a , 

meaning, that is with 


5 . THE BIANCHI IDENTITIES 

Bor every the full Biomannian tensor satisfies the Bianchi 
identities: 


E^+R^+E^ = 0 . 


( 1 . 5 . 1 ) 
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This form is more general since it is supposed to be valid for 
any energy-momentum tensor. 

The left and right-hand sides of these equations are of entirely 
different character from one another. The left-hand side is Bin- 
stein’s tensor of the gravitational field. It has a distinct geomet¬ 
rical interpretation, being at the same time a metrical tensor. The 
right-hand side, however, is a purely physical tensor, representing 
the distribution of energy and momentum in space-time. 

Einstein always regarded this mixture of physics and geometry 
as fundamentally unsatisfactory, as a temporary device to be 
replaced in the future by a unified field theory, in which all physical 
fields would have their geometrical counterparts. To this problem 
Einstein devoted over thirty years of his life, looking for field equa¬ 
tions that would give solutions representing matter free of singular¬ 
ities. The result of this labour performed by him and his many 
collaborators seems disappointing to most physicists. 

Eor the present, therefore, we assume (1.4.19) and therefore 
accept a dualistic view; we assume the existence of matter which 
determines the geometry of the Riemannian continuum. Or, con¬ 
versely: from the knowledge of the geometry we find the distri¬ 
bution of momentum and energy. 

Let us go back to the case of point particles and the definition 
(1.4.18) of the energy-momentum tensor. This definition means 
that outside the world-lines we have: 


©^ = 0 


(1.4.20) 


that is, the field equations for empty space. Only on the world- 
lines does the field become singular.' This means: we assume the 
field equations for empty space and the existence of some world¬ 
lines on which the field equations for empty space break down. 

Eor further considerations it will be more convenient to use the 
energy-momentum tensor for point particles in a slightly different 
form. We can rewrite (1.4.18) by means of (0.22) in the form: 


8 Jifc 


’g# = 


BtiJc 

~ 


* A 
A-l 


d8 A 

(0) ~A' 

d£° 


A 

d£ a 


ds A ds 


d * <5(x-|(® 0 )). 


(1.4.21) 



42 


GRAVITATIONAL INTERACTION AND THE THEORY OF MOTION 


I.S 


Since they are tensorial identities it is sufficient to prove their 
validity in a special coordinate system. Let us take a locally Galilean 
coordinate system, that is, one for which, at a given point 

and = 0, j^j = 0, that is a coordinate system for which at 

a given point of the Christoff el symbols vanish and the g' s take 
the Galilean values. Thus at such a point and in such a coordinate 
system we can easily find the left-hand side of (1.5.1) from the 
definition (0.12) of the full Biemaunian curvature tensor. Assuming 
that the derivatives of the third order of the </’s are continuous we 
find that in the chosen coordinate system and at the point of the 
left-hand side of the last equation vanishes identically. Because 
of the tensorial character of this equation and because of 1 is an 
arbitrary point, the left-hand side must be zero always and every¬ 
where. 

Multiplying (1.5.1) by tf* and putting /x = a, we obtain what 
are also called the Bianchi identities: 

= 0. (1.5.2) 

When referring here or later to Bianchi identities we shall mean 
equations (1.6.2); they play an important role in G. B. T. and there¬ 
fore we shall prove them here independently of (1.5.1). They are 
identities, by which we mean that they are satisfied for arbitrary g’B 
as long as these have continuous third derivatives. 

One could ask: what is the source of the Bianchi identities? 
The answer is: they follow from the covariant character of Ein¬ 
stein’s tensor Q a? and from the existence of a variational principle 
leading to Q afi . 

We shall prove the identities in such a way as to exhibit the 
reason for their validity, which is: covariance with respect to all 
transformations. 

Once more we write the expression for W t \ 

3 a 2 

*”' cr 9 s - 


( 1 . 5 . 3 ) 
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This is obviously a scalar; therefore if we vary the coordinate 
system, then the corresponding variation of Wj must vanish. 

We shall assume an infinitesimally small change in the coordinate 
system co a -»• x a , given by: 

x a = x a - v a (1.5.4) 

■where is an infinitesimally small vector field. 

The general transformation law for the g 's is: 



. . , . dx' M da}' r 

(1.5.5) 

From this, 
we have: 

neglecting everything of order higher than the first, 


d'afifa) — 9ap(.®)-\- fyaP (®) 

(1.6.6) 

where 

= 00^1,3+ 

(1.5.7) 

or 

^9afi = Vatf “I - V(l\a • 

(1.5.8.) 

Similarly, we 

obtain 



<5 tf* = 

(1.5.9) 


Bemembering (1.4.6) and (1.4.12), we may write SWf in the form: 
<*2 

SWl = -^f 

crj a 

where 6Q a is given by (1.4.11). Let us now take a field rj a vanishing 
with its first and second, derivatives on the two hypersurfaces <r x 
and cr s . The surface integrals in (1.5.10) vanish and we obtain 

-is 

*1 


W — 


(1.5.11) 
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* w >- - 4 -/ -£iS (!•»•«) 

The field rj a being arbitrary in the region between a x and <r 8 , we 
must have 

* 

. ©“^ = 0 or ©/ ;(3 = 0. (1.5.18) 

These equations are, as our proof shows, independent of the 
field equations. The proof also shows the character of the Bianchi 
identities. They are a consequence of the invariance of the gravita¬ 
tional action, that is, they are a consequence of the tensorial 
character of the theory which is covariant with respect to all 
transformations. 


6. THE EQUATIONS OF MOTION ARE A CONSEQUENCE OF THE 

FIELD EQUATIONS 

Wo return to our fundamental problem. We shall show here 
that the equations of motion are a consequence of the field equa¬ 
tions. This is possible only because the field equations are non¬ 
linear (which does not mean, of course, that every non-linear theory 
gives the equations of motion). Indeed, in the Newtonian theory 
of gravitation, where the field equations are linear, as well as in 
Maxwell’s theory, the equation^ of motion are independent of the 
field equations. In Newtonian theory, we remember, we could find 
the field for arbitrary motion of the particles which were understood 
as sources of the field. The equations of motion which we obtained 
(in the Newtonian case) from the variational principle wore inde¬ 
pendent of the field equations; knowledge of them was unessential 
for solving the field equations. 

The situation is, however, entirely different in G. E. T. The 
field and motion are intrinsically connected with each other. The 
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equations of motion are a condition of integrability for the field 
equations. 

We shall now give the mathematical details' justifying these 
remarks and those of Section 1. 

Once more we write down Einstein’s field equations for a system 
of point particles: 


©^ = 


8jz1c 


N 


f 


<5(4) 


-4-1 


A dp dp 
- **£*:**• 


(1.4.17) 


We form the covariant derivative ;i j of the left-hand side. 
Because of the Bianchi identities we must have as a necessary con¬ 
dition for their integrability 


N oo 

-ts*. 1 8( *> (a> ' 

'A —1 —oo 


* dP dp 

■f )~r~ i7~ d8A 
ds A ds A 


L 


( 1 . 6 . 1 ) 


Or, briefly, and more generally: 

+ { “ } = 0 , ( 1 . 6 . 2 ) 

where' is the energy-momentum tensor forming the right-hand 
side of Einstein’s field equations. We are here, however, mostly 
interested in the case of point particles in which S“< 3 is defined by: 


N 


A A » 

-4 c ^ ci£ a d £p 

A**l —oo 

Thus in this case we have from (1.6.1) and (1.6.2): 


(1.4.18) 


~¥ ~ ^ m (o) f ds A (< 5 ( 4 ) (» - $ )\p 

.4 — 1 -oo ' " 


A A 
° 


dSj 


+ 


(1.6.3) 
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This expression may be put into a different form. We see that: 

. , {. 4’iP d , , {.iP *P 

f )i» *r A = - ~p Wx ~ f ui 


A' A 

d£ a d * d l8S a \ , 

A 

Using this, we may change (1.6.3) into the form 


(1.6.4) 


/ /(®) <3 ( 4 ) [co — i)ds 


since the integral of 


d idS 


(t>)‘ 


is equal to zero for every finite x. 


Furthermore, expressions of the type fd w vanish everywhere with 
the exception of the point at which a <4> does not vanish. Thus we 
may write: 

A £ A 

f{x) <5 (4) (a;-|) =fd (t) (x-£). 

We can therefore write (1.6.3) in the form: 

Or in the form: 

4®**- U™«>) f ^3 W (®-f)fl a y=0 (1.6.6) 


A—1 —oo 


where, as in (1.3.20): 


£a_4 2 r,T7\dr df 

- + W^ &;• 


(1.3.20) 
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We see that the integrability conditions (1.6.6) are fulfilled if 


A 

= 0. (1.6.7) 


But the converse theorem is also true. This means: from the vanish¬ 
ing of the divergence of the energy-momentum tensor, one can 
deduce the equations of motion. Indeed, let Aa s he a finite three- 
dimensional region of space-like character (that is a piece of the 
space-like hypersurface <x B ). We assume further that this piece is 

B B 

cut by the line f ° = f a (s B ) and only by this line; the point of 
intersection is, say, s B = s B . If = 0, then we also have 


s 


ds B 


a B~*B 


J" dx.n e Z a ^.p 



Because of (1.6.6) this means: 


( 1 . 6 . 8 ) 


B 

dr 

ds B 


/**•*»> / 

Aojg —oo 


ds B <5, 


3 °(*) 


B B 

(®-f )£“(**) = 0. (1.6.9) 


Taking into account the properties of the ^-function, the above 
equation means: 

b 3 

m {t) Q a [8 B ) = 0 . ( 1 . 6 . 10 ) 


As this equation must be fulfilled for every s B , we see that indeed 
the vanishing of the energy-momentum tensor implies the validity 
of the equations of motion. 

In the proof, we used the covariant notation. By giving up the 
covariant notation, we simplify the proof. We can write (1.6.6) 
using d instead of 
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Integrating over a three-dimensional region V around the 

B 

point f, we see immediately that 

B n 

1 d£° r 0 ' b b 

— -j- dxZ aP ;P = w (0) fi a (a? 0 ) = 0. (L#.12) 

&8j$ p. 

Thus, again we obtain the equations of motion as consistency 
conditions for the field equations. 

We have obtained the equations of motion in two ways: firstly, 
by varying W f +Wj (or really only Wj) with respect to the world¬ 
lines; secondly, by finding the integrability conditions for the field 
equations. What has been said up to now may be summarized in 
the following table in which <5.../<S... denotes the variational 
derivatives: 


d 8 Wt 

—j-(W/-f Wj) =— -j- = 0 Equations of motion 

8 

—j ( Wf+Wi) = 0 -> Eield equations (1.6.13) 

og 


Thus we see the essential difference between Newton’s and 
Einstein’s theories of gravitation. To find the world-lines in G-. B. T. 
we have to know the field. But to find the field we have to know 
the motion. It cannot be arbitrary. 

What is the way out of this difficulty? We shall see that there 
are two ways out. The one mentioned in Section 1 is useful in 
principle rather than in practice. It consists of changing Einstein’s 
equations so that the motion is arbitrary. Then, imposing proper 
restrictions upon the motion, we obtain the solution of the equa¬ 
tions. The other, which we mentioned in Section 1, consists of 
devising an adequate approximation procedure. We shall have 
much more to say about each of these methods later. 
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I. 7 


7. FORMS FOR THE EQUATIONS OF MOTION 
OF THE SECOND KIND 


As we have seen, in the previous section, equations of motion 
of the second kind are closely connected with the Bianchi identities. 
These identities follow from the invariance of gravitational action 
with respect to infinitesimal variation; this we have proved with 
the help of formula (1.5.13). 

Now, making use of the field equations in (1.5.11), we obtain: 


fdx^dg^ =0 or =0, (1.7.1) 

where according to (1.5.8) and (1.6.9): 

-dg^ = (1.7.2a) 

dgafl = + 9afiU< if- (1.7.2b) 

If we introduce the last values into (1.7.1) we see that because 77 
is arbitrary, (1.7.1) is equivalent to 

%% = 0 or = 0. (1.7.3) 


Thus, instead of (1.7.3) we can equally well regard equations 
(1.7.1) as our equations of motion, just because of the arbitrary 
character of the jj ’ s . 

Let us introduce into the first of the equations (1.7.1) for 

the known expression (1.4.23): 

, N 

Jc ^ o VI 4, 

-^£<0 = 2;iU<5£VV (1-7.4) 

A -1 


Therefore (1.7.1) becomes: 



A-l s'O 


(1.7.5) 


Since the ^’s are continuous functions of the a>’s only (not 
depending on the |’s), we have: 

A A A A A A A 


&9al> = g att V^+g^^la+ga^ if- 


(1.7.6) 


Motion and relativity 


4 
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Since 


^ * drf 

Mo - (1-7.7) 

A 

we have, because of the arbitrary character of the ?)’s which, as we 
assume, vanish at the end of the time interval: 

/7 A r—' A A ^~ A A 

“T-g-(Ma/J ^ 10 ) — i'i w fl , /fl<!af ,, |of V |0 — 0, A = 1,..., N. (1.7.8) 
ax 


Let us first show that these equations are entirely equivalent 
to those gained by the use of the variational principle, that is to 
(1.3.16) 


A A A 

d /A< i&\ 1 J^de d r 

ds A \ 9al, dsJ 2 9Ha ds A ds A 


= 0 . 


(1.7.9) 


Before we prove this equivalence, let us make one remark: 
these equations are typical for the two formalisms which are used 
throughout our work: the s formalism, if « is an independent 
parameter, as in the last equation, and the t formalism, if t is the 
independent variable, as in (1.7.8). The s formalism connected with 
the use of 8(^(x) has the virtue of revealing immediately the covar¬ 
iant character of these equations; the t formalism, connected with 
the use of <3(x), does not have this virtue. It is, however, easier to 
handle in actual calculations. 

The proof of the equivalence of the last two equations is simple 
enough. Indeed, introducing into (1.7.9) instead of ds A : 


A 



dx° 

~A T ’ 

M®°) 


(1.7.10) 


(which is identical with (1.4.22)), we easily find that (1.7.9) becomes 
(1.7.8). 
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The reverse statement is also true, but the proof of it is more 
subtle. It would seem that by introducing into (1.7.8) the new 
parameter s A defined through the last equation we do obtain (1.7.9). 
This is true, but the remark does not make the proof complete. 
The question arises: is the ds A introduced through the last equation 
identical with the ds A introduced through 


A A 

r->d£ a as* 

9aP d8 A dS A 


= 1? 


(1.7.11) 


We shall prove that it is. In fact, differentiating the last equation 
with respect to s A , we obtain according to (0.18) 


d /A <lj a dp \ 

ds A v afi ds A ds A J 


AAA 

J- as* ar a i* 

9af,lt ds A ds A ds A 


+ 


A 


A A 

d^ a 2 £ a 


ds A ds A 2 


= 0. 


(1.7.12) 


Obviously, as we mentioned before, the substitution of the 
new parameter ds A defined by (1.7.10) turns the form .of equa¬ 
tion (1.7.8) into the form (1.7.9). But in it, ds A is the “proper time” 
since it satisfies (1.7.12). This can be seen in the usual way by 

A 

multiplying (1.7.9) by d£ a jds A . In this way we obtain an equation 
entirely equivalent to the last one, from which it follows that 

A A 

A dS a at* A 

< '“'&7*I = oonst (1 - os) 


is an integral of the differential equations for this world-line. 

A 

In the definition (1.7.10) of ds A we have used the constant 
Obviously it can always be so chosen that the constant in the last 
equation equals one. This completes our proof. 
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Both these equations of motion can obviously be deduced from 
a Lagrangian: the first one,- (1.7.8), from a variational principle 


<5 J Ldx° = 0, 

L — £ f g aP /tf a ,0 P\o <5(x- IW; 


(1.7.14) 


the second one, (1.7.9), from 


8 (~ 2™'( 6 ) c f d * A ) = °» 

AAA 

ds A = J g afi d£ a d£ / ’d (4) d%. 


(1.7.15) 


In both cases, the g’s should not be varied with respect to the f’s. 

In exactly the same way in which we passed from (1.3.16) to 
(1.3.20), that is by raising the indices, we can obtain from our two 
equations, again assuming the rule for tweedling the products: 


A A A f _ 1 A A 


(1.7.16) 


and, as in (1.3.20), we obtain: 


-+H— 

si W ds A 


A A 

at" d£ p 


(1.7.17) 


Hote that equation (1.7.8) could also be obtained from 

$4, = 0 (1.7.18) 

by putting into this equation the expression corresponding to that 
in (1.7.4) and following an argument similar to the one employed 
in the previous section where the equation 


= 0 


(1.7.19) 


was employed. But obviously we wish to get the same equations 
of motion whether we use the lower or the mixed indices in the 



I. T FORMS FOR THE EQUATIONS OF MOTION OF THE SECOND HIND 53 

field equations, from which we deduce either (1.7.18) or (1.7.19) 
respectively. These two forms of equations are equivalent only if 

A A A 

(1.7.20) 

A A A 

<n>, q! = ? I>, el- . d-7.21) 

Thus only if the. tweedling laws for products are assumed, are the 
field equations ot G. B. T. self-consistent. 

We shall introduce one more form of the equations of motion 
which is of some practical value, namely: in equations (1.7.16) 

A A 

there are 427 unknown and /x. From them we can eliminate the 

A 

/S s and obtain 327 equations for the 327 unknown f fc . 

Putting a ess 0 and a = a in (1.7.16), we have: 


A 


t»+i {£} = 

(1.7.22) 

A A 

A 

AA A ( -1 A A 

(1.7.23) 


A 

Eliminating / x ]0 from the second equation with the help 

of the 

first, we have the three equations for the .4’th partiole: 


A 

A 


^ °i° # + ({“„ 


(1.7.24) 

This we can write more explicitly though less symmetrically: 

A A 

A A 


^°|oo + {o a o}~{oo} f °l® +2 {^} £ n |« -2 {o»} £ a |o£*l<> + 


A 

A 


+{™} 

}fm mi / ^ l jt a j t n A 

10 ff 10 — $ 10S’ |0f |0 — u - 

(1.7.26) 
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A 

Thus the mass p(x°) can he made to disappear from the equa¬ 
tions in which ®° is used as a parameter. This fact suggests some 
remarks of a more general nature. 


In Newtonian Theory, the equations of motion of the first and 
second kind had one characteristic feature: it was possible to cross 
out from both sides of the equations the mass of the particle whose 
motion was considered. The mass on the left-hand side, appearing 
as the coefficient of acceleration, was the inertial mass, whereas 
the mass on the right-hand side, which was always multiplied by 
the mass of another particle, was the gravitational mass. This fact 
was the mathematical expression of a truth known since the times 
of Galileo: that all bodies fall in empty space with the same 
acceleration. This was confirmed again by the more precise exper¬ 
iments of Boland Eotvos which showed that the inertial mass and 
the gravitational mass are numerically equal. 


In our presentation of Newtonian Theory, this physical truth 
found expression in Section 6. In writing the Newtonian action Wj 

A 

we used the same coefficients m for the inertial constants as for 
the coupling constants with the gravitational field. That is: in 


A • • A A 

2 imps* we used the inertial constants and in — dx8{x—$)q>, 

A 2l 


the gravitational constants. There was no a priori reason to 
regard the w’s as identical in these expressions. We wore forced 
to do so only because of the experiments of Galileo and EotvOs. 
But we can well imagine a world in which the inertial and grav¬ 
itational masses are different from one another. 


In G. B. T., however, the equivalence of gravitational and 
inertial mass is embodied naturally into the foundations of the 
theory. Indeed, unlike Newtonian Theory, we have in G. B. T.: 


W x = 


N 



A 

m (0) c 


; 


n 


A A A 

Lords ' 1 ) 11 
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and there is no place here for the distinction between the inertial 
part and that representing the interaction with the field. Thus, 
there is no place for a distinction between inertial and gravitational 
mass. As we shall see in Chapter II, in the first approximation, with 
respect to 1 jc, the expression W x in G. R. T. consists of two parts: 
one of them representing the inertial part, the second the inter¬ 
action. But in both parts we do have identical masses, not because 
of a special assumption, but as a consequence of G. B. T. 

One of the results of the equality between inertial and grav¬ 
itational mass is that in the equations of motion of the second kind 

A 

the mass m (0) does not appear at all. What does this mean physic¬ 
ally? It means that the acceleration at a given point depends only 
on the field at the point where this particle is situated, and not on 
its mass. The fact that the equations of motion of the second kind 

A 

are externally independent of m (0) is the mathematical expression 
of the equivalence principle in G. R. T. 

The same is of course true for equations of motion of the first 
kind, that is for the equations of motion for a test particle. Then, 
of course, the law of the tweedling of products is fulfilled and the 
entire tweedling procedure becomes a trivial matter of putting the 
f’s in place of the ®’s. Thus our equations of motion of the second 
kind, changing into equations of motion of the first kind, rigorously 
become those of a geodesic line. They are entirely independent of 
the mass of the particle, not only externally; the field itself is free 

A 

from any dependence on Am, as long as it is small enough not to 
disturb the external field. 

In the natural explanation of the equality between gravitational 
and inertial mass, we see the great and essential superiority of 
G. R. T. over other gravitational theories. 

The difference between the gravitational and inertial mass 
appears only if we consider equations of the third kind. More about 
this in the last chapter. 
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8. THE EQUATIONS OF MOTION IN GRAVITATIONAL AND 
NON-GRAVITATIONAL FIELDS 


Up to now we have considered gravitational equations with 
the energy-momentum tensor 


A -1 —oo 


We sha D now briefly consider the more general case, in which 
some fields are added, say the electromagnetic or the mosonic field, 
or both. Let us assume, only for the sake of brevity, that one 
vectorial field is added, denoted by A a . (Of course, everything that 
is said here remains true for tensor or spinor fields.) Then the 
gravitational equations change their form. Besides the energy- 
momentum tensor belonging to the particles, there also appears 
a tensor belo nging to the field. If we call this additional tensor 
and the entire tensor jE^, then Einstein’s equations take the form: 


©^ 


8tcTc 


(£<*+©<*)__CK 


( 1 . 8 . 2 ) 


We do not need to care about the actual form of It is sufficient 
for. us to know that is a tensor formed of the A a , perhaps of 
their covariant derivatives and it may also depend explicitly on the 
metric tensor and on the velocity vectors. Thus, symbolically, we 
may write: 

B* = (1.8.3) 

The metric tensor and the equations of motion are determined 
by (1.8.2), if the external field A a is known. As to the JL’s, they must 
be determined by further field equations. Since wo do not care 
here about their specific form, we may write them symbolically 
in the following way: 

O a {A, g, I') = 0, (1.8.4) 


where the 0„ are some operations to be performed on the A, g 
fields and on the velocity vectors. Our only assumption concerning 
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the operators O a is that equations (1.8.4) should be covariant 
with respect to arbitrary transformations. Usually both equa¬ 
tions (1.8.2) for the gravitational field and (1.8.4) for the non- 
gravitational field can be deduced from a Lagrangian, the first 
by a variation of the metric tensor, and the second by a variation 
of the external field A. 

Obviously, because of the Bianchi identities wo have: 

= 0. (1.8.5) 

Now there are two possibilities: 

1 . 6 ^ = 0 , 2 . # 0 . 

The first possibility may be a result of the field equations for 
the A’s. To have a definite example before our eyes let us think 
about the electromagnetic tensor without sources. Then G°^ = 0 
because of Maxwell’s equations. In this case we have, for the equa¬ 
tions of motion: 

/$°V X = °> A = 1,2,..., A, (1.8.6) 

a 

exactly as before. In this case the equations of motion are influ¬ 
enced only by the change of the metric field, a change caused by 
the existence of the A field. But in their structure the equations 
remain unchanged; they are the equations of the “geodesic line”, 


where the g^, are replaced by g^, g a ^ e . 

The second possibility gives us different equations of motion. 
They are now: 


<£“*5 = $**+©*, = 0 . 


(1.8.7) 


An interesting case arises when the vector density is 
a linear function of the d’s: 


% = c‘£i- 


( 1 . 8 . 8 ) 
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This is so, for example, in the case of electromagnetic fields in which 

A A 

the particles of mass w (0) have a charge e. In such a case, character¬ 
ized generally by (1.8.8), we have for the equations of motion: 


m. 


(«) 


jT\dPiP A, 


A = 1,N. (1.8.9) 


These are the general equations of motion if (1.8.8) is assumed. 
In particular, in the case of an electromagnetic field, we have the 
Lorentz equations of motion. 

Usually, (1.8.8) will be valid in the case of point particles, as 
a consequence of the field equations (1.8.4). Thus in such a case 
we may say that the proper equations of motion are a consequence 
of the field equations (for the JL’s) and the gravitational equations. 
But of course the same equations of motion are valid in an arbitrary 
coordinate system, even if the gravitational field vanishes, that is 
if the gravitational constant h -*■ 0 and R^ = 0. In this case 
the equations of motion (1.8.9) are usually assumed, because the 
field equations for the A’s usually become linear and, as wo know, 
a linear theory cannot give the proper equations of motion. They 
have to be added to the field equations. We see the advantage of 
introducing the gravitational field. Even when it goes to zero, a trace 
remains in the form of the equations of motion/ Otherwise such 
equations would have to be assumed separately. 

It oan easily be shown that the K a vector in (1.8.8) must be 
a space-like vector. This can be seen by multiplying the last equa¬ 
tion by 


Safi— ( 1 - 8 - 10 ) 

di A 

The left-hand side of the resulting equation must vanish. This 
is so because of the definition of ds A and because of the considera¬ 
tions in Section 7. Therefore the right-hand side must vanish too. 
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Thus we have: 


A 



A 

dr 

d%A 


E? = 0 


( 1 . 8 . 11 ) 


which proves our point, since the velocity is a time-like vector. 


8. EQUATIONS OP MOTION IN DIFFERENT COORDINATE SYSTEMS 

We shall answer the following question: what happens to the 
field equations and to the equations of motion if we pass from one 
coordinate system to another? As in the last section, we assume 
the existence of the gravitational field and of the physical field 
characterized, as before, for the sake of brevity, by one vector field 
A a . How we pass from a given coordinate system to a starred 
coordinate system, that is from 



(1.9.1) 

to the expressions 


x* a , g*<*(x*),A* a {x*), f*«, 

(1.9.2) 

where 



<v a = a?°(a>*), x* a 

= »*“(»)> 



dx* a 

(1.9.3) 

A - A ° dx° ' g ~ r M dx a ’ etC - 

Similarly the expressions 


f *\a« 

i for 

(1.9.4) 

go over into J 


(1.9.5) 


where the starred Ohristoffel symbols have the obvious meaning. 
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But the field equations of General Relativity Theory cannot 
depend on the choice of the coordinate system. This means that, 
■whereas in the old coordinate system we had field equations that 
could be written symbolically in the form: 

9 W, n, 

c 

O a (x)[A{iD),g(x), I'] = 0, 

in the new coordinate system the equations can be written in the 
form: 

n, 

WlAMrt-o. 

All this is almost trivial f yet the fact that in the new coordinate 
system the the E* and the operator O a depend on the new field 
components and act in the same way as they depended and aoted 
before in the old coordinate system on the old field components, 
is not always well enough understood. 

From this it follows that the equations of motion in the new 
coordinate system are 

E* (A* (x *), g* (»*), £*')* = 0 (1.9.8) 

where the star over the semicolon means, as in (1.9.5), that the 
Ohristoffel symbols have to be formed from the g*{x*). 

Again the result is almost trivial: the equations of motion are 
covariant with respect to any change in the coordinate system. 

Since the field equations and those of motion have the same 
form in all coordinate systems, then the question arises: are all 
these coordinate systems permissible? But to describe motion in 
any particular instance we must specify the coordinate system 
in some way. Yet by doing so we violate the spirit of covariance, 
so important in G. R. T. 

In principle there are three ways out of this difficulty: 


(1.9.7) 


(1.9.6) 
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1. The coordinate system is specified by certain added coord¬ 
inate conditions. Indeed, in general four coordinate conditions can 
be added to the field equations specifying the coordinate system. 
These in turn specify the four arbitrary functions appearing in the 
solutions of the field equations. 

2. No coordinate conditions are added. The general solutions 
contain four arbitrary functions. If necessary, we may specify them, 
but we are not restricted in our choice by any added coordinate 
conditions. 

3. All coordinate systems are permissible if they satisfy one 
condition: the' world is Euclidean and the coordinate system becomes 
Cartesian at infinity. 

As to the first choice: it is represented by Fock and his school. 
The coordinate condition added to the field equations is the 
“harmonic” one: 

(V^% = 0. tl.9.9) 

Obviously such a condition, if used exclusively, changes the spirit 
of G. R. T. since it distinguishes a group of coordinate systems. 
Indeed the harmonic coordinate condition gives especially simple 
solutions in the case of linearized gravitational equations in empty 
space. However, such a case of solutions of linearized equations 
seems to be of little physical interest. Otherwise the harmonic 
coordinate system is rarely used. For example, it is not used in the 
classical Schwarzschild solution, and, as we shall show in the next 
chapter, it has little to do with the problem of motion. Therefore, 
from the practical point of view its advantage is very small, if any. 

Yet from the theoretical point of view, if the harmonic condition 
should always accompany the field equations, this seems to us to 
be a step backwards from the ideas of G. R. T. The ideas of G. R. T. 
were, if not deduced, then made plausible by the simple and very 
suggestive ideal experiment concerning the falling lift. True enough, 
this showed the equivalence of two systems based on the equality 
of inertial and gravitational mass, only in a restricted time-space 
region. Yet it was important for the development of G. R. T. that 
two systems, though not moving uniformly relatively to each other, 
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are equivalent in a certain sense. The observer outside the lilt can 
say that the lift moves with accelerated motion in his gravitational 
field. The observer inside the lift may claim that his system is an 
inertial system towards which the earth moves with accelerated 
motion. All this reasoning becomes meaningless if we admit only 
harmonic coordinate systems. 

How, about the second possibility. We cannot describe motion 
without specifying a coordinate system. We must have one in order 
to check the results of our calculations with experience. But to 
bring in any special coordinate system is a violation of the ideas 
of G. B. T. 

What remains, therefore, is the third possibility: all coordinate 
systems are allowed, as long as the field becomes Euclidian at 
infinity and the coordinate system Cartesian. To make our point 
clear, let us consider the simple case of the planetary motion around 
a very large sun. Let us take a plane at infinity, or so far from the 
sun that we can neglect its gravitational field. There, at infinity, 
there is sense in talking about planes, perfect clocks, and rigid rods. 
Every event in the world is mirrored uniquely on the plane by 
a light-ray emitted from this event and reaching the plane perpend¬ 
icularly to it. Therefore, to every event there correspond four numbers 
T, X k in the infinitely removed plane. Thus the motion of a planet 
has its image on such a plane. IE the motion is periodic, then its 
image forms a closed line. Now we can turn the plane until the 
area circumscribed by the image of the moving planet becomes 
a m axi mu m. Such mirror motion, in such an infinitely removed 
plane, is objective in this sense: for its description no further speci¬ 
fication of the coordinate system is necessary. (We may regard 
such mirror motion as objective only if the original motion was 
also in a “plane”; that is if the observer at infinity does not find 
any Doppler effect.) 

Thus we see from our example, that, without referring to any 
particular coordinate system we may, at least in principle, describe 
the motion objectively. If we call the coordinates on the “parallel” 
infinite plane T, X k then we have 

X k = X k (T) 


(1.9.10) 
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as the objective description of the motion. In an arbitrary coordinate 


system we have: 


1u 

II 

%* 

(1.9.11) 

where, theoretically at least, we have 


T = t -j- oo 

(1.9.12) 


since the light-rays take infinite time to reach the infinitely 
distant plane. 

Now we can define a convenient coordinate system in which 
we describe the motion (1.9.11) mirrored on the infinitely distant 
plane. It will be such that the motion on it is almost objective; 
that is, the difference between and X k and between AT and At 
is below the experimental error. We shall return to this problem 
later, perfecting mathematically the ideas described here only 
generally. 


10. ON THE METHOD OP SOLUTION OP THE FIELD EQUATIONS 
WITH THE HELP OP THE DIPOLE PROCEDURE 


If we wish to solve the problem of motion for point particles 
and find the gravitational field, that is, if we wish to find f°(l) 
and g a p{a:), then we must deal with the following system of equations: 


y? a f a tip dp 

2 j m («) J <^< 3 ( 4 ) («-•£)■ 


A-l 


Us A d*A ’ 


£T(Sa) = 


d 2 P jyj dp d£ _ 

&8a Wj ds A d8A 1 


tojL = (g^dptip) 1 ' 2 . 


( 1 . 10 . 1 ) 

( 1 . 10 . 2 ) 

(1.10.3) 


Equations (1.10.1) are those of the field; they determine g^ 
as functions of the world-points. 
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A 

Equations (1.10.2) are those of the motion; they determine £“ 

A 

as functions of s A , or £ k as functions of a>°, if -we introduce into them 
the corresponding change of the parameter. 

But we also know that the two systems of equations are not 
independent of each other. Because of the Bianchi identities = 0 
we have, as was shown before: 


ds A d w (x-1)^A). (1.10.4) 


N co 

( Vto (0 ) I d8 A 6 w (<c-() 

Mol — oo 


a dr de* 

d8 A ds A 


N oo 

= ^ ™( 0 ) J 


A =1 


Therefore, we repeat: equations (1.10.2) are a consequence of the 
field equations forming the condition for their integrability. 

Thus, we cannot solve the field equations with arbitrary motion. 
That is: we cannot find the g 's as functionals of arbitrary motion 
of the sources appearing on the right-hand side of the field equa¬ 
tions (1.10.1). Such a solution will exist only if the motion satisfies 
the equations of motion, that is (1.10.2). 

Therefore, the usual procedure of solving the field equations by 
assuming arbitrary motion, a useful procedure in all linear field 
theories, is of no avail here. 

How can we avoid this difficulty 1 ? We have said before that 
there are two ways of doing so; firstly: the dipole procedure, which 
is of principial rather than of practical value, and which we shall 
describe presently; secondly: the approximation procedure with 
which we shall deal in Chapters II and III. 

Modifying Einstein’s equations we assume that besides the 
moving point particles the sources of the field consist of an addi¬ 
tional, auxiliary field S)°. This field has to be chosen in such a way 
as to make the motion of particles arbitrary. Thus we shall be ablo 
to solve the non-Einstein equations in which the artificial field 
appears devoid of any physical meaning; we shall be able to solve 
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these generalized non-Einstein equations for arbitrary motion. 
But in the case in which the motion is the right relativistic motion, 
then the additional 2)° field vanishes. Or, conversely: if the additional 
field 2)“ vanishes, then the equations become the Einstein equations 
and the motion becomes the right motion, satisfying equations 
( 1 . 10 . 2 ). 

We call the additional, artificial field ©“ the field of dipole 
potentials. 

Let us, therefore, consider the field equations of General Rela¬ 
tivity Theory with some additional tensor 

©* = - ^ ($<*+©<*) = - ^ «*, (1.10.5) 

0 0 


where, we assume, as in Section 8: 

& = f ds A d { *')• (1.10.6) 

A —00 ^ 

Then, the equations of motion are: 


= ^ iP + = 0. (1.10.7) 

We ask: is it possible to find the <S a|3 and the field equations for 
the D vector in such a way that the equations of motion would 
be fulfilled automatically for any arbitrary motion? Indeed, this 
is possible only if becomes zero as a consequence of the field 
equations for 2)“. This is precisely the case if the field equations for 
2)“ are: 


©"*,-28*©,- - 2m (0) c* J ds A d w {w-f)S2°(8 A ), 


— 00 

A,., dr 

aM = "iA + M 


ds A ds A 


and 


( 1 . 10 . 8 ) 


(1.10.9 


Motion and relativity 


+ 2> fta — f* 25%. 


5 
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Indeed, because of Section 6, we can rewrite equation (1.10.8) 
in the form: 

X a *j-R a( “&t l = ( 1 . 10 . 10 ) 

According to the known rules of tensor calculus, we hare 

= 2 )^+ 3 )%-®*/+^ 

= = 0 , ( 1 . 10 . 11 ) 

which is identical with the field equations (1.10.10). 

What are the characteristic features of this theory? 

Firstly: in this theory the motion of particles is arbitrary. 
Indeed, using Bianchi identities in (1.10.5), we obtain the field 
equations for the additional field ®°, which are satisfied anyhow. 

Thns the system of equations (1.10.5-8) has a solution tor 
arbitrary motion and we can find from it: 

( 1 . 10 . 12 ) 

= ®"(«) cVe, (1.10.13) 

Here the square brackets, as usual, denote functional dependence 
Secondly: let us put ®° = 0 into equations (1.10.5-8); then 
those generalized field equations become Einstein’s equations, and 
the equations of motion (1.10.2) must be satisfied. Thus, at the 
end of our procedure we may return to Einstein’s field equations, 
fi nd i ng the motion from the condition that the additional dipole 
field vanishes. 

Thirdly: let us assume that the field equations £2* = 0 are 
fulfilled; that means that we hare chosen the motion, which is 
arbitrary in the generalized equations, so that it is the right motion 
according to G. E. T. Does it follow from this choice of the motion 
that the field ®“ vanishes? This is not necessarily so, because the 
field equations for ®° are then: 

= 0 . ( 1 . 10 . 8 ) 

These homogeneous equations may have solutions different from 
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zero. Only if we also assume that as a solution, of 35“ for the case 

A 

13“ = 0, we take the solution 2)“ == 0, do we have: 

= 2> a (a>) [f, I, ...,£] = 0. (1.10.14) 

Let us assume that we have solved the generalized system of 

A 

equations for g a ^{x). Let us put these g ’s into 13“ (s^). If we wish, 
since the motion is arbitrary, we may assume the equations of 
motion 

i3° = 0. (1.10.15) 

Since in this case we have agreed to admit only the solution 3)“ = 0, 
then the field equations become the Einstein equations. 

Thus we have formulated the general theory, according to which 
we can find the equations of motion of the third kind. There will 
•be more on this subject, in the next Chapters. 

Why do we call this additional field the field of dipole poten¬ 
tial i The answer is simple. Looking at the field equations (1.10.8) 
for 35“, we see that its solution will have the form: 

vi 13“ 

®°~2j 1-If' (1.10.16) 

Therefore, in the additional field appearing in (1.10.5) for g# we 
have derivatives of 25° — expressions which behave as though they 
represented “gravitational dipoles”. But gravitational dipoles do 
not exist in reality. By getting rid of this artificial field we return 
to Einstein’s equations and the proper equations of motion. 
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THE APPBOXIMATION METHOD AND THE EQUATIONS 

OF MOTION 

1. GENERAL REMARKS ON THE APPROXIMATION METHOD 

Id. the first chapter we gave the rigorous theory of the equations 
of motion in G. B. T.; we discussed the connection between these 
equations and those of the gravitational field. But what we said 
on this subject was rather formal. We did not consider any physical 
problems that could be compared with observation, at least in 
principle. 

Because of the non-linearity of the equations we are not able, 
for example, to give a rigorous solution of the two-body problem 
which, in Newtonian theory, is as easy as the one-body problem. 
We cannot find the field; therefore we cannot find the. motion. 
All that we can do is to use an approximation method which will 
allow us to find the field and the motion approximately. Such 
a method is more than adequate from a practical point of view. 
It will lead us to new effects which can barely be confirmed even 
by the most precise observations. 

In the case of non-linear equations the only method of approxi¬ 
mation used is that of the small parameter. Boughly speaking, the 
method consists in developing all quantities which appear in these 
equations in power series of this small parameter, then solving 
successively the equations formed by the coefficients of those 
series. 

What quantity should we choose as the small parameter in 
G. E. T. ? It would seem that a good choice would be the gravita- 
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tional constant k which plays the role of a coupling constant 
between the geometry described by the Einstein tensor €T and 
the physics described by the energy-momentum tensor I 1 ”. How¬ 
ever, this is not so. We demand that our parameter give us New¬ 
tonian motion in the first approximation, post-Newtonian motion 
in the next. However, using A; as a parameter, we obtain uniform 
motion in the first approximation, and in the next — Newtonian 
motion plus forces depending on the velocity of light. Such a method 
is certainly unsuited to the description of planetary motion, where 
the velocities are much smaller than the velocity of light. 

Thus another possibility, the one that we shall choose here, is 
to take as the small parameter 

A = —; (2.1.1 

e 

c = velocity of light. 

This method turns out to be very effective for the investigation 
of the problem of motion in G. R. T. 

Up to the present day the Newtonian theory, in which the velocity 
of light does not appear at all, forms the basis of astronomy. In 
comparison with this theory, G. R. T. is much more complicated 
from the technical point of view. Yet from the principal point of 
view, it is simpler. It does not require the strange concepts of action 
at a distance, or an inertial coordinate system, and the equality 
of inertial and gravitational mass does not appear in it as an artificial 
assumption. Yet the practical usefulness of G. R. T. demands that 
the following two conditions be fulfilled: 

1. Newtonian theory follows from G. R. T. as a special limiting 
case. 

2. G. R. T. can, without introducing any arbitrary constants, 
predict new events which do not agree with Newtonian theory. 

From what limiting case of G. R. T. can we obtain Newtonian 
theory? The answer is obvious. Newtonian theory, which is not 
a proper field theory and in which the constant c does not appear, 
can only follow from G. R. T. in the limiting case oo; or in 
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other words, “action at a distance” means a propagation with 
infinite speed. 

Thns we expect that for 

A = — -> 0 (2.1.2) 

o 

the equations of G. E. T. become those of Newtonian theory. In 
practice, we can regard e as infinitely great if v[e (where v is the 
characteristic speed of the body) is much smaller than one. On the 
other hand, if we do not neglect v jo, we may be led to new phenomena 
outside-the frame of Newtonian theory. 

We believe that these remarks justify the choice of A = 1/c 
as the parameter of our approximation procedure. Thus we assume 
that all our quantities are analytical in A. That is, if / is such a 
quantity, we shall write: 

/=/ + /+/ + /+••• (2.1.3) 

0 12 3 

We denote the order of each of these quantities by a subscript, that is: 

oo 

f=Y* n fn, / = r/„. (2.1.4) 

»=• n 

Besides the obvious rules (e. g. that the product <py is of the order 

k i 

k+l) there is one which is less obvious but essential to our later 
argument. It is connected with the differentiation of / with respect 
to time t, where 

a ? 0 

t = —r = Aa?°. (2.1.5) 

c 

We assume t to be of the order zero in A. Therefore: 

, a/ _ 1 df df 

0 7| ° da>° c dt dt 

must be of the order one. Or generally, we could write: 


( 2 . 1 . 6 ) 
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la words: differentiation of / with respect to x° raises the order 
of / by one; we shall omit the subscript 1 under a zero denoting 
differentiation with respect to x° and sha l l write simply 

St 

Differentiation with respect to co k does not change the order of /. 
This means that all quantities are assumed to change slowly with ®°. 


2. ON THE DEVELOPMENT OF THE METRIC FIELD 

We shall deal here solely with isolated material systems. Far 
from any matter, we shall have a vanishing curvature of the time- 
space continuum. Thus we can assume a flat time-space at infinity 
in which we can still fix the coordinate system in an arbitrary way. 
However, for the sate of simplicity, we shall assume a Cartesian 
coordinate system there. Thus, as far as measurements are concerned, 
there are the same conditions at infinity as in an inertial coordinate 
system, and as far as the metric is concerned we have for = ®*®*: 

= Vafii iVoo — lj Vok — b, rj ik = —da,). (2.2.1) 

r-*oo 

The transition from the Riemannian metric g aP to the Cartesian 
or Galilean metric rj^ at infinity plays an essential role in the 
theory of measurement as previously outlined (I. Section 9). We shall 
return to this subject later (V. Section 1). 

We shall now write the Riemannian metric in the form: 

= ‘‘lafi + ( 2 . 2 . 2 ) 

where the h a p go to zero as the distance r = (a® x*) 112 goes to infinity. 
How we wish to develop h aP into a power series in X. Since for c oo 
we have ->• 0, wo must start with at least the order one, that is 


&O0 = ha0~\~ hap 


( 2 . 2 . 3 ) 
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But by examining the structure of the field equations of G. R. T. 
more deeply, we can say more. We can show that 

^oo = ^om = 0. (2.2.4) 

i i 

(Later we shall learn to restrict still further). 

Indeed, we can show the validity of the last equation by 
comparing the equations of motion for a test particle in Newtonian 
theory with those in G. R. T. In the Newtonian theory we have 
for a test particle of mass Am under the influence of a gravitational 
field <p: 

h • it 8 

d(-Am) f dt(o 2 -ipf + ip) = 0, f (2.2.5a) 

h 

The equation of motion for a test particle in G. R. T. is: 
r ds 

( -Am) sj —dt = 0, ds 2 = ( Vafi + h afs ) ardt*. (2.2.5b) 

Developing, the square root with respect to A we obtain from the last 
equation: 

h . . 

(-d™)d/ at{c*~ i®+ wh 00 + ch^e+ r h - o. (2.2.«) 

<1 

Therefore: 

fc 2 h m -+ cp or like A 0 ; 

->at least like 1 [c, or like A; 

\nn ^ n £ n -»■ at least like 1/c, or like A. 

The last correspondence proves the validity of (2.2.4). Therefore 
we have: 
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Aoo — ho 4“ ^oo 4 ^oo 4 • • • ? 

2 3 4 

hn == ^'0n 4 4 ^oti 4- • • > 

2 3 4 

= 4 * hmn 4 

12 3 


(2.2.7) 


as a condition for the consistency of Newtonian Theory and G. E. T. 


3. ON THE THREE LEVELS OP OUR REASONING 


In principle, we picture material bodies whose motion we wish 
to describe, as singularities of the gravitational field. This leads 
to the use of the good <5 functions, which help us to avoid the process 
of renormalization. We shall say that if such point particles are 
assumed, we reason on the level of field singularities. 

To deepen the method of using the <3 function we shall have 
to ascend to the level of continuous distribution. That is, we shall 
assume continuous fluid matter. This problem is of some interest 
in itself. We may, and shall, descend by a limiting process from this 
level to that of field singularities. 

Finally we shall often use the still more general level of an 
arbitrary energy-momentum tensor, especially for discussions in 
which its specific form is of little importance. 

Thus on the level of field singularities, we have: 


_ 8 tcIc vi a r ■ a 

W = =- 2 j J 


N 


A 


A 

u v 


A -1 


ds A As, 


(2.3.1) 


Ai 

■ We recall that E'(s A ) is the world-lino of the A’th singularity; 

A 

m (0) .is its rest mass; because of the use of the good <5 functions, 
ds A is the regularized time-space interval of the A’th particle. 

On the level of continuous distributions we have: 


®~ = [(*+ ■££«"]■ <■ 


fy>\ 
a 


.3.2) 
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Here the density g is a function of the pressure p; u a (x) is the 
Eulerian field of the velocity vector. These quantities satisfy the 
conditions expressing the conservation law of matter: 

g afi u a % p = 1, (Qu a ). a = 0 (2.3.3) 

(if we assume the first of these-equations then the second one can 
be deduced from the conservation law 7'"., == 0 by contracting it 
with %). 

On the level of an arbitrary energy-momentum tensor we have 

- S'" = -8 nT”{A a , g ^). (2.3.4) 


Here JL° stands for the dynamical quantities describing the 


physical field; like the on the first level and g, p, u a on the second. 
We shall assume, as on the first two levels, that 7'°’(JL, g) vanishes 
outside a sufficiently large volume. 

We ask: with what orders should we start the development 
of the 9 r/0 ' tensor ? We shall answer the question first on the level 
of field singularities. There we have (see I, 4) 

”a AAA 
A~1 


a fem, ol dx° 


(2.3.5a) 


We shall assume £* to begin with the order zero: 

(2.3.5b) 

0X2 

A 

The function ft, starts with the order two. This we see from (2.3.5a). 
We also see it from the Newtonian equations of motion. There we 
have (in the case of a two-body problem) 
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Since the lowest order of both sides of this equation must be the 
same, we conclude: 


Therefore we see: 


p — ft+/i+p+... 
2 3 4 


. ST 00 = 9 00 +J 00 +7 00 +..., ) 

2 3 4 


cj-On _ cron 

7 mn = 


7 #n +7 #B +..., 

3 4 

7 mn +.... 

4 


(2.3.7) 


(2.3.8) 


In a similar way we obtain the same result on the level of con¬ 
tinuous distribution, and we assume that on the level of the arbi¬ 
trary energy-momentum tensor we have, exactly as in (2.3.8): 

7°° = 0 (A 2 ), 7°” = 0(A 3 ), T™ = 0(A 4 ). (2.3.9) 

Since 7°^ is a function of A a and g e „, the last equation determines 
the lowest order in X of A a . 


4. THE APPROXIMATION METHOD AND THE COORDINATE SYSTEM 

What is the relation of a coordinate transformation to the 
approximation method? That is: what conditions must the coordinate 
transformation fulfil so as not to upset the properties of the metric 
field, either at infinity, or in the order with which its development 
starts? We write the transformation to the starred coordinate system: 

x* a = x a +a a {o>,X), (2.4.1) 

and we assume that this transformation becomes an identity for 
r -*■ oo, that is 

lima"(a, A) = 0. (2.4.2) 

P-+O0 

Then from the formula: 


9afi — V 


(2.4.3) 
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we deduce, neglecting the products a-ft and a-a: 


h'oo — ^oo + 2a°| 0 , 

ft On = ^0n+® |n M |0> 

ft«in ~ ftmn ® \m ® |n • 


(2.4.4) 


From, the last equation, the answer to our question follows imme¬ 
diately. Since, according to (2.2.7), the ft’s start with ft 00 , fton> Km, 

2 2 1 

we see that, if we wish the ft*’s to start the same way, we must have;. 


a° = a°-|-a 0 + rt # +... 
2 3 4 

a n = 

12 3 


(2.4.5) 


We can easily generalize this result. Let us take for (2.4.1) 

®*o = ao+ao, a,*fc = x k +a k ( 2 .4.6) 

Z+l z 

then the first Ji^ affected by such a transformation are: 


\nn ? n ? ^00 * 

Z Z+l Z+2 

Indeed we have: 

^win = h rnn ® in & \m > 

Z Z Z Z 

^ 0 n = ^ on ® |o > 

Z+l Z+l Z+l Z+l 

fcoo ~ ^ oo “t“ - 

Z+2 Z+2 Z+2 


(2.4.7) 


(2.4.8) 


Let us recall one of the most important theorems of Biemannian 
geometry, which concerns the necessary and sufficient conditions 
for a space to be pseudo-Euclidean. Thus we ask: when can 
a coordinate transformation change the metric form 

g a? dx a dx ? into rj a? dx a dx^i 


(2.4.9) 
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The necessary and sufficient condition for this to happen is 
known to he: 

R^a = 0 (2.4.10) 

that is, the vanishing of the great Kiomannian tensor. The proof 
that this is .a necessary condition is trivial and can be dismissed with 
the re mar k that g a p = T)„p annihilates the great Eiemannian tensor. 
The proof that it is a sufficient condition is much more troublesome 
and many different proofs can be found in books. At least one of 
them is based on the fact that (2.4.10) forms the integrability 
condition for the existence of four transformation functions by which 
is changed into rj^. 

Now, in connection with our approximation procedure we may 
ask a similar question: what are the necessary and sufficient condi¬ 
tions for the existence of a transformation a 0 , a k that turns the ten 
functions Ifl 1 

\nni \n> ^00 (2.4.11) 

Z l+l 1+2 

into zeros? 

Let us introduce the linearized curvature tensor: 

J — i(^«/J|;w+fyt«|a/3 — ^a* w)* (2.4.12) 

Then the necessary and sufficient condition for'a coordinate trans¬ 
formation to annihilate h mn , h 0n , h 00 is: 

l Z+l Z+2 

S man b — 0, 8 mon b ~ 0, $ mo, 0 b — 0. (2.4.13) 

l ' Z+l Z+2 

The proof that (2.4.13) is necessary is trivial. The proof that this 
condition is sufficient is more complicated as it is similar to the more 
general case whore no approximation procedure appears. We shall 
omit it here. 

Having the h#, we can define the h a? through the equation 

00 

g «e = n^+h^, h al> = £h a$ . 


(2.4.14) 
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We can find the successive values of h 00 from those of 7L» (fc < l) 

l k 

since 

a^gpy = %■ 

Thus we find that the h afi start with 


h 00 = -K, h m = h 0m , hT = —hnn, h W = -h 00 , 

2 2 2 2 1 13 3 

h 0<> = (fc 00 ) s —&oo. (2.4.15) 

4 2 4 

The above formulas will be sufficient for all our practical calcu¬ 
lations. Merely for the sake of completeness we quote the general 
formula: 


= 2<-ir w-», 


kp 2 p-1\ 


r, 


(2.4.16) 


where the sum is to be extended over all combinations (of positive 
integers) satisfying Tc 1 -\-lc i -\-...-\-'k v = l. 


6. THE APPROXIMATION METHOD AND THE FIELD EQUATIONS 

We shall now introduce the general scheme of dealing with the 
equations of G. B. T. by the approximation method. It will be 
slightly more convenient — at least in the general scheme — to 
treat as our fundamental quantities, not the h’s, but the y’s. Let 
us therefore start with the definition of the y’s. We write: 

V^g f* - Q ap = ri af + v* (2.5.1) 

which equation gives the definition of the y’s. They will be regarded 
as our principal expressions characterizing the metric field; from 
them all others can be derived: * 


g = \gj = I^T 1 = l^+y^rV (2.5.2) 

from- which follows 


g = l»r+ya 


( 2 . 5 . 3 ) 
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therefore 

V —g = l + £j?a 0 y o|S + non-linear expressions. (2.5.4) 
From the last equation and the definition of V* follows: 

W* = y afi —iy aP y t <,'y ea + non-linear expressions. (2.5.5) 

We see, because of (2.4.15): 

y »o = y oo + y oo +y oo +< _ # ( 2 . 5 . 6 a) 

12 3 

y om = (2.5.6b) 

2 8 

= y™+y™+y™ + .:. (2.6.6c) 

1 2 3 

We shall now consider the fundamental equations of G-. R. T. 

©^ 4 - 871 ^ = 0. (2.6.7) 

We shall start by calculating explicitly the linear part of that 

is the part that contains only the second derivatives of the y’s and 
those only linearly. The Einstein tensor density S'" can be put in 
the form 


W = —L= (or or + ©** (2.5-8) 

2 V—g 

where ©*'" is a function quadratic in the first derivatives. 
Therefore, the linear part of ® >u can be written as 

L(®n = JF"**, (2.5.9) 

where 

K luvP — K Ma ’ rfi = \{n w y**+r\rty n —r( iv y , *—ri al> y ia ). (2.5.10) 

It is interesting to note the symmetry properties of the ex¬ 
pression K^. We see that is skew-symmetric in the indices 

fi, a and skew-symmetric in the indices v, (H. We see that is 
symmetric with respect to the simultaneous replacement of n by v 
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and a by /8- Finally satisfies the relation 

= 0. (2.5.11) 

Thus we see that the expression KT* has all the symmetry prop¬ 
erties of the great Biemannian tensor. 

We write the fundamental equations of G. B. T. in the form: 

+ = 0, (2.5.12) 

where the symbol N[8) means all the non-linear terms in 8. 

Let us now discuss briefly the transformation properties of 
Since a? is a tensor density, we have: 

dx 

0*" r = g a VV*V det ^p 

Into this equation we introduce 



(2.5.13) 

dx 

“a/- 1 

(2.5.14) 

Neglecting the products aa and ay we conclude: 

y* ,4V = y MV +r J ^a\ + v av a\^r l ^a\, 

(2.5.15) 

or explicitly: 

y = y +a | 0 — a |a 

(2.5.16a) 

*0n _ on „0 i n n 

y ■= y (l | 0 

(2.5.16b) 

Al *mn _ ,.mn ~n \ *mn n 8 j_ 

y — y |» —ft i m^r 0 a \sto ft jo- 

(2.5.16c) 

We shall not npset the orders with which we start the 
develop the a’s in the same way as before, that is: 

y’s, if WO 

<i° = ft 0 H” ft 0 +. •.» 

2 3 

(2.5.17a) 

a = a a +.... 

(2.5.17b) 


1 2 



II. S 


THE APPROXIMATION METHOD AND THE FIELD EQUATIONS 


81 


After these introductory remarks -we pass to the chief problem 
of this section "which is to adapt the approximation procedure to 
the field equations of G. E. T. 

We rewrite (2.5.12): 

K aa - ob lab +N(®°°) + 8tc 7 00 = 0, (2.5.18a) 

K oa ' nb [ab +K oa ’ nt {a o +^(©°”) + 8rt'7 an = 0, (2.5.18b) 

+ir M -^ |06 +/r' ia ’”° |a0 +z OTO - m0 100 +f(®™)+ 8>n7 mn 

= 0. (2.5.18c) 

Writing (2.5.10) more explicitly wc have: 

XT** = i ( - fTy"*- 6 mb y ™+ #"*•/* + 6°* y mn ) , (2.5.19a) 

K 0*,r* = |.( _ fTyOb + gObyOn^, (2.5.19b) 

K 0 a,ob = y 00 ) . (2.5.19c)- 

Thus K ,m ’ nb and -SC 00 - 06 start with the order one and K 0a,nb with 
the order two. But, since 7 mn starts with the order four, 7°" with 
the order three, and 7 00 with the order two, if we start this time 
from the bottom w r o have in the lowest order: 


K ,na ’ nb lnf> = 0, K oa - ob [at) = 0, if 00 -’ 


0. (2.5.20) 


The general solution of the first of these equations gives 

■mn _ h m _L h n — A i,m h 8 
7 — (> |«i 0 0 |s> 


(2.5.21) 


whore the b' s are arbitrary fimetions. Comparing this with (2.5.16c), 
we see that y mn can be annihilated by a proper choice of the coor- 

i 

dinate system. We assume this to be done, which leads to y mn = 0 

i 

and if 1 " 0 -* 6 starting at least with the order two. Because of this 
result, wo obtain for the differential equation, (since S’ 00 is of the 
order two): 

Jr 0a '°% = i/V = 0* (2-5-22) 


Motion and relativity 
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Similarly, the third of the equations in (2.5.20) gives 

= *V (2-5.23) 

2 2 

This also shows that y°" can he annihilated by a proper time 
s 

transformation. Thus the result is: merely by restricting our coordi¬ 
nate system we can always choose 

/“* = /<> = /» = o. (2.5.24) 

112 

Let us now go one step further, again starting from the bottom 
of (2.5.18). We find: 

JP" 8 -" 6 ,,* = 0, +8ti7 00 = 0. (2.5.25) 

2 2 2 

Thus now, as before, we conclude: 

y mn = 0. (2.5.26) 

2 

But this is no longer true for y 00 . It must satisfy the equation 

2 

£r°V+8 je7 00 = 0. (2.5.27) 

2 2 

Therefore y M #0. Now about y 4 ". Since the non-linear terms must 

2 3 

be of at least order four, we may omit them. We then have: 

E 0a ’ nb lab +Z 0o ' n0 |0a +8^7°" - 0, (2.5.28) 

3 2 13 

which also gives us a value for y 0n which is different from zero. 

3 

Turning once more to (2.5.18c) we see that the first y mn which must 
be assumed different from zero is y™, since K m0 ‘‘ nb [bo , K m0 - n0 m , N( ® mn ) } 

and 7 rrm start with the order four. Thus in a properly chosen 
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coordinate system we may assume 

y°° = y 0# +/ 0 +y 00 +..., (2.5.29a) 

2 3 4 

y on = /•+?**+. ? (2.5.29b) 

3 4 

y wm = y mn +.... (2.5.29c) 

4 

Because of (2.5.5), this induces in the &’s the development 

h w = h°°+h m + h°°+..., (2.5.30a) 

2 S 4 

h 0n = A 0n +7i 0n +..., (2.5.30b) 

3 4 

h mn = h” m +Ji mn +h mn +... (2.5.30c) 

2 3 4 

where 

A 00 = — A« 0 = iy°°, *”*" = -femn = i^y 00 , = fe on = y°". (2.5.31) 

2 222 2 2333 

Therefore we know how to start our approximation procedure 
The next question is how to prolong it. That is: we assume 


00 00 
y i y i 

2 3 


.,y 00 }y mn ,y mn , 

l-l 4 5 


9 


y ™*- yOn 

l-l 3 


J y J 

4 


y 01 * (2.5.32) 

: 


as kaown. The question is how to find y nm , y 00 , y°". We shall look 

i i i+x 

for the y’s in the order in which they have just been written out. 
Let us write the field equations of G-. B. T. in the 7’th order starting 
from the equation corresponding to (2.5.18c): 


2T’* a,n6 |a & +K m ’ nb [0b +K ma - n ° +Z“°> n# l0 o +N( @ mn ) + 8 *3™ = 0. (2.5.33) 

l I-l 1 I-l 1 1—2 2 i 

The y mn are only present in the first expression. All others are 

known. Indeed JV(©"* n ) can contain the y’s mostly to the order 7—2. 
And now, a few words about ST””. In the case of point particles 

we have according to (2.3.5): 


N 

v-i JL 




A 

p = 




A 

Tt W( 0) dar 
o 2 dsj 


(2.5.34) 
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How, we have: 

/z = ^+^4-^+... (2.3.7) 

2 3 4 

However, if we start the A’s according to (2.5.30), we see that ^ = 0. 

3 

Therefore: in 7 wm only expressions in the y’s to the order l—i can 
z 

.appear. Similarly in 7° n only expressions to the order 1 — 3 and fin- 

i 

ally in 7 00 only those to 1—2 can appear. Furthermore, because 
z 

= 0, the development of 7°^, if we develop in it only the /z’s, is: 

8 

joo = '.7 00 +'7 00 +..., (2.5.35a) 

2 4 

yon, _ cj-ttn -j.... ? (2.5.35b) 

3 5 

7 mn = J mn +J mn +.... (2.5.35c) 

4 6 

We assume that the same is true on the level of an arbitrary 
momentum tensor. That is: we assume that depends on the 
y’s in such a way that only the known y’s appear in it. Besides, 
we assume that generally '.7 00 = 9’ 0n = J mn = 0. Therefore wo can 

3 4 5 

rewrite the differential equation for y mn in the form: 

i 

= Q nn , (2.5.3B) 

• i i 

where ■Q mn is a known function of the known y’s. Because of the 

i 

symmetry properties of K ma - nb , we have 

= (2.0.37) 

Therefore 1 

0"*,. = 0 (2.5.38) 

i 

is a necessary condition for equations (2.5.36) to ho solvable. 

We shall discuss this condition in the next section. Wo shall see 
that it is connected with the equations of motion. However, at the 
moment, we shall simply assume that (2.5.38) is valid. 
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Prom (2.5.21) we know the general solution of the homogeneous 
equation (2.5.36). Denoting by y* mn any particular solution of 
(2.5.36) wo have the general solution: 

. mn *mn , %.m . in aim is /o s oo\ 

y —y +0 |»+» |»»— O 0 (2.5.39) 

Z Z I I z 

We may look at this equation in two different ways: either we 
regard the ft’s as arbitrary functions added to y mmn so as to obtain 
the general solution, or we regard them as a transformation from 
the starred to the un-starred coordinate system. Indeed, a com¬ 
parison of the last equation with (2.5.16) makes this point clear. 

Since y* mn is any particular solution, it is convenient to assume 
that 

/"% = 0. (2.5.40) 

z 

Indeed, if (2.5.40) is not fulfilled we can find a new solution 

y** mn = y* nm + b* m , n + b** m - S™ b** a , (2.5.41) 

z z I I 

so that 

y**™ = 0. (2.5.42) 

z 

To do this, it is sufficient to solve the Poisson equation 

- ~y* m V- (2-5.43) 

z z 

Thus, we can always assume that (2.5.40) is fulfilled. Then the 
differential equation (2.5.36) reduces to a Poisson equation 

= 2Q mn , (2.5.44) 

and we find the general solution by adding the solution of the 
homogeneous equation K ma ' nb lab = 0. 

Having thus disposed of the y™ 1 wo turn to the y #0 ’s. The proper 
equation for them, according to (2.5.18a), is 

= -2V r (© 00 )-8jr 7 00 
i i z z 


(2.5.45) 
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where i2 0# is a function of the already known y’s. According to 

i 

(2.5.19c) we have 

A 00 ' 06 ,^ = i/V-irV (2-5-46) 

X l l 

But the y" 6 are already known from the previous step. We found 

i 

them as the general solution of (2.5.36): 

= y** + &%+& 5 „-<5°X (2.5.47a) 

i i i i i 

y ** {b = 0, y* 06 ,*, == 2i3“ 6 . (2.5.47b) 

l l x 

Introducing this into (2.5.46), we have: 

iy°V = iJV+f 0 - (2-5.48) 

Or, if we call y* 00 the solution of 

i 

/ 00 18a = 2fi»°, (2.5.49) 

we find 

y ° # = y*o° + 6 . (2.6.50) 

i z i 

which, according to (2.6.16a) can again be interpreted as a change 
in a coordinate system. 

Finally we have to consider the equation of the order l+l for 
y°*. This is: 

4+1 K M ’ n \ ab +K 0a ‘ M m +N{®' >n ) + &7i3'' M =0. (2.5.51) 

i+i i+i i+i i+i 

Again, according to (2.5.19b), only JE 0 *** 6 contains the unknown 

Z+l 

quantities y 0 *. On the other hand, according to (2.5.19c), 2r 0a,n | a0 

l+l 

contains the quantities y mn and y 00 just calculated. We denote: 

i i 

Z+l Z+I z+l 


(2.5.52) 
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Then, introducing the values for y mn and y 00 already calculated, 

i x 

that is (2.5.47) and (2.5.50), we have, because of (2.5.19) and 
(2.5.18b): 


£ (y°*+& n |o)i«—£(/“+&“<>) 
2+1 2+1 2+1 2+1 


| an 


= Q^ + iy* 00 
2+1 2 


| 0 »- 

1 


(2.5.53) 


Here the necessary condition for these equations to have 
a solution is the vanishing of the divergence of the right-hand side, 
which contains functions already known. That is: 


i3 0n ln +i2 00 [0 = 0. (2.5.54) 

2+1 2+1 


We shall discuss this condition again later but for the moment 
we shall assume that it is fulfilled. 

Let us now designate by y* 0 * 1 any particular solution of 

2+1 


= £°* + £/ 0 V (2-5.55) 

2+1 2+1 2 + 1.21 

Then 

/* + &% = + (2.5.56) 

2+1 2+1 2+1 2+1 

whore b° is an arbitrary function, is again a general solution of 
2+1 

(2.5.53). To the last equation we may easily add the condition: 

y * 0a =y *°" + /o° 0. (2.5.57) 

2+1 2+1 2 1 

£E this condition is not fulfilled, we can find a different particular 
solution 

y**0» = y*o» + 6 *0 (2.5.58) 

2+1 2+1 2+1 

fulfilling this condition. It is only necessary to choose 5*° so that 
it is a solution of Poisson’s equation: 


J,*0 _ *o« 

0 inn-y |n 

2+1 2+1 


y*°° 

l-hl 


10- 
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For a y* 0a fulfilling (2.5.57) the equation (2.5.56) again becomes 
a pure Poisson equation: 

y* m lBS = 2 O'™. (2.5.59) 

z+i tz+i 

As before, we may interpret (2.5.56) as a coordinate transformation. 

Let us summarize briefly what has been done in this section. 
Firstly, we have shown that by choosing a simple coordinate system 
we can start our development with y 00 , y®*, y mn . Then, having 

2 3 4 

started it, we can prolong it indefinitely, assuming that the inte- 
grability conditions, as expressed by (2.5.38) and (2.5.54), are 
fulfilled. We can find particular solutions step by step by solving 
only Poisson’s equations. Then we can find the general solutions 
by adding, appropriate arbitrary functions. Their addition may be 
interpreted as passing from an especially chosen to an arbitrary 
coordinate system. Of course, we do not know whether such an 
approximation procedure converges; nor do we know whether it 
behaves decently at infinity. 

We have shown that we' can always find a coordinate system 
such that the y’s start according to (2.5.29). But we can show 
more; we can show that by a choice of an appropriate coordinate 
system we have: 

r oo = y oo + r oo _|_ yoo 4 .... f (2.5.60a) 

2 4 5 

y°“ = /*+/* + /”+..., (2.5.60b) 

3 5 6 

y mn = y” m + y mn + y mn +. .. . (2.5.«0c) 

4 6 7 

This can easily be shown on the level of field singularities. 
Indeed because y = 0, the corresponding 7 00 , V7°”, ‘7”*" vanish. 

3 3 4 5 

Moreover N (® 0# ) vanishes, because it is at least of the fourth order. 

3 

2i r (© 0tt ) also vanishes since the only contribution to it would have 

4 

to come from expressions like y 00 | 0 y 00 , n which are of the fifth order. 

3 2 
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H (®" m ) also vanishes since the only contribution to it would have 

5 

to come from expressions of the type y° m | 0 y 00 |n which are of at least 

4 2 

the sixth order. If we wish to obtain the post-Newtonian equations 
of motion, then for practical reasons only the disappearance of y 00 

3 

is of some interest. 


6. ON THE TWO FORMS OF THE EQUATIONS OF MOTION AND THE 
INTEG-RABILITY CONDITIONS 

We shall now return to the general theory. We remember that, 
again on the level of field singularities, the equations of motion 
for the J.’th particle were 

i dx° , 

< 3 -°> 

A 

Q 

A 

where Si was a small neighbourhood surrounding the JL’th singular¬ 
ity. This form is derived immediately from the symbolic equation 

= 0 ( 2 . 6 . 2 ) 

which is a consequence of both the field equations and the Bianchi 
identities. We shall regard the last equation as that of motion on 
the level of continuous distributions and the arbitrary energy- 
momentum tensor. 

We shall now introduce a different form of the equations of 
motion historically prior to this one; in fact it was from this form 
(which we shall introduce now and return to in the last chapter) 
that all the development of the equations of motion in G. B. T. 
originated. 

We write the field equations of G. R. T., according to (2.5.18), 
in the form: 

©"* +8jt9> n = + 

+.K>‘°’ no m -f 8 7i\7“ n = 0. (2.6.3) 
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Nov let ns differentiate the last equation with respect to as". 
Because of the structure of the first two expressions on the right- 
hand side, their derivatives disappear and what remains is: 

+E tia - n \ an0 +X f ‘°- M {M0 +N(® tm ) [n = 0. (2.6.4) 

These are four equations gained from the field equations which we 
assumed to be satisfied; they owe their form to the structure of the 
field equations. 

Also because of the structure of the field equations as expressed 
by Bianchi identities, if the field equations are fulfilled, we obtain 
the equations of motion: 

+$"*<>+ {^} = 0. (2.6.5) 

The two equations are identical in their physical contents. 
Both form conditions that must be satisfied; both are different 
expressions of the equations of motion. 

Let us split (2.6.4), writing out separately the equations for 
ju = m and for /< = (): 

8rc7 mn |n + +E m0 - na lam Voo+A (©”"% = 0 (2.6.6a) 

8rc7'V+^’Vo+^(© W V = 0. (2.6.6b) 

If l had been written under the first equation and l+l unde- 
the second, these would have been.precisely the integrability condir 
tions mentioned in the last section in (2.5.38) and (2.5.54). Thus 
there is a very intimate connection between the equations of motion 
and the integrability condition. But there is one question which 
remains to be answered: how can we cut equations (2.6.6) or equa¬ 
tions (2.6.5) into pieces corresponding to different approximations? 
This question was the essential stumbling block in the work of 
Einstein, Infeld and Hoffmann. It seemed that cutting up the equa¬ 
tions of motion would give an infinite number of contradictory 
equations. To answer the question, or rather to find one of itB 
possible answers, we must deepen our approximation procedure. 

Since we are more familiar with the form (2.6.5) than (2.6.6) 
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it ■will be better to explain this procedure for the equation 7°^$ = 0. 
Also, to make the ideas as simple as possible, we shall remain on 
the level of field singularities. We know, then, that by integrating 
over a small volume surrounding the A’th singularity we can obtain 
equations in the form: 

A 

{i} A.= «. (SAT.) 

A 

(2.6.7b) 

How let us isolate from (2.6.7a) all expressions of the Z’th order. 
d 

These are in —y (ji £ ,„): 

(IX 


A A, 


A A , 


A A , 


A A,. 


A A , 


j“|0 s |0+(« f* |00+(“|0 ? |0+i“ f fc |O0 + »..+ £*„• 

3 1- 3 21—2 S 1-5 4 1—4 l-l 1 


( 2 . 6 . 8 ) 


In the second expression we shall have in the lowest order: 

with /“(»“,1°). 

2 2 2 

A 

Because y D0 depends on £ a , this will also contain contributions up 
2 

to the order l. They will come from the development of £ a inside y 00 . 

2 

A 

Because /xy is of the order 4, we shall have to develop £ a up to the 
2 2 

order l— 4. In other expressions (in 2.6.7a) we shall have to develop 

A 

£“ up to an order lower than l— 4. 

A. 

Therefore, the highest order expression appearing in £ is of 
the order 1—4. Therefore, wo can now regard the equation 


A- 


a 




|0 


= 0 


(2.6.9) 
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as defining £ k . This means: the equations of motion can bo cut 

Z-4 

into different orders, each of which allows us to calculate motion 
with an increasing accuracy. The equations of motion of order Z 
allow us to calculate motion up to the order Z—4. Thus, if we put 
Z = 4, we have the equation for £ s , that is the Newtonian motion. 

O 

Putting Z = 6 we hare the post-Newtonian motion. 

2 

Now, what about the equation of the order Z+l: 


M,o+ f ‘ = Of (2.6.10) 

i+i ’ 

This equation simply gives us fi! If we take the order Z+l =3 

z 

we obtain simply (i = constant and for Z+l = B we obtain the 
2 

first additional, non-Newtonian; expression. 

Thus we see that we may cut the equations of motion and obtain 
the motion with an increasing accuracy. We must assume the same 
situation to be true on the level of the general energy-momentum 
tensor. There is the dynamical characteristic A a which plays the 
role that £“ does on the lower level. It may also be found by the 
approximation procedure with an increasing accuracy. 

All that has been said with respect to the equations of motion, 
applies just as well to our integrability conditions. They determine 
f* and /z. This can be seen from (2.6.6a) if we write Z underneath: 

Z-4 z 

87t7” l ”'| n +£” ia * n ' 0 |Ort0 +i wl0,nO | On0 +jr mo,n0 | n0 o + N (© m,l )/n = 0 (2.6.11) 

l l l l 


In 7 mn the highest order of ! fc is Z—4: 
z 



( 2 . 6 . 12 ) 


Similarly in the y’s appearing in the and ®’s, I” must bo devel- 
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oped up to the order J — 4. Thus we have exactly the same situation 
as before. 

. Summarizing: the integrability conditions are equations of 
motion cut like the field equations. They allow us to find the mo¬ 
tion with an increasing accuracy. 

The procedure outlined here is almost unnecessary if we do 
not go beyond the post-Newtonian motion. In general it is some¬ 
what “messy” since it mixes up the development of the y’s with 
that of the f’s or, in the general case, with the J.’s. 

In the original paper of Einstein, Infeld and Hoffmann a dif¬ 
ferent procedure was used. The authors of this book later perfected 
this procedure and called it the dipole procedure. From the practical 
point of view it is of little use, but the theory built around it is 
more elegant than that first given, because it avoids the development 
of the motion itself. We shall sketch the use of the dipole theory 
only briefly in words rather than in detailed calculations. 

The dipole method was outlined in the last section of the first 
chapter. We recall that the idea of this method consisted in the 
introduction of an artificial vector field D a . It was introduced in 
such a way that the Bianchi identities were identically fulfilled. 
This means that the equations of motion were identically fulfilled 
for arbitrary motion. In this way, by introducing the artificial 
field D a , by abandoning Einstein’s equations, wo tore apart the 
connection between field and motion. We could find the field for 
these non-Einsteinian equations and the motion remained arbitrary. 
Thus nothing stood in the way of cutting the equations and finding 
the field and the D’s step by step. Thus we could find step by step 
the /**, y #0 , y°”, D k , Jl°. Let us assume that we stop here. Then 

Z Z Z-{*1 Z Z+l 

we can change the solution of the non-Einsteinian equations into 
that of the Einstcinian equations by putting 

./>*+ +I> k = 0, (2.0.13a) 

4 5 Z 


/>«+./>« 4 ...+/>» = 0. 


(2.0.13b) 
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This means that only at the end of our approximation procedure 
do we find the equations of motion, which had been arbitrary until 
then. Indeed, from the last equations it follows, according to (1.10.11) 
that: 

• • +9% = 0, (2.6.14a) 

4 5 l 

9^+9 %+... +9% = 0. (2.6.14b) 

3 4 1+1 

Thus by this method we obtain the equations of motion at once 
when our calculations of the field are finished. 

In principle, therefore, two methods of procedure are possible. 
The first consists in developing the masses and the motion in series: 


fi — p-}-..., 

2 4 


I* - !*+£*+.... 

0 2 


(2.6.15) 


By cutting the equations of motion into pieces, we can find con¬ 
secutively . 

ft, H, ft, •••; I*, £*, f*,.... (2.6.16) 

2 4 5 0 2 3 

We shall say that such methods consist of the use of the approxi¬ 
mation procedure of a certain order. 

In the other method, as explained for the dipole procedure, 
we do not split the equations of motion. We regard 

ft, I* (2.6.17) 

as arbitrary, which, at the end of our calculations, we find at once 
up to a certain order, that is, we find: 

iu , f*. (2.6.18) 

—►n+2 —►n 

The arrows denote that the /j’s and f’s are not of a certain, order 
but up to a certain order. 

As far as our calculations of the Newtonian and post-Newtonian 
motion are concerned, there is little difference between these two 
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methods. We shall proceed in the following way: we shall start 
with the [Newtonian equations of motion: 

7^ = 0; 9' OT % = 0. (2.6.19) 

3 ’ 4 ’ 


Then we shall proceed to the post-Newtonian equations of motion: 


= 0 , 

3 5 

9 m, ^+ c r mv . v = o. 

4 * 6 


( 2 . 6 . 20 ) 


However, into the V and 7 we shall introduce the Newtonian 

5 i 

solution already known from the previous approximation steps. 

These remarks finish our general theory of the approximation 
method and we may now pass on to actual calculations. 



CHAPTER III 

THE NEWTONIAN AND POST-NEWTONIAN 
APPROXIMATION 

1. THE NEWTONIAN APPKOXIMATION 

In this chapter we shall find the actual equations of motion, 
using the approximation method already described. We shall start, 
of course, with the Newtonian approximation, and shall find both 
the Newtonian and post-Newtonian approximations on the three 
levels outlined in the previous chapter, starting with the level of the 
general energy-momentum tensor and then descending to the level 
of continuous and finally, discontinuous mass distribution. Thus 
the Newtonian equations of motion, that is those of the lowest 
order, are: 


'.7’ 00 ,o+-7 # “|« = <>, (3.1.1a) 

2 13 

9 m *+V m v >+ih m tf" = 0. (3.1.1b) 

3 14 2 2 

Therefore, all that must be known in order to write out these 

equations of motion is h 00 . Rut, according to (2.5.31) and (2.5.27): 

2 ' 

= -*y #0 , y°V = -i«*y M . (3.1.2) 

2 2 2 2 

If we assume that 7i 00 vanishes at infinity it follows that: 



(3.1.3) 
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From the ' in S"' 00 we understand that the arguments are x k and a? 0 . 
2 

The arbitrary additive harmonic function vanishes and Q denotes 
the entire space. Therefore, the Newtonian equations of motion 
are: 

= 0, (3.1.4a) 

2 3 


7%+7 m V-- Us?' 00 ? 00 

3 i 4 J 2 2 

1 a 



o. 


(3.1.4b) 


Now we descend to the case of a perfect fluid characterized by 
the energy-momentum tensor: 

7- - 4*^ [(*+-$-[*) *•*'—|r «*] (3.1.5) 

and we assume that 


P=P(e) or 6 = SiP) (3.1.6) 

is given. 

Following our general scheme, we are faced with the problem 
of establishing the order of the first expression in the development 
with respect to A. Since u k starts with the order one, u° with the 
order zero, V—g with the order zero, then it follows from the form 
(3.1.5) that q starts with the order two (as p did before) and p with 
the order four. Thus we have: 


S w = <r, S w = <ro n , S™ 1 = ffV m v n -\-nd mn , (3.1.7) 

2 2 3 21 4, 211 4 

where: 



Therefore the Newtonian equations of motion are: 


+ = 0, (3.1.8a) 

2 1 2 1 


(° r ® m )|0+‘ ( <y ® TO ®")|n+ n \m — 
2 11 2 11 4 



1 



0. (3.1.8b) 


Motion and relativity 


7 
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When, multiplied by c 3 the first of these equations becomes precisely 
the continuity equation; if multiplied by c 4 , the second becomes 
precisely the Newtonian equation for a perfect fluid. 

Now let us descend to the case of point particles. We can do 
this in two ways: either substituting directly into the equations 
of motion 

r * = 2^10^10, (3 - 1 - 9) 

or, as a limiting case of (3.1.8), by putting 

a = JT 1 fj, <5 (x — |). (3.1.10) 

A 

We shall choose the second way, though it is slightly more com- 

A 

plicated. Then by integrating equation (3.1.8a) over Q, and because 
the integral of the divergence can be changed into a surface integral 
with integrand vanishing at the surface we have: 


/i, 0 = 0, fj, — constant. 
2 2 


Then for the same reason, the second integral gives: 


A 

a 


-2 


AS d 

w— 


b7a d£ m I—| 


B \ /,f» 


(3.1.11) 


#) 

(3.1.12) 


which are again the Newtonian equations of motion of point parti¬ 
cles under the influence of the Newtonian gravitational field. Here 
the ' on S reminds us that the sum is to be extended over all .B’s 
different from A. This result does not depend on whether or not 
we used our good <S’s, since the contribution from the term B = A 
is always zero. This we see by interchanging x and x' for A = B 
in (3.1.12). 
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2. THE GRAVITATIONAL FIELD FOR THE POST-NEWTONIAN 
EQUATIONS OF MOTION 

To find the equations of motion in the next, or post-Newtonian 
approximation, we shall have to know: 

cjoo, 9-W Cj™ (3.2.1) 

4 5 6 

since, as was assumed in (2.5.35), 

gr 0 o = 9-o» = y™* = o. (3.2.2) 

3 4 5 

As far as the gravitational field is concerned, we shall have to 
develop 

{"}?" aid (“J?" (3.2.3) 

up to the fifth and sixth orders. To do this we must know: 

Am»> ^onj ^00 • (3.2.4) 

2 3 4 

Therefore we shall calculate them in turn. 

We know from (2.5.31) that 

hmn = -£<TV° = <r B fcoo = —2<5” m fdx'j"' 00 |x — x'l -1 . (3.2.5) 

2 2 .2 ■'2 

For A 0 ” wo have, because of 

3 

A 0 * - /” = h w , (3.2.6) 

5 3 3 

and because of (2.5.14) and (2.5.28): 

iW-iW= -s^+^Vo- ( 3 - 2 * 7 ) 

$ 3 3 2 1 

The condition of integrability of this equation is, as can be seen 
from (3.1.1a), 

sr ## 10 +s% = 0, 

2 13 


(3.2.8) 
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which, is the zero equation of the Newtonian equations of motion. 
Assuming now that 

y\ = 0, (3.2.9) 

we have from (3.2.7): 

iy'on m = -8^°“. (3.2.10) 

3 3 

The general solution is, by Chapter II: 

Von = 4 f dx'-^ - —- +a 0 ,„, (3.2.11) 

3 a l x ~ x l » 

where a 0 is an arbitrary function, if necessary fulfilling certain 
conditions of continuity and vanishing at infinity. 

It is much more complicated to find fo 00 . The simplest way is 

4 

not through the y’s, that is not through the Einstein tensor 
but through the Bicci tensor density W. Let us write once more the 
principal equations of G. B. T. 

< X a/3 -ig afi 9l = -8 jrT* (3.2.12) 

Prom these we find 

* = g^afi = 8n7, 7 = g+T*. (3.2.13) 

Therefore, instead of (3.2.12) we may write: 

W* = -8 (3.2.14) 

Let us once more write R m explicitly: 

Por brevity we write: 

K=<P= -2 fdx'T 00 |x-x'|-i. (3.2.16) 

* o 2 
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Therefore up to the order two we have, because of (3.2.5): 

V—S = 1 + H*oo—ftes) = 1 —(3.2.17) 

Since 

9t M = (3.2.18) 

we have, for 9? 00 of the order four: 

W° = -S<pB 00 +R ol) . (3.2.19) 

4 2 4 

Thus the problem of finding W 0 is reduced to that of finding JR 00 

4 

up to the order four. A straightforward calculation gives: 

^00+^00 = -h<P\*>— iftooi*s+^o 4 |o*-|9 ;, |oo + i9’i«7 , i«—(3.2.20) 

2 4 4 3 1 

Because, according to (3.2.11) and (3.1.1a): 

^©9|os 8=8 2^|oo+® o|om (3.2.21) 

we can write 9? 00 in the form: 

4 

9l 00 = —i(^oo — 2a 0 j 0 — + (3.2.22) 

4 4 

Now, on the right-hand side of (3.2.14) we have 

9t 00 = -8jrC7 00 -i^ 00 7) - -4tt (9 00 +7 8$ ). (3.2.23) 

4 4 '—v—' 4 4 

4 

Therefore our differential equation becomes: 

(^oo—2a 0 | 0 — -J-^ 2 ) = ^100+W|aa“t“87r(9' 00 +'7’ S5 ). (3.2.24) 

4 4 4 

Because of the form of 9?, (3.2.16), we obtain the following solution 
of the last equation: 

%oo = -2f<2x>-xT^' M 4^ 

4 J 44^ 22 

+2/dx'/dx"[x, x"]7' 00 9"' 00 +2a 010 , 


(3.2.25) 
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where 

r_ / —//l _ ^ _i_^_I ^ 

L ’ ’ J- |x-x'||x'-x"| |x'-x"|l*"-*l 

This formula reveals the stnieture of h 00 and all the relativistic 

4 

influences upon it. It consists of four expressions. The first one is 
the usual contribution coming from the distribution of matter in 
the non-hTewtonian approximation. The second would also have 
appeared in a linear theory. It results from the fact that we began 
with Laplace instead of d’Alembert equations. The third expression 
is due to the non-linearity of the field equations in G. E. T. The 
meaning of the fourth expression is obvious. 

We should now like to find out how h 00 behaves at (spatial) 

4 

infinity, (if a oto goes to zero there). If we assume that all matter 
is contained in a sphere with a finite radius, then all the expressions, 
with the exception of the second, obviously go to zero. Thus it is ne¬ 
cessary to examine more closely only the expression 

-/dx'|x-x'| 7 '«o |00 . (3.2.26) 

But we remember the integrability condition for y M which we have 

s 

assumed is fulfilled: 

-7 09 10 = 7%. (3.2.27) 

2 13 

Therefore we may write (3.2.26): 

|x| Jdx' 7 % 0 --^ J®' fc 7'°Vx'+ 

+ 2R/ < * xV * a, '* 7 ' < V~- j ^“7 jdx'x' k x' l \T m m +... (3.2.28) 

The first expression is zero, because S’* 08 vanishes beyond a finite 
sphere; the third and the fourth go to zero like l/|x|. The second 
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expression can be written: 

- ^ Jc7 ,w ^ fc) |80 5x'+ (3.2.29) 

Again the first of the expressions written above vanishes and the 
only one that must be discussed is: 

jT ok [0 dx'. (3.2.30) 

Now, besides the validity of the integrability condition (3.2.27) 
we also assume the validity of the Newtonian equations (3.1.4). 
Therefore, again because 9"°^ vanishes beyond a finite region, we 
have: 

jV 0 Vx' = J dx'Jax"T 00 7" 00 ^ =0. (3.2.31) 

We see that h 00 vanishes at infinity but only if the Newtonian 

4 

equations are fulfilled. 


3. THE EQUATIONS OF MOTION IN THE POST-NEWTONIAN 
APPEOXIMATION 

We assume that the Newtonian equations are fulfilled. Indeed 
the zero equation must be fulfilled, because it is a necessary 
integrability condition. Also the other equations must be fulfilled 
if we wish h 09 to vanish at infinity. Therefore the post-Newtonian 

4 

equations of motion may be written: 

9 % = 0, 9”’*. /} = 0. (3.3.1) 

5 6 

In accordance with what was said in the previous chapter these 
equations must be understood as equations of the fifth and sixth 
order for dynamic variables in the post-Newtonian approximation 
of which 9' < * is a function. Besides, .the 9' 0 ' 5 may also depend on the 
metric tensor. 
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We start by •writing out the first of these equations: 

3 2 

5 5 2 3 


= (7«®+^S r o'>) |0 +(v7 08 +^7 #s ) |s +^ 1 * l 7 08 = 0 

4 2 1 5 3 3 

(3.3.2) 

or, since 

<p = -2 J 

(3.3.3) 

we have: 


e 

® 

1 

* 

|x-x'| /|0 \s j |x —x'l ]|« 



[ dx'7' r >«7 0s ( l 1 J =0. 

J \\x-x\l\s 

(3.3.4) 

Now as to the second of equations (3.3.1): 


C7«l0 i C7«M 
J |o~r^ | a 

5 16 

+{»lr + Wr + 2 ft}r + {™}r = 

2 4 3 2 

(3.3.5) 

Here the Ohristoffel symbols are: 


looj 

2 

— &9V> {oof — ^omio + iwimj 

^ 4 ^ ^ 3 1 


2 i®*j 

3 

| = ^08\m h(m\8 

P 3 3 


\nm\ 

2 

5=51 ^naVim <9?ria9V)* 

(3.3.6) 


In place of the ft’s we must introduce their expressions (3.2.11) 
and (3.2.25), then after taking account of the Newtonian equations 
of motion, we obtain, remembering again that (p stands for an 
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integral according to (3.2.16): 

l a \ = ( joo-qtj™-* f dlx / '7 00 .7 /a0 |x~x'r 1 )|o + 

6*5 3 J 2 3 

+ (7 ab -<p7 ab -i f dx'7 b0 : 7 ,o0 |x-x'|- 1 )i6 + 

6 4 * 3 3 

+ kia(-7 0# +'7**)- f dx' ( - ) [7 M (7 , “+7'")-47 /U 7 u ]+ 

4 4 J \ |X — X | /ja 2 4 4 3 3 

+Jdx'J dx"[x,x',x"] la 7»»'7'<"> .7" 00 -£ J dx' |x-x'| |o 7 ,( >0| 0 0^ 0# , (3.3.7) 
where again 


[X, X', X"] = 


|x-x'||x'-x' 


+ 


X— X X”— 2 


+ 


X" —x X—: 


(3.3.8) 

Since nothing was assumed about the special form of 7°^, these 
are the general equations of post-Newtonian motion. We note one 

interesting feature: a 0 does not appear in it at all! This is so 
s 

because it appears in and in h 0Q , but not in the combination 

3 4 

iKoim — ^om|o and not in the combination h onlm —h mn which are 

4 3 3 3 

the ones that appear in our equations. 

Before we descend to the other two levels we would like to 
bring (3.3.7) into a slightly more convenient form. We see that under 
the integral sign we have only the 7’s with the exception of the last 
expression in which 7' 00 | 00 appears. We can, however, change it 
in the following way: 


7 % 0 7 #0 = i7'%7'%-7!%7 M l0 

= -{7 m {a .7'>%-7'\7% 


(3.3.9) 


The s' derivatives can be shifted (with a change in sign) into 
the function under the integration sign. A similar procedure, though 
slightly more complicated, can be carried out with the r derivative. 
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Thus we obtain: 


-i/«!x'|x- X '| |a 7'%o^ 00 = -*(/^x-*V 7,W7 ° # )lo+ 

-i(f dx'\x-x’l la3 ,7' oe 7° r ) lf +$ f dx'\x-x'\ laa , r T oa '7 or . (3.3.10) 


Therefore, finally, we may ■write our general equation (2.3.7) 
in a slightly different form, more convenient for practical use: 

( 7* 0 —4 fdx'\x-x'l-*T oa< 7 00 +2 fdx , lx-x , l-iT 00 7 0a + 

'S*' 32*^ 23 

-i f dx , lx-x , \ lab s ,0t sr <><>) |0 + (9' a6 -4 J dx'lx-x'\- 1 T oa '7 ob + 

3 2 v 6 3 3 

+2 f dx'lx-x'l-^'ooT^-i f dx'\x-x'\ lae S ,0Cl :7 0b ) ll> + 

- f dx' f - 1 — ) [9"°° (7°° H-ST M ) + (9"°° +7' S8 )7 00 - 49" 0 *7 M ] + 

J \ X— X /2 4 4 4 42 23 


4 - fdx' [ dx"[x, x', x"] |a 9' 00 7 , 0 # 7" 00 + 
J J 2 2 2 

+ ij < 2x'|x-x'| |a6 , c 9" o6 9' oc = 0. 


(3.3.11) 


Now we wish to descend from the general level to that of contin¬ 
uous and discontinuous distributions. We start with the continuous. 
According to (3.1.5) our energy-momentum tensor is in this case: 


This can be written 


where 


7 af = croV-(^+y^)7r 


(3.3.13) 
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We develop (3.3.13) into a power series: 


and we obtain: 


a — or+cr-f..., 

2 4 

V = V + V-f-..., 
1 3 

— 7Z-f- Jfj~.. . t 
4 6 


(3.3.15a) 

(3.3.15b) 

(3.3.15c) 


7 00 = o-n, 7°“ = av n + ov n , 

4 44 5 4123 

7 ™ = ov m v n + ov m v n +av m v n +ad mn . 

6 411 231 213 6 


(3.3.16) 


The zeroth equation of motion gives, according to (3.3.4): 

(cr n—{dx'\x —+ lav n -\- av n — fdx' |x—x , [ -1 CT'ofO n ] 4- 

' 4 4 V 2 2/1° \4 1 2 3 J 2 21 /|» 

(3 - 3 - 17 > 

Now we put into (3.3.11) the expressions (3.3.16) and we obtain: 
lov a -\rov a - f-2 fdx'\x — x l \~ 1 oo'(v a —!2v' a ) -f- 

'4123 221 1 

— i f dx' |x — x' I !o6 , go '?/ b ) + (ov a v b + <™v+ ov a v b +• + 

J 22 1 /I® '41 1 23 1 21 3 6 

+ 2 f dx'\x—x'r 1 o'ov b (v a —2v' a ) — i f dx'jx — x'\ lae .o ov' e v b \ + 

J 2211 i J » *i i h 

~ {ax' , -~| [o'ff+ocr' — o'n—<m'+ cr'oiv’tf+v^v'* — 4?>' s ©*)] + 

J \|X“"X|/|a 24 24 24 24 2211 11 11 

+ f dx ' fdx"[x, x', x’'\a<to' O "+i fdx' |x~x'|,^ 0 or' cro' b v c = 0. (3.3.18) 
J J 2 2 2 J 2211 

These equations determine a and v* if ox is given as a function 

4 3 

of a and v. They form a transition to the more important problem 
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of discontinuous particles. Tet the problem of continuous distribu¬ 
tion may be of some interest in itself; perhaps in the theory of great 
nebulas where the internal velocities cannot be assumed to be small 
compared with the velocity of light; or perhaps in other cosmological 
applications. 

We shall now pass to the more important problem, to the 
application of the general equations (3.3.4) and (3.3.11) to the case 
of point particles. Generally we have there: 

'.7 00 = £ (3.3.19) 

If we develop the p into a power series and the £’s also, wc have, 
since 

7 00 = - yP$J* +yU, d = d(x-i). (3.3.20) 

4 fr 2 2 a * 

Both here and later we shall leave out the indices written under 
each quantity if they are the lowest indices. Thus, for instance, p 
stands for p or, permitting an error of a higher order, for p+p. 
We have: 8 24 


'J tn 

s 


„„ OO OO 0 

7 w =y(p6-pd,n, 

i 7T 4 ! 




000 00 0 c coo 

4 2 3 


(3.3.21a) 

(3.3.21b) 


'r n *=Y (^! m 10 r,.+^rJ n io+ 
« <sr * » 


aoo a 


ooo m c oo c c cr 
+ p6£ % £\-pd\, f 8 f To f,o) - (3.3.21c) 

3 2 


Introducing this into (3.3.2) and (3.3.7) we wish to find 
f 7 0a . a dx and (7 m, . a dx 

A 5 A 6 

Q Q 


(3.3.22) 
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where Q is a small region surrounding the A r th particle. Therefore 
we are faced with integrals of the type 

f dxfdx'f(x, x’) d(x-i)d(x' — i). (3.3.23) 

A a 
a 

Since Q is the entire region, the result is: 

/(#, #)<W (3.3.24) 

Thus this is the simple general rule of calculating the integrals 
(3.2.23). But special attention must be paid to' the case B = <7. 
Namely if / is a singular function, e. g. 

/ = |x-xT\ (3.3.25) 

then we have an expression 


it 


I-II" 1 =t 2 b = oo for A = B. 


(3.3.26) 


But it is here that our argument about the choice of good <3’s 
intervenes. They annihilate this type of contribution. However, 

A B 

if /(|, |) is finite for A = B, then, of course, it must be taken into 
account. This is especially true in expressions of the type 



B,C 


ABO 

fABTAa 


(3.3.27) 


Here we see, for example, that the expressions B = C must not 
be ignored, and we obtain: 



ABO 

fABfAO 



ABO 

w 

fABTAC 



w* 

r S AB 


(3.3.28) 


where S' means that the sum is not to be taken for B = C or A — B 
or A = 0. 
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After these remarks the calculation' is fairly straightforward. 

A A 

We start by calculating from (3.3.2). Integrated over £, because 

4 

of the Newtonian equations of motion, and because the divergence 
expression does not give any contribution, this expression gives: 


AB 

A Vl' UU AA A 

M=\ (3.3.29) 

4 ^ r AB 

We must put this expression into (3.3.21). (The arbitrary additive 

constant can be absorbed by /i.) The rest is a straightforward cal- 

2 

culation which finally gives us: 

[ AA / AA" A . vY , \ A 

M$ p+(W ,°f,o+3 2_i MM r ASJ f°|o+ 

VI' AB B VI' AB B t 

A 4tW T AB£ |0~i 2^J W r AB\i a f 8^ l°]|0 + 

B B 

VI' AB A B 

B 


YY-48.-4. -4 , J? JR .4 B 

/, w(jf 10 f |0 + 2^ |of |«~4{ |0 f®o)(^Js) 1 fa + 

B 

Yl' ABA B VI' AB A B 

+ i 2j l° f W r AB\faf4 r+ %2j MMiM+M)( r AB) l fa + 



1 

^BO 


1 1 , 1 1 \ 

-1-.4 

r BO r CA fa A r AB I \t a 


= 0. 


(3.3.30) 


These are the final equations of motion for N particles in the post- 
Newtonian approximation. 

Obviously in the expressions of the sixth order which do not 

AAA A 

contain the f’s it is immaterial whether we take f® or £®-f*1®. The £® 
2 2 2 
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appear only in the expressions 

/“f ioo- £ MMif(r£) l f g f a + | 8 (r^),f 0 f s ]. (3.3.31) 

All other expressions are already known from the Newtonian 

approximation. Therefore we have 3JN r differential equations to 

determine the 3 N unknown f’s. 

2 

We mentioned before the difference between the equations of 
motion of a certain approximation and those up to a certain approx¬ 
imation. Here we followed the theory and practice of equations of 
motion of a certain approximation. The theory up to a certain approx¬ 
imation requires the use of dipole fields. However, as we also 
mentioned before, the story is essentially simpler if we do not go 
further than the post-Newtonian stage. There the only integrability 
condition that we met was 


t 7 < % = 0 (3.3.32) 

5 

which defines /x. Thus with the proper choice of /x this integrability 

4 ,4 

condition is fulfilled. Therefore, ignoring the Newtonian equations 
of motion, we can write the equations of motion up to the post- 
Newtonian order. They are: 

7 m %-K7 m % = 0. (3.3.33) 

4 6 

Here we have the unknown functions: 

t* = P+P (3.3.34) 

0 2 

and by integrating (3.3.33), we obtain at once the motion up to the 
post-Newtonian order. Since in 17*^ it is immaterial whether we 

A A 6 

put f* or f *, we can regard this expression as known if the solution 
« 

of the Newtonian equations is known to us. Thus we can write down 
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the equations of motion up to the sixth order, introduc ing 

v 00 +r°=Y$+Y$, s = d(x-§). 

2 4 ^ 2 " 4 

Then we have: 

K A , AA A AB \ A jIJB J? 

/“+£/*£ |o£ |0+3 2j PPjb)? |0-i 2j W r AB\t*f*t a \*+- 
B . B 

_ B -I VT- 45 A ~B B A 3 

-42, ~ 2 j wa +1^10^10 +5^10f*,o-4^,0r,o) 

2 ? 2 ? 

1 w AJBA s K~V AB A B 

(rAB)lta + i2 j W* \o£*\or ABl tafsfr + i2j + P)(Xab) ,f« + 

B JS 

v 1 ' ABO l 1 1 , 1 1,1 1 \ 

+i > AW* — — + — -h -= 0. (3.3.35) 

Sfi ' Tab Tbo Tbc r ° A r OA ‘Tab! fa 

These equations can be derived from a Lagrangian. Indeed, 
these expressions have the form J 


where L — L+L and 

4 6 


( £ ,f>>.o L f a ~ 0 


(3.3.36) 


L=\Zihh+\2kiw+ 

A A 

+1Z K)-* Z Or- i *>Pnh+ 

A * B A,B 

^ABiicfb^w^o+i M^ab— ; w{(*+/*)?!»+ 

A * a# Xfi 

l ri 7 AJBO 

~* a£jo ABTBcJrrB ^ r ^ A + r OA r AB) • (3.3.37) 

Obviously this is not a satisfactory way of obtaining a Lagrangian 
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since its existence appears accidental and is not connected with 
the theory. We shall come back to this problem in the next chapter. 

There is one more question which we should like to discuss briefly 
and generally without going into particular calculations. While 
finding the equations of motion we have used our good d’s and 
therefore we did not have any contributions from the self-terms. 
This was because, according to our rules, we put 

(»aa = / dx J<Zx'|x — x'| -1 (5(x — f)<5(x' — £) = 0. (3.3.38) 

A 

a 

The question arises: would the result have been changed if co^ 
were not set equal to zero 1 To find a rigorous answer to this question 
we should have to investigate in detail the transition from a con¬ 
tinuous to a discontinuous distribution. This is a rather troublesome 
and lengthy procedure which we shall therefore omit here. Its 
result is almost obvious. Let us assume that we start from drops 
of dimension Z. Then the post-Newtonian equations of motion of 
the center of mass of each of these drops are, with the accuracy 
0(1), identical with those derived here in (3.3.35) or in (3.3.33). 
For Z 0, there are some additive singular expressions too, of the 
order 0(l~ l ). However, they give a contribution only to the mass 
of order /z 2 /Z. Since is the gravitational radius of this drop, then 
we obtain a correction to the mass of the order 




gravitational radius’ 
dimensions 


(3.3.39) 


This quantity is not observable, because it appears only as an addi¬ 
tional constant. 

Coming back to the question of what would have happened 
if °>aa were not equal to zero, or, in other words, if we had used 
different <?’s and not the good ones, then the answer is: we would 
have to use the renormalization procedure in a way that is similar 
to its use in quantum electrodynamics. Since it can be shown that 
the theory can be renormalized, we see that the use of our good d’s 
avoids the process of renormalization. 


Motion and relativity 


8 
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4. THE CONSERVATION LAWS POE A SYSTEM OP PARTICLES 

From tli© ©xistenc© of a Lagrangian 

L=L+L (3.4.1) 

4 6 

n,Twi from its invariance with respect to certain transformations, 
w© can deduce certain conservation laws. And so: 

The conservation of gravitational energy E can he deduced from 
the invariance of the Lagrangian L with respect to the transformation 

a;' 0 = x 0 + a 0 , a 0 = const. (3.4.2) 

As a definition of the constant E, we obtain: 

2i’ v L A -L = B. 

A l f |0 

Introducing (3.3.37) here in place of L, we obtain: 

|0^ |0 —4 £ |0? |o) + 

A AJB A 

+* y / w^„?,«r2i+ 

A,B £5 

r AB j j} A B 

-1 2j W r ABtf>$ 1 * 10 ^ 10+1 2j «*(/*+/ 0 » Sb + 


A^B 


A# 


+ T £ Wt*{ r AJirBC+ r BO r OA + r OA r lB) = E . (3.4.3) 

AJifl —* 

Up to the fourth order we find the Newtonian conservation law: 

= e. 

A A £b 

A 

if we add JjP = -® as the constant of integration, 

A 2 


(3.4.4) 
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The conservation of gravitational linear momentum can be 
deduced from the invariance with respect to the transformation 


f°+a a , a a = const. 
As a definition of the constant P fc we obtain: 

S L «* -**• 

A I® 

Introducing the proper Lagrangian, we obtain: 

(1+if*|of*io—4 mvab )io + 

A B 


(3.4.5) 


- i h (P-h ?,o rji = P a - (3.4.6) 

'ts — 

The conservation of gravitational rotational momentum J 1 ** 1 
can be deduced from the invariance of L with respect to the trans¬ 
formation 

where M a b is a constant orthogonal transformation matrix. We 
obtain as a definition of the constant J lab] 

V (LA P-Lfa P) - J™. (3.4.7) 

A l» I® -*» 

Introducing the proper Lagrangian, we obtain: 

(1+3 2V^n+KV !*|o) ( P\o P - £|of*)+ 

A B 

-i Yw(P,<f‘-hh'ii+ 

AJ} -* 5 
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We obtain the conservation of the gravitational mass center 
in a slightly different way, since L is not invariant with respect 
to a Galileo transformation. Beoanse of Newton’s equations of 
motion, we have from (3.4.6): 






= P“ 

—>5 


from which follows: 


(3.4.9) 


= PV+<? a (3.4.10) 

A B -+5 

whore Q a is constant. 

There are some interesting physical conclusions that can be 
deduced from these equations. We see that the gravitational total 
mass as defined by the conservation laws is different from the total 
inertial mass, as previously defined. Indeed, let us call 

_ A A 

2 (/*+/*) — the total inertial mass up to the fourth order, 

A 2 4 

E — the total gravitational mass up to the fourth order. 

4 

Then we have, as we see from (3.3.29) and (3.4.3): 


AB 

2 s -*2 2 

A * A A# AB 


4 


\— = /«(<?) = B. (3.4.11) 

A AJ3 TaB \ 4 4 


It is the total gravitational, and not the total inertial mass, that 
is conserved. It is the center of gravity of this gravitational, and not 
inertial, mass which moves uniformly. 

We shall return to this problem in the last chapter. Then we shall 
also discuss an additional, more general problem: can we, independ- 
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entity of the Lagrangian and of the approximation method, formulate 
the conservation laws in a field language, for an arbitrary energy- 
momentum tensor. What this reasoning has shown us is the intimate 
connection between these conservation laws and the equations of 
motion. But the equations of motion can be expressed in a field 
language. The same must be possible for the conservation laws. 
It is this problem to which we shall return in the last chapter. 



CHAPTER IT 


THE VARIATIONAL PRINCIPLE AND THE EQUATIONS 
OP MOTION OP THE THIRD KIND 

1. FORMULATION OF THE PROBLEM 

la Section 3 of Chapter I, we discussed the equations of motion 
of the first and second kinds starting from a variational principle. 
Let us recall briefly how we formulated the equations of motion 
of the second kind. 

We introduced 


ds A - (tpafdh 11 * 

where 

(4.1.1) 

i^ = fdxd(x-t)g aP 

(4.1.2) 

and then we varied the integral 


^ A. ^ A A A 

W = - )° / Cg^di^ 

A—*l 

(4.1.3) 

with respect to f°, where <51° vanished at the ends of the interval. 
The result was 

» ffj , , A . A A. 

v A c\ia d * & 

(4.1.4) 


from which the equations of motion could easily be deduced. 
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All this has already been done. Here -we should like to draw 
attention to one essential assumption which accompanied this 
reasoning. It was: the have to be treated as given functions, 
not to be varied with respect to the ff’s. This means: using the 
definition (4.1.2) and varying with respect to the |’s, we only 
took into account its dependence on the £’s through the d’s. But 
in fact we know that the g'a depend on the f’s and their time 
derivatives: 

Safi = 9^(AA h (4-1-5) 

if we assume, as we do here, that only the first, time derivatives 
of i k enter into the g'a. 

We remember also that the tweedling process consisted in 
replacing the ®’s by the £’s and getting rid of the singularities. 
a 

Thus g,# meant: 

9aP — J7a/j(*°> £*» £ k \01 •••> £*> f |o)- (4.1.6) 

Als o we recall from (0.18), or from what is said in. Appendix 2: 

A A ± ■ A A ■ A 

ffafilO = ffafilO + ffafilS f* 0 = 9afi\, f \|0 > (4-1-7) 

A ± 

9jfa = 9jfr + 9a?[B Sab- ( 4 . 1 . 8 ) 

We now formulate our problem. We are looking for a La- 
grangian 

L = L (®°, l\ h h) 

such that a variation of the action with respect to the f’s, van i s h i ng 
at the end of the interval of integration, would give the proper 
equations of motion of the third kind, that is 

x"0 

d f Ldx° = 0 (4.1.10) 

x'o 

would give the equations of motion of all the particles. 
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Such, a variational principle, if it exists, we shall call Fokker’s 
variational principle since he was the first to introduce it into 
electrodynamics. 

It would seem that a good way to start would he to imitate 
the way we obtained the right equations of motion of the second kind. 
Let us therefore, as the first attempt, start with a La gr a.p g ig.ri L': 


L' = L' (x°,£ k , £ fc , 0 ) = 4o) c . (4.1.11) 

A 

Here, the gj s have to be treated as functions of the f’s, the £, 0 ’s 
and ®°. 

We form the equations of motion from such a Lagrangian: 



If, as before, we introduce: 


■a Jcm (0 . duo 0 



(4.1.12) 


(4.1.13) 


this gives us the following equation, because of (4.1.8): 
(PSaft £ fi \o)\0 — kf*9aP\8 I a |of^|0 + 


hK+i&Ci Khl =o. (4.1.U) 

^ JL I® 


We must distinguish here between the expressions in the first and 
in the second row. If there were only expressions of the first row, 
everything would be splendid. We would then have obtained the 
three equations of motion written under (1.7.8), the fourth being 

the definition of as expressed in (4.1.13). But the disturbing 
factor is the two expressions below. Therefore £' is not the proper 
Lagrangian of the Fokker type. Of course we do not know before- 
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hand even whether such a Lagrangian exists. For example, in the 
case of electrodynamics we know that such a Lagrangian does not 
exist for advanced or retarded action. Thus our aim is to find out 
whether and under what conditions we can get rid of the expressions 
written out in the second lino of (4.1.14). 

If we wish the Fokker Lagrangian to he invariant, or rather 
d(Ld®°) to he an invariant, and also to depend only on s»°, the I’s 
and the f| # ’s, then there is little choice in this matter. Let us assume, 
as we did before, that 

»)• (<LU5) 

Since we wish to have only second derivatives of f*(a>°) in the 
equations of motion, we should avoid higher derivatives of g aP as 
much as possible. Thus our choice is essentially restricted to 


L” - f V-gGdx, 

0(3) 


(4.1.16) 

(4.1.17) 


where the integral is to he taken over the entire three-dimensional 
space. If we write down the variational principle 

s"o 

8 j dca 0 J V —gGdx = 0 (4.1.18) 

X'O 0 ( 3 ) 

we see that the left-hand side is an invariant, since V^gG differs 
from V—gR only by expressions that can he put into the form of 
a (four-dimensional) divergence. 

First we vary V~~gG with respect to g^ • This gives. 




+ 



dg ap \ a 


n s 8 g a + 


rdV^gG 
>1 a ^i° 


«"o 
*'<> * 


(4.1.19) 
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Here 2?^ means a two-dimensional surface at infinity. Now the gr’s 
vary because we vary the f’s and the || 0 ’s in them. Therefore 

1 

B B 

Since we assume that the <3|’s and their time derivatives vanish 
at the end of the time interval, we may set the last integral in 
(4.1.19) equal to zero. But it is different with the surface integral. 
Let ns simply assume for the moment that the surface integral 
vanishes, that is we assume: 


/ 


dV-gG 

&9ap\S 


n s dg a? d^ = 0 


(4.1.21) 


where dg^ is (4.1.20). From the contents of Section 4, Chapter X 
follows the well-known formula: 


V9aP \ Og a fi |y /,y \ 

= =8 n2 it a \Qp\* t. (4.1.22) 

• - A 

Therefore, taking into account our expression (4.1.20) for <5 g^ we 
obtain: 



A-^A A 






(t*9a0 |fj 0 f ' “|0 f *o)|o}d£ # + 6' °|0 f ^0 &£* 

A# 


*"0 
*'o * 


(4.1.23) 


Again the last expression vanishes because of the arbitrary if. 
We have thus again obtained the wrong equations of motion, namely 
a a ^ A 

2 rtf'*- 2 (A* Eif%* *of Vio “ (4.1.24) 

A A 
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Our long argument is almost complete. We see that, if multiplied 
by —fa these are exactly the expressions that caused difficulty in 
(4.1.14). Thus, by a linear combination of these two Lagrangians 
we can obtain the right equations of motion provided that (4.1.21) 
is satisfied. 

Therefore: we can obtain the equations of motion of the third 
kind from a Lagrangian: 


£ = - m (0 )<KMVV /2 + f 

A ' " 0 ( 3 , 

where 

In order that 



(4.1.25) 

(4.1.26) 

(4.1.27) 


may give us the right equations of motion, we assume that the 
variation of and vanishes at the end of the time interval and 
also that the equation 


*2 = 0 (4-1-28) 

is satisfied. We obtain the differential equations of the second order 

A. 

for the 3 N quantities f*(a>°). 

If we introduce the /i’s we do it through 


I 


A 


n = 



(4.1.29) 


which we regard as the definition of the inertial mass appearing 
in the definition of 




A 

d. 


(4.1.30) 
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Thus, if only (4.1.28) is satisfied we can find the Lagrangian that 
we guessed in the previous chapter, using the definition (4.1.25) 
and checking, of course, whether (4.1.28) is satisfied. To this task 
we now turn. 


2. THE LAGRANGIAN UP TO THE SIXTH ORDER 


The calculation of the Lagrangian up to the sixth order is a very 
straightforward affair. But one amazing fact arises in it: we do not 
need to know A 0 « in order to calculate such a Lagrangian. It is 

4 

sufficient to know only h 00 , which, as we know, are very simple. 

2 3 

Indeed, the only complication in finding the equations of motion 
of or up to the sixth order from the vanishing of consisted in 
finding h 00 . This we do not need to do any longer. The only thing 

4 

we must know about fc 00 is that it must vanish at infinity like 1 [r. 

4 

But (though we know it from Chapter III Section 2) we can find it 
without any explicit calculations of h 00 . Indeed h msa can depend 
only on 44 

^00100 = 9=1005 9=ia9=;s 5 (4.2.1) 

2 XX 


But <p\ a <p\t, <P\at<P for r -> 00 are of the order at most —4 in r. There¬ 
fore their contribution to is at most of order —2 in r. Now, about 

4 

the contributions from <p ]W . Bor r 00 we have for <p: 


A A 


(7), 


+... 


(4.2.2) 


Because of the Newtonian equations of motion <p !00 is at most of 
order —3, therefore the contributions coming from this expression 
are at most of order —1 in r. We shall very soon see that this faot 
insures the vanishing of the surface integral (4.1.28) which was 
a condition for the existence of such a Lagrangian. We may remark 



IV. 2 


THE LAGRANGIAN UP TO THE SIXTH ORDER 


125 


in passing that and h m are also of order —1 and all their 

2 3 

derivatives are, because of the Newtonian equations of motion, at 
most of order —2 in v. 

We shall start by calculating V —gG up to the sixth order. A cal¬ 
culation which is straightforward, though slightly troublesome gives: 

V-gG+V—gG = —&>, m 9V + WV9V»"~ 

—- —r- a i 

^Ol\in^om\l ”f~ <P\m ^00| m • (4.2.3) 

3 3 3 3 4 

We see that this expression is invariant with respect to the change 
7*0 m “*■ hQm~^~ftQ\m7 h 0Q -> /i 00 + 2a 0 | 0 (4.2.4) 

3 3 4 4 

which is not astonishing in view of the fact that a 0 , as we know, 
does not appear in the equations of motion. Therefore we may 
assume that satisfies the coordinate condition 

3 

y° 8 \8 ^ 7*os|a ^ y°°|o = 2A 00 j 0 . (4.2.5) 

3 3 2 2 

We can rewrite (4.2.3) in the following form: 

V—gG+V—gG = | W|*> — |9 ,S ?V>+ ih«i\ihtn\m J r 

J l 3 3 

-JW'oa + ?’| S ^'o«+W° |0 + WV, (4-2.6) 

3 3 4 

where W° and W” are: 

3 4 3 3 3 3 

(4.2.7) 

From (4.2.3) we see that the function under the integral (4.1.28), 
the vanishing of which was a condition for tho existence of our 
Lagrangian, is of order at most —3 in r. This is so since dg# is at 
most of order —1; therefore the integral vanishes; therefore the 
Lagrangian np to the sixth order exists. 
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Furthermore we see from the last two equations that forming 
the integral of V—gQ, we can omit the last two expressions in 
(4.2.6). This is so for W° because the variation of and vanishes 
at the end of the time interval. This is so for W”* because W™ is at 
most of order —3 in r and therefore its surface integral vanishes. 
Thus, putting: 

V tt, A 0o|s » = -16 n V i (4-2.8) 

A 3 A 

we can write: 

f+f" = i hhI i '—t + 'E£ )dI ~ t(2 

4 6 A 


~ A - 1 r 1 

+ i^oaf°|o + Poo] + g7 J Ki\ih m \mdx\- (4.2.9) 

* ’ * °(3) 

The last expression in the above equation requires special attention. 
We have: 


r 7 ,# ° 




(4.2.10) 


Therefore we can write: 


f ^ fl h mtm dx = 16 f dx'dx"dxr , *7” m (----) (4.2.11) 

j J \|x— x'||x— x"|/|JW< 


Let us first perform the integration 


C= f 

J |x-x'||x-x"| * 


(4.2.12) 


From this follows: 


C„, - -4*J = - -j^Z^TT. 


(4.2.13) 



XV. 2 


THE LAGRAHGIAN UP TO THE SIXTH ORDER 


127 


therefore 


G = — 2 ji|x'— x"|. 


(4.2.14) 


If we put this into (4.2.11) we can perform the other integrations 
simply as described in the previous chapter (Section 3) and we 
obtain: 


J = — 32 tz l |of TO |o- (4.2.15) 

AJB 


Therefore: 


L"+L" = y t id-&+froat a » + lko'} + 

4 6 4 

-l y ^r AB ^ m pJ m ± (4.2.16) 
Having finished with this part of the calculation, we now find 


l'+l' - - X ™ {0 )C(gJ a J%) va 

4 6 A 

where m ((l) is connected with p, according to (4.1.13) 


A liC® -A -4. 

w (o) — ~7“ (&tf£ a |o£%) 


(4.2.17) 


(4.2.18) 


Therefore 

L'+L’ 

4 6 


A 0* -4 

m (0) = — jU = const. 

aJ 2 


(4.2.19) 


(4 -4 w <4 -4 

yjlXWf’V , i.+fr-*5 s -KfW+ 


+Jflf fl ,« f a |o+^»f°|o + Poo)}- ( 4 . 2 . 20 ) 
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We see from (4.2.9) that the h 00 from L' and L" cancel each other. 

Therefore -we do not need to calculate h 00 explicitly. We also see 
from (4.2.20) that replacing * 

Jwby ^o«+<*oia and Ka by lioo+2« 0 , 0 (4.2.21) 

3 3 4 4 

changes the Lagrangian by 

c 3 ri A ^ c 3 ( ;£\ 

AL=== Sj ( a o\8^ 8 \o + ^o\o) = ^)l° C 4 - 2 - 32 ) 

vrhich does not play any role in the variation of L and therefore 
in the equations of motion. 

Omitting now the constant factor c®/&, we obtain from (4.2.20) 
and (4.2.16): 

A A 


M<P 2 f*ha£ a \o-$£ / J ’/ J ' r AB^fi>$ a l oi b lo . 

(4.2.23) 


A,B 


Here we introduce in place of <p and h^: 

3 


A „ B "■ B B 

V = “ 2 Z = 4j?pr2h $ a |o (4.2.24) 


and we finally obtain 


i - -SUiZthh+i z&a+i £ A?vV+ 


A 

-t 2 j /“l“(A‘ + /“)»*i'i-| 2 j ^ r scf 1 + r BcPCA + 

A P a tT« 


A,B 

/ ABO 


A > B A % B,Q 

if A B B A 


+ r ai y Jp)—2 J}' A</* y iUV“,o+i 

^ its 


+ f |o£ |o) r AB ~i i“A ir AB|f a f 6 # < ’|O^ 6 | 0 - 
4.B 


(4.2.25) 
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This Lagrangian coincides (up to a trivial constant) with the one 
guessed in the previous chapter. 

Thus we have not only found the Lagrangian of the Fokker 
type leading to the proper post-Newtonian equations of motion, 
but we have also formulated the general conditions under which 
such a Lagrangian exists. We saw that to calculate it we do not 
need to know h 00 explicitly. It' was sufficient to know that h 00 van- 

* 4 

ishes like Ijr as r -*■ o o. 

Perhaps a historical remark would not be out of place. When 
the equations of motion were calculated with the help of certain 
surface integrals for the first time in 1938, we needed the knowr 
ledge of: 

^OO > ^mn > ^on, ^00; ^mn j ^On • (4.2.26) 

2 2 3 4 4 5 * 

The first three were for the Newtonian approximation, the much 
more complicated last three for the post-Newtonian approximation. 
Later, we found that using the principle 7 = 0 we needed only. 

^00 1 hon } ^00 • (4.2.27) 

2 2 3 4 

Now we find that to formulate the post-Newtonian equations of 
motion, we need as much knowledge of the metric field as we 
needed in 1938 for the Newtonian equations of motion. Or, today 
-for the post-post-Newtonian equations of motion we need as much 
knowledge of the metric field as we needed before for post-Newtonian 
motion. Today we could calculate the post-post-Newtonian equa¬ 
tions of motion with not too much difficulty if there were any phys¬ 
ical sense in doing so. 


3. GENERALIZATION 

So far we havo considered the motion of spherically symmetric 
non-rotating bodies represented by singularities of the single pole 
type. In this section, we shall show how to generalize the action 


Motion and relativity 
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principle to include the case of rotating bodies which will be 
represented by the pole-dipole singularities of the gravitational 
field. 

We saw from (4.1.22) that the variation with respect to g^, 
if certain boundary conditions are fulfilled, gives the equation 

8 j V^g Gdx = - / dg^dx. (4.3.1) 

Therefore the variational equation 

8 f (i V~^g G-^ng a p-7 afi )dx = 0 (4.3.2) 

(where c = 1 and the gravitational constant Tc = 1), when 7 af does 
not depend on 'g^, gives the equation tor the gravitational field. 

Let us now use (4.3.2) as the definition of our Lagrange equation 
and vary its Lagrangian with respect to introducing as an example 
the case of single singularities — that is 

N A A A 

‘-7 a ' = 2vi a s'8, r = r,o ( 4 . 3 . 3 ) 

Ami 

then we do obtain the right equations of motion, identical with 
(1.7.8) if we assume p as a known function of time. Therefore, in 
general we may regard the result of the variation of (4.3.2) with 
respect to f* as giving us the right equations of motion, assuming 
of course that 7°^ does not depend on and that p is a function 
of time given by the added equation 7 0fi ;? = 0. We may thus also 
expect to obtain correct equations of motion in the case of pole- 
dipole singularities representing rotating bodies and described by 

N A A A A 

— (4.3.4) 

A -1 

4. A FOKKEB-TYPE LAG-BANGIAN FOB BOTATING BODIES 

Our aim is to obtain the Fokker-type Lagrangian leading to 
post-Newtonian equations of motion of rotating bodies. The bodies 
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are now characterized by a set of parameters t and t Tafi . According 
to the general considerations of Chapter m, the developments of 
these parameters for (a, /?) = (0, 0), (a, /}) = (0, r) and (a, j9) = (r, s) 

A 

should start with the orders 2, 3 and 4 respectively. But t r0 ° would 

2 

influence the Newtonian equations of motion through h 00 . To avoid 

2 

A 

this complication, we shall start t m with the third order. 

The energy-momentum density (4.3.4) does not lead immediately 

A A 

to the correct Lagrangian. The coefficients t ap and t ral> must first 
be expressed in terms of the f° and their time derivatives. This 
can be done by means of the equations 

s-*, = s- |0 +a«v+s” {»} +25:”-{“}+s:"{ r “} 

= i[(t- l .+f“{4+2>(«} >{;}) 4 

[A n A A A A 

-[t ra \ 0 +t a0 r^-t aT + 

+< M0 W+^{; o }+^{“})t + 

+ (t ,a0 1 S [0 -t ra8 )ts], ( 4 . 4 . 1 ) 

which by integration over the three-dimensional neighbourhoods 
of the singularities lead to 

A A 

s sv* - 4+ 2 * , '{”}+ 

a 


A A A A 



(4.4.2) 
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A 

JV- f r )Z e = r“° 10 +f o0 | , 1 o-1 ar + r 00 { o a o } + 

o 

A A 

+ 2 * U {«} + ‘‘"‘{j}-°’ < 4 - 4 - 3 > 

J = r a0 l , |0 +‘t s “°l r ,0—f ras —= o. (4.4.4) 


A 

a 


Equations (4.4.2) for a — 0 give in the third order: 


t* |0 = ^ 0 = 0. (4.4.S) 

3 3 

Similarly (4.4.3) give 

t 0r =it r i 0, (4-4.6) 

3 2 

and (4.4.4) give 

t rw +f aor = 0 . (4.4.7) 

3 3 

Introducing the notation: 

fros_fsar = firs ( 4 . 4 . 8 ) 

3 3 3 

from (4.4.7) we have 

t m - t sor = \8 n . (4.4.9) 

,333 

A 

Here S r * should be interpreted as the Newtonian internal angular 
momentum of the J.’th body since it is equal to 
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/ [(« r -f r )£ M -(tf a -|S3: 0r ]dx. (4.4.10) 

A 3 3 

a 

In the fourth order and for a = m equations (4.4.2) give the New¬ 
tonian equations of motion. In the fourth order and for a =m 
equations (4.4.3) and (4.4.4) give: 

t™ = t (4.4.11) 

4 2 

= 0 (4.4.12) 

3 1 

and 

t rai = i(S ra t\ 0 + S rt f\) . (4.4.13) 

4 3 3 

Equations (4.4.12) mean that the Newtonian internal angular 
momenta are conserved, as should be expected. We may summarize 

A 

the results obtained thus far as follows: the values of the in the 
approximation considered are unaffected by the rotation and are 

A A . 

the same as for non-rotating bodies; the quantities t™ and t r * are 

A 

expressed in terms of constant parameters ^"representing rotation. 

Let us go on to the next approximation. In the fifth and sixth 
order equations (4.4.2) define 

(t 00 £ W |o)| 0 —t 0 ” l '| 0 = (/« f OT |o)|o — t<m \o (4.4.14) 

—s— 1 5 1 —2—' 6 

through the terms of the order 2, 3 and 4. But the expression 
(4.4.14) would be exactly zero only in the Newtonian case in which 
jom _ However, in the case of rotating bodies, 
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this equals a spiaJl quantity: 

2 0m = (t°° t‘ * 0 ) 10 . ~ % (4.4.15) 

6 6 '~6~' 

which we can calculate explicitly from (4.4.2), though this would 
he quite troublesome. Going further, for a = 0 we can write the 
equation (4.4.3), after differentiating it with respect to time, in the 
form 

A 

H.+ « 0r +2*'”{.o} |0 = 0, (4.4.16) 

6 2 1 

which gives 

A A ^ t 

t m = S T8 £ s + f(f or -f #o | r |0 )d«. (4.4.17) 

4 8 <0 5 — r 1. 

Therefore, f r0 ° is determined by equations (4.4.2) and (4.4.3.) when 

4 

a = 0. Then we may use equations (4.4.3) for a = to as determining 
t mr if t m0 is arbitrary. Equation (4.4.4) determines «"**, and so on. 

« 5 6 

We also observe that condition (4.1.21) is satisfied since additional 
terms in the gravitational field are of the order r~ s . The additional 

term in A „ is 
» 

•«r A A 

■Kor = -2 V S^r" 1 ),.. (4.4.18) 

3 aTi 

It follows that the post-Newtonian Lagrangian for rotating bodies 
does exist and has the form 

/(— &jdx. (4.4.19) 


In this Lagrangian we are looking for new terms connected with 
the rotation. These are: 
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•4 ^£1^ .4 ^ 

^ Of 4“ i ^ 00 T&oojr 4“ ^ f0S fe oa|f 4™ \ t r8t Tl s $| f ) 4" 

3 k 4 2 '3 3 4 2 


32sr J 

•N .4 .4 .4 .4 

— — i (t 01, h Q r + t r °° ^O0(r + t r ° 8 A o«|r 4" t rat h^ r ) . (4.4.20) 

2“J 33 2 3 [3 * 2 

Neglecting terms proportional to the squares of angular momenta 
■we finally obtain, for a system of two rotating bodies: 


12 1 /l\ 12 2 /]\ 

2 ^^ , |o(-),2 8 —2^-r„( 7 ) i i g + 

2 1 2 /I \ 2 1 1 /I \ 

-wr,„ j-)i, +2^-r 10 ( 7 )^ + 


(t tr —t 00 £ r lo )dt (4.4.21) 
s —:— 


which we shall consider in the next chapter. 



CHAPTER V 

THE ONE AND TWO PARTICLE PROBLEMS 

1. ON THE QUESTION OF MEASUEEMENT 

We shall specialize our Lagrangian to the case of one particle, 
then to two particles and we shall find the motion in the chosen 
coordinate system, characterized by the approximation method. 
But before doing so we must know the physical meaning of our 
calculations. As we well know, the one-body pibblem — that is 
that of a small planet travelling around the sun — can be solved 
rigorously without too great difficulty. The result is that the planet 
moves along an ellipse, as in Newtonian theory, and that the ellipse 
very slowly rotates. But what do these statements mean? If they 
refer to a special coordinate system, then there is little physical 
meaning in such statements if we do not characterize the coordinate 
system within which they are valid. Yet, if we have to characterize 
the special coordinate system, then we violate the spirit of General 
Relativity Theory; its statements should be valid in any (or almost 
any) coordinate system. Some authors consider that the only way 
out of this dilemma is to introduce certain definite coordinate 
systems. This seems to us (as we mentioned in Chapter I) essentially 
a step backwards from the ideas introduced by G. R. T. 

Yet there is a simple way out of these difficulties. We have 
merely to admit coordinate systems which are Galilean at infinity. 
That is, coordinate systems in which for r -> oo the metric form is 

ds i = rj^dx'dctf. (5.1.1) 


There, we can speak about ideal rigid rods and about ideal clocks. 
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But physical events do not happen there. They happen in the 
Riemannian space with arbitrary coordinate systems satisfying only 
the condition that they go over into Galilean coordinate systems 
at infinity. The Galilean system is inertial and, of course, there 
exists a family of such systems whose members are connected by 
Lorentz transformations. However, we shall speak only about one 
member of the family. In the case of the sun, it will be the Galilean 
coordinate system in which the sun is at rest; in the case of a two- 
body system, it will be the coordinate system in which the center 
of gravity defined by (3.3.4) or one of the two bodies is at rest. 

How can we connect the events happening in our arbitrary 
coordinate system with an event in the Galilean coordinate system? 
Let us call r and £* the time and space coordinates of a particular 
event. An observer in the Galilean coordinate system learns about 
the event through a light-ray sent to him and reaching him at the 
point S k and at the moment T. In this way, every event in the 
Riemannian space is projected into the Galilean coordinate system 
very far from the event, where the influence of the gravitational 
field can be neglected. Let us assume that there are a number of 
events taking place around the sun at rest, all of them lying in 
a plane. (For the moment we omit the difficulty connected with 
the meaning of defining events lying in one plane.) Let us call this 
plane the (a?, y) plane. From each event a light-ray signal is sent 
so that it reaches the (X, Y) plane in a Galilean coordinate system. 
The light-rays are so sent that they all fall perpendicular to the ( X , Y) 
plane which is parallel to the (®, y) plane. (Again we omit the dif¬ 
ficulty connected with the meaning of parallel planes.) The descrip¬ 
tion of the pattern that the r, f l , f® events form on the (a?, y) plane 
depends on the choice of the coordinate system. The pattern of 
the corresponding signals T, S 1 , S 2 on the Galilean (X, Y) plane 
permits an objective description, since we know the construction 
of the rods and clocks by which the times and positions of these 
events can be ascertained. Thus our task is to find the cbrrespond- 
ence between the r, f* and the T, S h events, in order to be able to 
pass from the subjective to the objective description of the events. 
If we change the coordinate systems internally, keeping the Galilean 
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unchanged, then r, f* changes, but T, S k does not. To give a me¬ 
chanical example, imagine a rigid plane from which wires go out, 
ending in a rubber plane. We may change and deform the rubber 
plane, change the position of the wires, as long as we keep the upper 
.plane firm and the ending of the wires in it. The comparison with 
our example is obvious: the rigid upper plane stands for our Galilean 
coordinate system. The rubber is our coordinate system in the 
Biemannian manifold; the wires are the light rays. 

We shall therefore start our mathematical analysis by finding 
the connection between the r, event and the corresponding 
T, S k event. Since all test particles, independent of their mass, 
move along the geodesic, so also does the light-ray. This is partially 
an assumption, but a well-founded one. It is a null geodesic along 
which the light-ray moves. Since ds = 0 and there is no place for 
the concept of mass in it, it will be best to adopt a form of the equa¬ 
tions of motion in which ds does not appear; a form in which the 
time t is the parameter. For the sake of brevity from now on we shall 
assume: 

.. dx k fr 

c = l, k = 1, os 0 = t, «r = ---== ar |0 . 

at 

The equation of a geodesic line with t as a parameter is: 

(iiP)+A { ^ } x a if = 0. (5.1.2) 

These are three equations which will allow us to find the three ®’s; 
the fourth one which determines A is: 

gjifo/* = 0 (5.1.3) 

which states that the geodesic line is a null line. 

Let us start by completely neglecting the gravitational field. 
Then we have: 

?+i^ = 0, i k i k = l. 

A 


(5.1.4) 
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From this follows: 

So, x k = P+n k {t-r), n k n k = 1. (5.1.5) 

A 

This is the -world-line of a light-ray that was emitted from a point £ k 
at the moment r in the n k direction, in a Euclidean space and in 
a Galilean coordinate system. Indeed, we have: 

P = n k , P = 0, P — S k for t = t. (6.1.6) 

Now we wish to take into account the gravitational field. Obviously 
our approximation procedure, based on the smallness of the velocity 
of the moving bodies compared to the velocity of light is successful 
only if such a velocity is smaller than that of light. Obviously we 
cannot use this method for light-rays. Therefore, for our next approx¬ 
imation we shall assume: if P is the point on the world-line of the 
light-ray then its order is not changed by taking its time derivative. 
Therefore for the world-line of the light-ray we have: 

x k = x k +P+x k +..., 

0 1 2 

P — x k +P+x k +..., . (5.1.7) 

0 12 

P = P+P+P+.... 

0 12 

All our other assumptions remain unchanged. This means: the 
differentiation with respect to time of all other quantities increases 
the order by one. Therefore, in the next approximation, since 

fl>00 = 9 00 1 = Vi hmn = ffmnP 4m = <4n nV (5.1.8) 

2 2 

we have from (5.1.2) and (5.1.3): 

® fc +(logA)| 0 ® fc = <P\m® m P—<p\k, (6.1.9a) 

2 2 0 0 0 

PP S3 <p. 

2 0 


(5.1.9b) 
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From this we see that x k = 0, and equations (5.1.9) allow ns to 

I 

find the world-line of the light-ray in the aoxt approximation. 

We shall assume, without much loss of generality, that we are 
dealing only with a one-body problem, that is: 


<p = -^-, r 2 = ^x k . 
r 


( 8 . 1 . 10 ) 


The extension of the calculation to the many-body problem would 
cause very little change, since in this approximation the theory is 
a linear one and the motion of the bodies would affect only the light- 
rays in the next approximation, in which wo are not interested. 

We regard as our zero approximation that referring to the 
absence of the gravitational field, that is the one in (5.1.5). 

Multiplying (5.1.9a) by x k , we obtain: 

0 

(toggle = -9 = (5.1.11) 

0 

and because of (5.1.9a) the equation for x k is: 

2 

= 2<P\ m n m n k -<r ik . (5.1.12) 

2 o 0 

We shall now make a special assumption about the vector n k . 
We shall assume 0 


**»* = 0 (5.1.13) 

that is, that in the zeroth approximation the light-ray is sent 
perpendicular to the vector pointing from the sun to the point from 
which the light-ray is sent. We make this assumption for two 
reasons: firstly, from the point of view of observation this is the 
only case in which we are really interested; secondly, although the 
genera case in which (5.1.13) is not assumed could be solved almost 
as easily, still the formulas are much more clumsy. 
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Thus out simple problem now is to solve (5.1.12), assuming 
(5.1.13). Instead of (5.1.12) we write, because of (5.1.10): 




2 



x*n 8 n k 
o o 



(5.1.14) 


Since p is of the order two, we can introduce into this equation 
the expressions (5.1.5) for x k that is, the Galilean solutions. We 
shall now write, because of (5.1.13): 


r* = «*+(*_ T )*, £*!* = «* 

and for (5.1.14) we have a simple differential equation: 


~k 


2p 


(5.1.15) 


(5.1.16) 


(a a +(«-r)f/* T) ?* 

By integration of the last equation wo find: 

■k_ in k _ 2 / * , £*(«-*) \ 

i s (o® + (t —r) 2 ) 1/a \o + « a / 

where b k is the arbitrary constant to be determined by some added 
s 

conditions. As such we wish to assume: 


2 


(5.1.17) 


n k = 0 for t = oo. (6.1.18) 

2 


This means: far from the gravitational field the light-ray has the 
prescribed direction, which is perpendicular to £*. That is, for t = oo: 


therefore 


n 

2 





(5.1.19) 



2^ 

(<**+(<— t )*) 1/1 yt 


f*«-T) 

a® 



(5.1.20) 


For t = r we have: 
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Therefore n k +n k must have a small component in the | fc direction 

•os 

at the time t = r, in order to be perpendicular to the £* direction 
at time t = oo. This is the famous formula for the deflection of 
light-rays. For t = — oo it follows from (5.1.20) that the component 
in the £ k direction would be 

= ip. (5.1.22) 

2 

Now we wish to integrate .(5.1.17) again, under the condition 

*V T = o ( 5 . 1 . 23 ) 

2 

because £* is the exact point on the light-ray for t = t. We obtain: 

s o a a 2 a \ a j 

r 2 = l 8 f i +(i-r) 2 = a 2 +($-r) 2 . (5.1.24a) 

We see that this integral of (5.1.20) does indeed satisfy the conditions 
(5.1.23). Now for t -► oo we have: 

x k {t) -2[M k \og —- ■ ~ T ^ . (5.1.24b) 

2 d 0 d 

On the right-hand side of this formula we have omitted all ex¬ 
pressions which go to 0 as t -> oo. We are essentially interested in 
the projection of the events £ k , r, into the far-away Galilean system, 
that is for t-*- oo. Let us now, as in our introductory remarks to 
this section, call 

S k = x k +a? for t —^ oo. (5.1.25) 

0 2 

Then, because of (6.1.24b) and (5.1.5), we have: 

S k = £* (l + + n k (t- r) - 2/Mi fc log — *~ T ^ (5.1.26) 

\ d J o o d 

and 

S k = n k . 


( 5 . 1 . 27 ) 
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The formulas (5.1.26) and (5.1.27) are valid "with the same accuracy 
as (5.1.24b). They are the result of our calculations connecting the 
events in our coordinate (Newtonian) system defined by our approx¬ 
imation procedure with the events in the far-removed Galilean 
coordinate system. But we should remember that these formulas 

were obtained only under the assumption that n k is perpendicular 

0 

to We shall use these formulas when we try to interpret the 
motion of a particle in an objective way. In any case we see that, 
generally, because pja is small, the approximation 

S k = £*+»*(* -r) (5.1.28) 

o 

is a good one. Or, with a great degree of accuracy, we may inter¬ 
pret the events in our Riemannian coordinate system in a Galilean 
way, as long as they are not too near the sun. For example, if we 
consider only the earth and the sun with mass M then for an event 
on the earth we have: 


a «* 1.5-10®km, p ■■ 

TiM 

= - -ft* 1.5 km. 

o* 

(5.1.29) 

Therefore 



2p 

- ph 2 • 

a 

10 -8 . 

(5.1.39) 

As an example, let us consider an 

event for which 


c* 

II 

o 

h-* 

II 

f* = I s = o, 

(5.1.31) 

and since £*«* == 0: 



0 

= n a = 0 , 
0 0 

n* = 1. 

0 

(5.1.32) 

Then we have (compare fig. 1): 



iS 1 = $-f“ 2p , 3* = 0, 

. (5.1.33) 


Or, for a sufficiently great t: 


S 1 = a+ 2j«, 5* = 0, 3* = t. 


(5.1.34) 
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We shall return to the contents of this section in order to interpret 
the motion of particles objectively. 


-|i:-a+3u 


SUN 


/ 

J 

<'-a 


Fig. 1. 

Connection between a subjective and objective representation of an event 


2. ON THE MOTION OF A TEST PARTICLE IN THE FIELD 
OF A HEAVY PARTICLE 

We know that the problem of a test particle moving in the field 
of the son can he solved rigorously; however, using our general 
formulas here, we shall find the motion in the post-Newtonian 
approximation without making use of the explicit metric form 
or the equation of the geodesic line. 

We start from the Lagrangian and assume: 

£* = 0 , £* = !*, h — Hi n = Aft. ( 5 . 2 . 1 ) 

We shall consistently neglect all products of A/i . Then, dividing 
the old Lagrangian by Ap, for the new Lagrangian L we have: 

£+£= i£V 0 ,o+7-+M (5 - 2 * 2 > 

2 4 T AT r* 

It is now more convenient to introduce instead of a?° a new variable 
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so that 

Ldx 9 = L' dx' 9 . (5.2.4) 

Denoting the derivatives with respect to x'° by O', we have: 


•7+^ *V«V-.< s - 2 - s > 


l '+ l ' = - 10 >l° |o'+ 

2 4 & 

The expression in the bracket in the above equation is constant 
because of the Newtonian equations of motion. Ignoring it, we obtain: 


L'+L' = *fy«V+ f 


H 3/a* 


(5.2.6) 


This Lagrangian corresponds to that in classical mechanics with 
the potential function 

3^ a 
r r a 


V = —— 


(5.2.7) 


We have as its integrals the law of conservation of momenta, 
which in the polar coordinate system (<p, r) are: 


r a 9>i 0 / = J ~\~J 

1 3 

and the law of conservation of energy: 


= E+E. 


(5.2.8) 


(5.2.9) 


The problem is therefore reduced to the one well-known in 
classical mechanics. Its solution, written in a parametric form, is: 




1+ (T=7ijs) 


r = <*(1 — scosm), 

^ /l *4" s u 

2arctan 1/ --tan — , 

r 1 —e 2 


(5.2.10a) 

(5.2.10b) 


,x 9 = I 1 + (« — esina). (5.2.10c) 


Motion and relativity 
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The constants a and s are determined by the constants E and J. 
However, it is more convenient to regard a and e instead of E and J 
as given in advance. That (5.2.10) are solutions of the equations 
(5.2.8), (5.2.9) can he checked by substitution. In the case e > 1, 
we have to replace 

a by —a and u by — in. (5.2.11) 

We now wish, after the integration of our problem is finished, 
to go back to our original time ®°: 


= 



dx' a 


- 4 ’’’ + «(i -Ls.) V 7 ( 1+ lr) 


After integrating and putting for a>'° its effective value (5.2.10), 
we have: 

*Vt[( i+ £M i+ £H- (n - 3 - w) 

This result should be substituted instead of (5.2.10c). 

In the Newtonian approximation we have: 


r — «(1 —ecos«), 

<P — 2arctan|/iiLf tan-^, 

= ]/~ («—esinu). 


(5.2.14a) 

(5.2.14b) 

(5.2.14c) 


• ^ uestio11 considered in almost every text-book on G B T 
is that of the path of such a particle. This path and the perihelion 
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motion can be easily deduced from our moTe general formulas. 
We introduce an auxiliary angle ip: 


~ ,/l + e v> 

w = 2 arc tan 1 / -tan—. 

X 1 —e 2 

It can easily be shown that 


1 — ecosw = 


l-«* 


l+scosy 

Therefore from (5.2.10) we find: 

# 

a(l —« 2 ) 


r = 




(5.2.15) 


(5.2.16) 


(5.2.17) 


After one rotation, the perihelion has moved by an angle cr, so that 


which gives the famous perihelion motion: 

° = «(1 — e 2 ) ' 


(5.2.18) 


(5.2.19) 


The question now is: what is the physical meaning of all these 
results! Or, how will they differ in the Galilean coordinate system! 

Suppose that the observers at infinity are scattered in a plane 
parallel to that of the motion. They receive light-signals on 
a plane for which. Z = x* is very great, but constant. Let us call 
(®, y) the plane in which the motion proceeds and (X, Y) a far-away 
plane, but parallel to the (at, y) plane. 

If the light-rays fall perpendicular to the (X, Y) plane then 
putting X for t and denoting 


6 * = ft 


(5.2.20) 
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from (5.1.26) we get the connection between the |'s and the S'&: 

& = 5*= f*(l + ^-), (5.2.21a) 

F® = T-T-2^1og 2(I, ~ T) . (5.2.21h) 

Q 

Here |*(r) (& = 1,2) is the motion in the (®, 2 /) plane. This motion 
is observed on the plane 3 s = const. We ask about the connection 
between T and r, or rather about the connection between 

AT = T i —T 1 and Ax = x i —T l (5.2.22) 

referring to two events. We wish to find the connection between the 
objective rhythm of events denoted by T, and the subjective rhythm 
connected with the choice of the coordinate system. Since S 3 = 
= const we conclude from (5.2.21b) that 

AT-Ax- + ^ At = 0 . (5.2.23) 

T—x q dx 

Since T—x 00 we have: 

AT = Ax(l + ^^. (5.2.24) 

For our solar system 2/t/g is of the order 10 -a . The velocity dQ/dx 
is very small compared with that of light (remember the choice 
e = 1). Therefore, neglecting expressions of the third order, we 
can just write: 

AT = Ax. (5.2.25) 

We see, therefore, that for all practical purposes the time has the 
same rhythm on the objective plane as on the (a?, y) plane where 
the motion takes place. 

How, consider the r,<p coordinates. We have: 
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Therefore 

jja £* 

— = tan<P = — = tan<p. (5.2.26) 

Of the three elements describing the motion r,<p,x 0 two are the 
same in the (a>, y) as in the (X, 7) plane. What about r? Denoting 
the projection of r on the (X, Y) plane by R, we have: 

R = r|l + -~-J = r + 2/x. (5.2.27) 

Thus the change in r is a trivial one. It consists in adding to r the 
constant value 2/i. 

Therefore we see the great virtue of the accepted coordinate 
system in Biemannian space. Only by adding a small but constant 
value to r and by adding a large but constant value to a>°, can we 
interpret all our results objectively as those of observers in 
a Galilean coordinate system who receive the signals of the events 
through light-rays falling perpendicularly onto a far-removed plane, 
parallel to the one in which the events take place. 


3. THE TWO-BODY PROBLEM 

We know that in Newtonian mechanics the solution of the two- 
body problem is just as easy or as difficult as that of the one-body. 
In G. B. T. we can also reduce the two-body to the one-body 
problem, although it is more difficult to do so than in the case of 
Newtonian dynamics. This statement applies only to the mechanics 
of finding tho solution. It does not apply to the mechanics of finding 
the Lagrangian; here an essential difference existed between finding 
it for the one- and for the many-body problem. 

Thus, as in the case of the one-body problem, we start from 
our General Lagrangian which we now apply for two particles. 
Writing 

12 1 2 
r 2 -£*)(£*-£*) 


(5.3.1) 
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we obtain the following expression: 



1 2 



(f'V f*|o + f' 8 |o fV) H~ 



hho+ikhho) i +ikhho)+ 




(5.3.2) 


We shall choose an. especially simple coordinate system in which 
to solve our problem. It will he the system in which the Newtonian 
center of mass is at rest. This means that we assume: 

d+ft)X k = = 0 (A 4 ). -(5.3.3) 

2 2 2 0 2 0 

Or: X k , developed into power series, will start with at least the 
order two. Therefore the product of Xf 0 Xf 0 by the mass will be of 
order at least eight, and can be neglected. 

We now introduce a new variable and new constants: 

1 2 

£*— £* = v’t Po = /*+/»> ft = -^-- (5.3.4) 

Po 

This, together with (5.3.3), gives: 

«*—>=-4^- (5.3.5) 

ft ft 

We can now express L as a function of r\ - and its first derivatives 
with respect to time. Using 

r * = fe-hCf-h 


(5.3.6) 
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to obtain Liu, whin* ^ 

IM ’ mct we again denote by L; 

Z ‘ + ‘ ~ Hr+i ( 1 ~%) < ’’’ri‘+t<3/‘,+p)j-iy + 


Mo 

' r s r * ' 


(5.3.7) 

This form of L s«m, 0 + . 

ta t t, OM ' 1> ^ CaM ' » 0®°be e rZSd”o ttT Pli0ated «■»* 

hut we must make an anuronrift * * the form for one particle 
“ We sM *■ * 


== d$'° ^ 

-4 + *)« 


(5.3.8a) 


(5.3.8b) 


errz^srjss 




= j^o 


“« tortToSS “” lapl5r tte by an arbiter, wmt J 

f+f - i->v,v + f +i f hv , v _^j + (MJ1) 

^ lorn can be acM.vcd cm, mte certain J^on. „ tte 
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constants o^/S ,s,p. These are: 



/*o 


V. 3 


(5.3.12) 


Since the expression in brackets in (5.3.11) is constant because 
the equations of motion in the previous approximation are satisfied, 
we may write: 


3/4 


+^+- 

' 1 

1-3—)r' 2 


l fij 


(5.3.13) 


Thus the problem is identical with the one-body problem under 
the influence of a field with a potential function: 


V = 


<“o 

r' 



(5.3.14) 


Therefore the rest of the argument is almost identical with that 
in the previous section. For the path we find: 


r = 


a'(l-e a ) 


1 + ecos 

But from (5.3.8b) we find: 

r = ar'+a8fi 


3/*o 


\ i“o/ 


* 2 ) 


- (- 3 ' 


+£/<• 


(5.3.15) 


(5.3.16) 
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Therefore, introducing 


we have: 


r = 


•'('-SI- 

a(l — s 2 ) 


i+, -[( 2 -s<&iH 


+f 
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(5.3.17) 


(5.3.18) 


Ignoring the additional constant —, we find that for s < 1 the 

2 

motion is along an ellipse and that after one rotation the perihelion 
moves by an angle 


6^|t{ 0 

^ ~ a(l —e 2 )’ 


(5.3.19) 


which is identical with the case of the one-body problem, the only 
difference being that now ju 0 denotes the sum of the masses. 

What is the physical meaning of these results? To answer this 
question we must go back to the first section of this chapter. At 

i 

the point and the moment r a light-ray is emitted. But now the 
two particles move. Yet the contribution coming from their motion 
can be neglected, since it is of higher order than two. The other 
difficulty is that the mass emitting the light-ray is now finite. This, 
too, will be without influence on the light-ray, because of the 
spherical symmetry of the emitting particle. Also the fact that 

the potential is infinite at the point has to be ignored, since the 
assumption about the existence and character of the singularities 
is of a purely mathematical character. Indeed, the tweedling process 
was invented and introduced in order to remove this kind of 
difficulty from the analytical treatment. 

Therefore we may treat the problem of a light-ray emitted at 

i 

the point £®, at the moment r, in the same way as we treated the 
problem in the first and second sections. 
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1 2 

Let us slightly generalize the transformation (5.3.4) from f, f, 

to X‘, rf, hy assuming that the Newtonian center of gravity moves 

2 

so that at any moment in our coordinate system we have = 0. 
Indeed, such a generalized transformation is: 


t = X 8 +-~-rj s , 


and the condition 



Therefore 


2 

I 4 = 0 means 



(5.3.20) 


(5.3.21) 


+ (- r >-3.22) 

' n fi 

This shows the physical meaning of rf-. it is the coordinate of k 
in a system in which the second body is at rest. Therefore the motion 
of the first body relative to the second or vice versa is described 
by (5.3.18). The light-rays sent by the two bodies are slightly curved 
and reach a plane parallel to that of the motion. Those from one 
end, curved by the presence of the particle at the other, are shown 
in diagram 2. Therefore, if we denote by B the objective distance 
measured in the infinitely far removed system we have: 


R = 


«( 1 - 8 S ) 




+ -J- + 2 /“ 0- 


(5^3.23) 


These are trivial changes in formula (5.3.18) from which the 
additional /i/2 can be left out as well. Also, as follows from the 
end of the previous section, the rhythm of time is almost the same 
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in the Galilean as in our Riemannian coordinate system. Once more 
we see the virtue of the chosen coordinate system. It gives us almost 
an objective picture of the motion. 



2 1 
u r u 


Fig. 2 

Distance R judged by an observer in a G-alilean coordinate system 


4. THE MOTION OF ROTATING BODIES 


In the last section of Chapter IV the Lagrangian of a system 
of two rotating bodies was found. This Lagrangian contained new 
terms connected with rotation. We already know the corrections 
to the motion of spherically symmetric non-rotating bodies arising 
in the post-Newtonian approximation. In the present section, we 
shall derive the corrections due to rotation, using the Lagrangian 


i i 


2 2 




1 2 
HH 




(5.4.1) 


in which the only post-Newtonian terms are those accounting for 
rotation and in which the g’s from (4.4.21) are neglected. 
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Using the following notation: 

1 2 

11 2 2 

H 0 2. r = fi£ r +/uf, 

1 2 

itt 0 = p'+p, (5.4.2) 

1 2 
flfi 

f* ~~1 2 ’ 

M.+i u 

pS" = p8 n +pB n , 

we ean rewrite the Lagrangian (5.4.1) in the following form: 

L = torr+p[\j'rf+ - -isSy)• (5.4.3) 

Since we are considering motion in the system of the Newtonian 
center of mass, if we omit the constant factor, the following 
Lagrangian of relative motion is obtained: 

This Lagrangian is invariant with respect to the time displacement 
transformation (3.4.2) and the rotation 

v r = v+jry, 

$"* = 8 r *+M rt 8 te -M st S tr (5.4.5) 

with -M 1 * a constant bivector. Noether’s theorem gives: 

\rfrf— — = J( = const 
r 


(6.4.6) 
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and 

(5.4.7) 

where 

JT * = ’JV“ *f*f+ ^sr [8**1*V s - 8 t, *i t rj r ). (5.4.8) 


We shall now consider two typical particular cases. First, we 
assume that the motion is bounded and takes place in the pi a.™ a 
perpendicular to the biyector S ra . Written out in vector form, equa¬ 
tions (5.4.7) and (5.4.8) are: 

2 

J=-jjSxJ, (5.4.9) 

T 2 

J = yXv- — [Sxy)xy, (5.4.10) 

T 

which, in the case under consideration, gives: 

J ra = *f vf — rj*rf+ —S'* = const, (5.4.11) 

r 

or 

2 

J = -S= const. (5.4.12) 

r 


Introducing polar coordinates r,<p in the plane of motion we have: 


J — 


2ji 0 8 

r 


(5.4.13) 


where J and 8 are the magnitudes of J and S respectively and 
the sign + is for J and S pointing in the same direction — 
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the opposite. The right-hand side of the identity 

( du \ 2 . • • l 

—J +« =i?V(rV)- s , u = T 

can now be expressed as a polynomial in the variable 


V.4 


(5.4.14) 


(d^r) ' JrU * ~ ^ 2 ^ 2 -®+| 2 i“o+8 —+ (5.4.15) 

The solution of the latter equation in the Newtonian approximation 
is a Kepler ellipse 

jp / 2EJ 2 \ 1 I 2 

u=p [l+«cos( 9 J-y 0 )], P=—, e = |i+__.j . (5.4.16) 

The next approximation yields a rotating ellipse with the very 
small displacement of the periastron 


±8ji 


/“off 

pj 



(5.4.17) 


for each period of the Newtonian motion. 

Another case is that of Newtonian motion around a circular 
orbit r = a. Equation (5.4.7) now gives: 

2 

J=-^-SxJ. (5.4.18) 

Thus, the vector J, which is perpendicular to the plane of motion, 
processes around the constant vector S with angular velocity 



(5.4.19) 


which we can see immediately if we place the S vector in, say, 
the ® s direction and let the J vector process aro und it. 

All the considerations concerning the question of measurement 
axe also valid in the case of rotating bodies since the path of a light- 
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MOTION AUD RADIATION 


1. A SIMPLE EXAMPLE 

Tn this chapter we intend to show the close connection between 
the equations of motion and gravitational radiation. To do this 
we shall start with a simple, almost naive, bnt instructive example. 
Let ns assume (as we did in the pre-General Relativity era) that 
we are looking for. a field equation of the gravitational field which 
is consistent with Special Relativity Theory. The simplest idea 
would be to postulate a scalar D’Alembertian equation, that is one 
of the form 

N A A 

(6.L.1) 

A -»1 

Sere, we assume the mass to be a function of time. We also assume 

AAA 

m = (tt. 1.2) 

2 4 

A 

where m is constant and throughout the entire chapter we take the 
2 

velocity of light as one. To find the radiation emitted from the 
sources of the field which, of course, move arbitrarily, we must first 
find the retarded solution. We know that such a solution is: 


»w»,o - - 2;jv i, (MJ) 
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And for the advanced potential: 

w*.*> - -gf*h ±w-fo+*>i . t6 . a . 4) 

In the last two equations we have: 

B' 2 = (x ,a -x a ){x' a -x a ). 

Finally we have the standing wave solution: 

^ST = iK9 , BBT + 9 !, ADv)* (6.1.5) 

Obviously the radiation is connected with the retarded potential 
solution to which we shall devote our special attention. 

We develop the m ?s and the <5’s into power series: 

= m(t)-m{t)R'+$m{t)R' 2 +... (6.1.6a) 

S(x'-i(t-B')) = <5 (x - f(<))- V# + i # |00 B' 2 -i d mo R '*+... (6.1.6b) 
This gives, up to the filth order in <p EET : 

J ’ * A A A A 

[m(X -1 H-ifipo — l ;R 2 im )—m] (6.1.7) 

1 * 5 

where 

B 2 = {x 3 -f*){x a -i‘). (6.1.8) 

This is the usual approach to the field equations (6.1.1) through 
the retarded potential (6.1.3). But it will be more instructive to 
approach the problem from the point of view of our approximation 
method. 

To simplify the problem slightly let us assume for the moment 
that m = m(t). In other words, we shall not develop m into a power 

series. Then, according to our approximation method we have: 

A A 

(p\sa = 4 7t 2jind. (6.1.9) 


Motion and relativity 
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ray depends only on fc 00 and h„ and rotation does not affect these 
2 2 

quantities. 

We wish to add some remarks concerning the neglect of the g®”* 
defined by (4.4.15). Of course, they are of the same order as some 
other expressions in L. Therefore, they should not be neglected. 
But the situation is this: we have the equations 

J &*.jdx = 0 (5.4.20) 

which for 

( 5 - 4 . 21 ) 

Awl 

determine g® 1 " and f°® in a complicated way. The motion is still 

< 4 

arbitrary and determined by our equation 

djLdt = 0 (5.4.22) 

in which variation has to be performed with respect to f*. But, 
to determine X, we must know g®”* and this requires the solution 
of an integral-differential equation. We can, however, find a parti¬ 
cular solution of (5.4.20) and (5.4.22). It is: 



In this case, the entire reasoning so far presented remains almost 
unchanged. Only the definition of S'* in the last equation of (5.4.2) 
changes. This means that for a coordinate system in which the center 
of gravity is at rest, we have motion of the same character as 
considered in this section. 
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Therefore: 

A 

V! m 

X B 

(6.1.10) 

From fhis we find 

§ 

rt 

1 

i 

ii 
$ 

. (6.1.11a) 

and 

•T-, A A 

<p = -£ >j imR )\»o- 
* JL 

(6.1.11b) 

Or, if the order of 

<p jumps by two, we find: 


fST 

A ?= 0 ' ' 

(6.1.12) 


But we see immediately that this is tlie standing wave solution, 
since it is insensitive to the transformation t' = —t. Thus, the 
simplest solution by our approximation method corresponds to tho 
standing wave solution. To obtain the retarded potential solution 
we must start with <p and then again jump by two steps. "We must 

3 

start with a harmonic function and (6.1.7) tells ns what kind of 
harmonic function to take. We must assume 

<P = m. (6.3.13) 

3 3 

Then, again proceeding by steps, we find the radiation part of <p: 

n n 1 d a+1 * x , 

Vbad = 2 2 (21+1)1 d ^ 1 {mR )- (6 - 1,14) 

Thus the approximation method gives us a natural division into 
the standing wave solution and that representing radiation. From 
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this we have: 

9ret = ^st + ^ead? (6.1.15) 

^ADV ^ ^ST ( P'RaX> 9 (6.1.16) 

We are returning now to our problem of finding the radiation 
with the help of the retarded potentials as expressed — up to the 
fifth order — by (6.1.7). 

The energy-momentum tensor for such a field is: 

4nT„ fi =<p^-id^ W ' 9 . (6.1.17) 

Indeed, we find: 

inTj’i,, = <p\/<p ]a +<p' p <p ]afi —<p [t <P\ ta = — in JTmdgv (6.1.18) 

A 

Therefore, outside the singularities we always have: 

Tf v -0. (6.1.19) 

As usual, we define the flux of radiation by 


/ T a ‘n s d8 = P a ■ (6.1.20) 

£ 


where 2 is a very great sphere surrounding all the singularities. 
Wo wish to find the energy-momentum flux, up to the eighth order, 
in the case of our field equations. 

Wo return to our expression (6.1.7). We are interested in it 
only for very great values of the 2?’s. Let us call r the distance from 
a fixed chosen point which is the origin of the coordinate system: 

r 2 = ®V. (6.1.21) 

Then 


df a i d*f al a l 


a a 

£* tV — 


1 d*f 


6 dr * 


A, A n A n 


( 6 . 1 . 22 ) 
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Since r does not depend on time, we have from (6.1.7) and the last 
equation: 


-1 


w> 


-\MX 8 x , +\D tt +M 


(6.1.23) 


where 

M=]?m, MX 8 =^mP, D” = £ mpf*. (6. 1.24) 


We find for fyn, only up to the order r~ 3 : 

(6.1.25) 

We shall now find 

4jrT 0 m = <p^ m = -<P\o<P\m- (6.1.26) 

But we are only interested in the expressions that give contributions 
to the surface integral such that they are of the order l/?* 2 and such 
that their surface integral oyer the infinite sphere does not vanish. 
In the lowest non-trivial order we obtain: 


P. = / Tfr^dS 

JP m = J T m *n a dS = $M*X m , [M = M). (6.1.27b) 


We see from the last equations how intimately the radiation 
is bound up with the motion. The radiation would have vanished 
if we had dealt with the standing wave potential. But even for 
the retarded potential the radiation can be made whatever we 
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wish hy the proper choice ol the motion of the center of gravity 
and by the proper choice of M. Let us assume (for example) that 


X m = 0 


(6.1.28) 


which is even less than assuming uniform motion for the center 
of gravity. Then we have: 


p 0 = - |if 


i a*D* 
6 dtf 


Pjfc — 0. 


The entire vector vanishes if 


(6.1.29) 


M = --D w +trf4-/3, (a,/? = const). (6.1.30) 

4 6 


.D m =(£ m P P) [90 = 4T+ 2 £ m £ 9 f », (6.1.31) 


T = l£m£ a i a . 


(6.1.32) 


In the case of a Newtonian force, for example, we have: 


A a 




ABA JB A A B A B A 

mmf (f — I ) 1 V 1 ' m>m£ a {£ a — f*) 


-t2 


(6.1.33) 


ABB B A AB 

1 (£•—£ ) 1 V 1 ' mm 


l vr 

i2j 


1 V 1 ' mm -rr ^ . 

— “— = F = potential energy. 

AfB ^ 


Therefore 


D* = 4T+2F = 2T+2JS7 


(6.1.34) 
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where E is the (constant) energy. 

Therefore, in the case of a Newtonian force, P 0 vanishes if 

Jf = -%T. (6.1.35) 

4 

We developed T m n np to the seventh order and T™ up to the 
eighth. But the general idea is clear. We can annihilate the radia¬ 
tion by a proper choice of the laws of motion. The situation is more 
complicated in the case of G. R. T. There the equations of motion 
are determined by the field equations. Yet the general idea of our 
investigation remains the same: to establish the connection between 
the equations of motion and those for radiation. There is little sense 
in considering radiation without the sources; but motion of these 
sources and therefore the gravitational radiation, too, are determined 
by the field equations. They are also determined by the chosen 
coordinate system, the chosen boundary conditions, and the chosen 
solutions. We shall show that for the coordinate condition and 
boundary conditions that seem reasonably simple, the radiation 
vanishes. 


2. THE EQUATIONS OP MOTION IN THE POEM OP A SUBFAOE 

INTEGRAL 

We shall recall here how the equations of motion in G. R. T. 
were originally formulated in 1938, or rather eleven years ’ later 
in 1949. Although the formalism here is much simplified and 
although we are interested only in the general form of these equa¬ 
tions, yet the idea remains the same: to express the equations of 
motion through surface integrals enclosing the sources. 

We recall the Einstein field equations written in the form 

( 2 . 8 . 12 ): 

+ = K imP ^+N(® ia ) + 87c‘.T' a ' = 0 . ( 6 . 2 . 1 ) 

Here 

i w* = K (6-2-2) 
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is the linear part of © M ', being antisymmetric in ju, a and v, p and 
insensitive to the interchange of pa with vp. The expression W(©'°') 
is the non-linear part of ©^ which we shall denote by A m and 
is the energy-momentum tensor density. The relation of the y’s 
to the jf’s was explained fully in Chapter II. We rewrite (6.2.1): 

K>‘ a ^ afl -\-A' a +toff* = o. (6.2.3) 

We shall now prove a lemma which will be used several times later. 

Lemma. We have a function 9 r(, " )fcs , skew-symmetric in the 
indices Tc , s and having arbitrary Greek or Latin indices besides, 
represented by the dots in the brackets. Then we have: 

ftl~ )k ° la n k dS = 0 (6.2.4) 

jj 

if S is an arbitrary closed two-dimensional surface not passing 
through the singularities of the field. By n k we understand: 

n k — cos(® fc , n) (6.2.5) 

that is, the components of the “normal unit” vector to the surface E. 
The words “normal” and “unit” are used here only in the conven¬ 
tional sense, to designate the corresponding functions of the 
coordinates, which are implied by these terms in Euclidean geometry. 

The proof is simple enough. First we see that the integral 
(6.2.4) is certainly independent of the shape of the surface as long 
as the number of singularities enclosed by the surface does not 
change. This is so, since 

7 ( ->*V = 0, (6.2.6) 

and because of Green’s Theorem. Let us now write: 

7<- )to = (6.2.7) 

' or explicitly: 

gr(...)2» ^ ^(..Osi = ^•••) > 9*->“ = ■4- ) . (6.2.8) 

Then we can write the integral (6.2.4) in the form: 

= f e k3r .^ r{g n k d8 = f enrl„ sZdS. 

£ £ 2 


(6.2.9) 
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But this integral can he changed by Stokes’ Theorem into a line 
integral over the rim of the surface. But the surface is, as we have 
assumed, a closed one. Therefore its rim is of zero length. Therefore 
our lemma is proved. 

We go back to the field equations. We shall rewrite (6.2.1) 
putting fi = m: 

® mr = yl 1 "*'. (6.2.10) 

Let us take any two-dimensional closed surface surrounding 

j. 

the .d'th singularity. We call such a surface I. Then, because outside 
the singularity Einstein’s tensor vanishes, we have, if we take our 
lemma into account: 


fS mo ’%n m d8+ f A^nndS = 0, A = 1 , ..., JV. (6.2.11) 

A. A 

£ £ 

These are 4 N equations valid (again because of our lemma) for 
arbitrary surfaces, each of them enclosing only one singularity. 
Because of this arbitrariness of the surface, they cannot give us 
any relations between the space coordinates of the field. They can 
give us only relations between the coordinates of the singularities 
and their time derivatives. Thus these are by definition the equa¬ 
tions of motion for the N particles. From them we can obtain the 
Newtonian and the post-Newtonian equations of motion. This can 
be done with the help of the “new” approximation method, although 
the application of the method is much more elaborate than that 
presented in this book. 

Indeed, to find the post-Newtonian equations of motion in this 
way we should have to calculate y m , y m , y mn and y m and then many 

2 3 4 $ 

surface integrals! 

But, as we said before, we regard (6.2.11) by definition as the 
rigorous equations of motion for N particles. 
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3. THE EQUATIONS OF MOTION IN THE FORM 
OF A SPACE INTEG-RAL 


Let us recall now the form of the 4 N equations of motion for N 
singularities known from Chapter I: 

fS^dx = 0, A =1,..., N. (6.3.1). 

A 

a 

As we remember, this leads to: 


Ch +{ ap} °> A-1,..., IT. (6.3.2) 


The transition from the first to the second of these equations is 
bound up with the expression for Z** which is: 


and 


= y p"d 


A A+i 

r-ppr. 


(6.3.3) 

(6.3.4) 


It can be shown that it is sufficient to assume (6.3.3). Indeed, the 

A 

specific form of tf", as expressed in (6.3.4), follows from the field 
equations. This means: it is sufficient to assume that the energy- 
momentum tensor is a linear function of the <5’s in order to specify 
its form. 

To show this, let us start from the expression: 

j0Z afi . J dx = O (6.3.5) 

A 

Q 

which follows from the Bianchi identities. Here 0 is an arbitrary 

A 

function continnons on* the world line f. We can write (6.3.5): 

A A 


A 

Q 


(6.3.6) 
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This leads to: 



A^e+A^e^ = o. 

(6.3.7) 

We shall now calculate A afi . Because 



1 

11 

(6.3.8) 

we have: 




= (-r+r°f) 0 18 = o. 

(6.3.9) 

Since |° = 1, this can be written: 



A*% = v = 0. 

(6.3.10) 

Because 0^ 

is arbitrary, we have: 



A* = (-r’+^V) =0. 

(6.3. IL) 

Putting a = 

‘ 0 we find 



03 * 

* Mil 
o 
o 

II 

(6.3.12) 

Therefore generally 



f* = t ## f a £' , = /*?? 

(6.3.13) 


which, is the proof of our theorem. Obviously A a = 0 gives the 
equations of motion, that is (6.3.2). 


4. THE EQUATIONS OP MOTION IN THE FORM OP STACK 
AND SURFACE INTEGRALS 

We shall introduce here the third and last definition of the equa¬ 
tions, of motion. 

Let us recall that the first definition was based upon integrals 
taken over the two-dimensional surfaces surrounding one singular¬ 
ity; that the second was based upon integrals taken over a volume 
surrounding one of the singularities. In the definition which we shall 
give here, both the surface and the volume integral will appear. 
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We start again from the equation 

©'■'+8 rtST" = 0 (6.4.1) 

where 

©/*■ = (6.4.2) 

Because of the skew-symmetry in the indices vfi, we have 

(®p-_|_87ESr^’')| X , = (A^+ 8n 9n w = 0. (6.4.3) 

Here, obviously, the containing only non-linear expressions 
in the y’s is not a tensor. One would be inclined to regard such an 
equation as the expression for the conservation law, involving the 
“pseudo” energy-momentum density A** for the gravitational field. 
For the moment, however, we shall ignore such an interpretation 
and we shall use our last equation to define the equations of motion 

A 

with its help. This we shall do by integrating it over Q. Changing 

A 

the integral over the space derivatives into a surface integral over Z, 
we obtain: 

J A am n m dS = - -1- J (/1°“+ 87^7 0a ) dx. (6.4.4) 

A A 

x a 

Theso are again 4tN differential equations which we define as 
the equations of motion for N particles. The energy-momentum 
tensor does not appear on the left-hand side of these equations 

A 

since it vanishes at the surface Z. 


5. THE CONSEQUENCES OF THE DIFFERENT FORMS OF THE 
EQUATIONS OF MOTION 

We introduced the equations of motion in three different forms. 
Firstly: the ( Z ) form, according to (6.2.11): 

-JiT” °' vp w n m dS = fA mv n m dtt. 

A A 

X X 


(6.5.1) 
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Secondly: the (13) form, according to (6.3.1) and (6.3.2): 


J 

A 

Q 

(6.5.2a) 


(6.5.2b) 

Thirdly: the {IQ) form, according to (6.4.4): 


/ ATn^dS = —f(A 0a +8nff 0a ) dx. 

(6.5.3) 

A A 

£ O 


We can easily show that these three forms are 
each other, assuming that the field equations 

equivalent to 

©'"+8jtST' w = 0 

(6.5.4) 

are satisfied and assuming also 


■it 

II 

& 

II 

(6.5.6) 

Let us first investigate the connection between the (13) and 
the (1713) form. The (13) form was derived from the symbolic equa¬ 
tion: 

^ = ^+{“} s: ' w = o, 

(6.5.6) 

whereas the {IQ) form derived from the symbolic equation: 

8 *9* V +A* V = 0. 

(6.5.7) 

From these two equations follows: 



(6.5.8) 


as a consequence of the field equations and the Bianchi identities. 
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Lot us now put into the (X£f) form of the equations of motion, 
instead of 9'°“, its explicit expression. Then we have: 

A 

4r^' n “ —£■[]>«»«+/>*]. ««•») 

■4 A. 

2 Q 

or, because of the ( £2) form: 

i ( 6 . 6 . 10 ) 

A A 

2 a 

In these equations we see the connection between the “matter” 
aspect — on the left-hand side — and the “field” aspect — on the 
right-hand side — both concepts and their connection being well 
known since the formulation of classical electrodynamics. 

We turn towards the more important problem of finding the 
consequences of the equivalence of the (X) and {IQ) forms. We 
can write both these forms in one single equation, which will be 
fundamental to our further reasoning: 

J A^nmdS = - f (A^+Sn&^dx = - fK m0 ’ a/, l/} n m dS. (6.5.11) 

A A A 

£ Q £ 

Prom the last equation follows: 

J (A 0a +S7t7 0a )dx = J ET^^da+Cr (6.5.12) 

A A 

o r 

where the C a are constants of integration. Yet it can be shown that 
because of the field equations, these constants must vanish. Indeed, 
the field equations (6.2.1) and (6.2.3) give, for n = 0: 

4*+8*9* = -(6-5.13) 

A 

By integrating the last equation over Q we obtain: 

J (/l 0 *+87r9’ 0 ”) dx = f K^^^aS. 

A a 


(6.5.14) 
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We remember that JC W0,V ^ contains only linear expressions in 
the y 1 s and (6.2.2) allows us to write them out explicitly. They are: 

. (8.5.15) 


6. THE THREE LINEAR MOMENTA 


Because of the form of 7°°, we can write (6.5.14) in the form: 


~ J A^dx+fa = (6.6.L) 


We shall call the “vector” 

A 


& 


= P 


(IN) 


( 6 . 6 . 2 ) 


the inertial linear momentum of the A’th particle. This linear 
momentum depends on the particle, its velocity, and the fiold at 
the points through which this particle moves. This is so because 
the n's depend on the tweedled field. 

We shall call the “vector” 


A. J A - lx _ (8.0) 

A 

Q 

the field linear momentum around the A’th particle. It will depend 
on the region over which we integrate. 

We shall call the “vector” 
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the gravitational linear momentum of and around the JL’th particle. 
We can simply write (6.6.1): 

P(Q) = +P(F)- (6.6.5) 


A 

But PJg) is defined only through surface integrals of linear expres- 

A 

sions of the first derivatives of the y’s, whereas P*^ is defined 
through a space integral of a non-linear function of the y’s. 

We should like to know the physical meaning of these three 
linear momenta, at least in some special cases. To do so we must 
quote some results obtained before with the help of the approx¬ 
imation method. 

We shall make a special assumption: we are dealing with a two- 
body problem and we are only interested in the 0-component of 

i i i 

the last equation. We denote by /tqui), ft^, P(Q), the inertial, field, 
and gravitational mass of the first particle, corresponding obviously, 

since f° = 1, to P“in) , and P®^. Similarly, for the second 
particle, we shall write “2” above the masses. Now we write for 
any of these masses: 


fi = (t+fj, 
2 4 


( 6 . 6 . 6 ) 


where n is the mass in the Newtonian and n in the post-Newtonian 

2 4 

approximation. For (jl we have simply that in this approximation 
2 

the inertial and the gravitational masses are equal, both being 
constant, whereas the field mass equals zero. Now about the 
calculation of /*. The result is, according to (3.3.29): 

4 


i i . . n 

I + PM r u • 
4 2 2 2 


(6.6.7) 


Here r 12 is the distance between the two bodies. 
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Let us row calculate . A straightforward calculation, not 
worth repeating in detail, gives: 

A™ = -I9V9V, <P = —\ ~r m (6.6.8) 

* 3Ii r 

Therefore, in the case of two particles 

8&M(F) = | JV|«,<pdx—| 

4 1 1 

0 £ 

1 2 

= -1231 - 4 (6.6.9) 

From this equation and (6.6.6) we can find 

1 1 2 . 

/*(©) = /“+£/“£*£*— — —— ——— [<p\ a <pn a d8. (6.6.10) 

2 2 jj IF•y 2 

£ 

We see that the gravitational mass, like the field mass, depends 
on the finite region over which we integrate. Only in the Newtonian 
approximation is there no difference between the inertial and 
gravitational masses. 

7. THE EQUATION EOR GRAVITATIONAL RADIATION 

In Section 4 we obtained the equations of motion from the 
differential conservation law 

=0. (6.7.1) 

We shall use the same equation from which we obtained the laws 
of motion in the {EQ) form as a means to define the gravitational 
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radiation. The only difference between the point of view represented 
in Section 4 and here consists in the simple fact that now we do not 
integrate over the neighbourhood of the A’th singularity but over 
the entire space, which we shall denote by Q. The infinite spherical 
surface surrounding such a space we shall call 2. Inside the surface 2, 
we have particles moving according to the laws of motion deduced 
from those of the field and formulated in the three last sections. 
We assume that the particles never reach the surface 2. Thus, 
integrating (6.7.1) over the entire space, we find: 


j'A^nndS = - J(i"° + 87t'>Vx. (6.7.2) 

2 a 

We shall regard by definition the left-hand side of this equation, 
that is 


f A Mm n m dS (6.7.3) 

2 

as the flux of gravitational radiation. Thus equation (6.7.2), with 
the help of which we have defined the gravitational radiation, is 
closely related to the equations of motion. The only difference 
between (6.7.2) and the equations of motion consists in the region 

A A 

of integration: 2 and £2 instead of and £2. Equation (6.7.2) is the 
{2Q) form of the conservation law, in the same way as previously 
(6.6.3) was the (2£2) form of the equations of motion. 

As before, so now we can write the conservation laws in the {2) 
form: 

/ A» m n m dS - J (6.7.4) 

2 X 

or, confronting (6.7.2) and (6.7.4) we bog that 

f A^tindtS =-/ {A^ + MJ^dx = - / («- 7 6 > 

SO £ 

Motion and relativity IS 
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or, as before in Chapter YI, Section 5, we conclude: 

‘J(vl # "+8jr7 0# ‘)dx = fST^^dS. (6.7.6) 

Q 2 

Finally we may also write the conservation equations in tho (Si) 
form: 


f 


ii / 


* r 

= V - 


iri L* 


A A A 


All these equations are a consequence of the equations of 
motion, which in turn are a consequence of the field equations. 

We know that the conservation laws in the (J?) form, that is 
(6.7.4), are independent of the shape of the surface, as long as this 
surface does not pass through singularities. Therefore, we may 
form the two-dimensional surface, I enclosing all singularities as 
consisting of small spheres, each of them surrounding one of the 
singularities and connected with each other by very thin tubes. But 
because the flux through these thin tubes equals zero, wo can 
write (6.7.4): 

J [A^+E^'%)n m d8 = £ JIA’»+Br*# v )« m da = 0. (6.7.8) 

Z A-1A 

£ 

Similarly we may write, using the (IQ) form: 
x N 

Z l^mdS = - 2 J (>+8^°)dx. (6.7.9) 

AmtlA A-l A 

27 Q 

Yet, generally 

N 

J A m n m aS ^ / A^n^S. 

£ 


£ 


(6.7.10) 
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We regard the right-hand side and not the left-hand side of 
this inequality as the definition of gravitational radiation. This 
definition is not precisely the usual one, because A 1 "” 1 also contains 
derivatives of the second order. But it seems to me more consistent 
to use for the gravitational radiation this definition which is closely 
connected with the equations of motion. 

Similarly as before we define the total linear momentum: 

n 

P° { iN) = J (6-7-11) 

Q A -1 

the total linear field momentum: 

fA°'dx (6.7.12) 

ora a 

and the total linear gravitational momentum: 

Pfo - ~ / R m0 ^„n m d8. (6.7.13) 

Similarly as before, we can write (6.7.6) in the form: 


Ph (6.7.14) 

Here we face a situation that appears at almost every step in 
electrodynamics. Quantities referring to the world-line as P\i jy 
does, are, connected with those referring to the field, like Pfa 
and Gi¬ 
lt is instructive to compare the situation as presented here 
for the gravitational field with that of an electromagnetic field 
in Special Relativity Theory,, in which case we have the electro¬ 
magnetic energy-momentum tensor JE^ and the conservation law 
= 0. Prom it by integration over the entire space we find: 


(6.7.15) 
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The right-hand side is usually understood as the (negative) 
time derivative of the entire field linear momentum. For the 
electromagnetic field it is the only linear momentum that appears 
in the conservation law. In the case in which the loft-hand side 
representing the flux of electromagnetic radiation vanishes, 
we conclude that the field linear momentum is a constant, 
vector. 

Th^ situation is different in the case of a gravitational field. 
The flux of gravitational radiation is equal to- the (negative) time 
derivative of the sum of two linear momenta: that of the gravita¬ 
tional field and that of inertial matter. In the case in which the 
gravitational radiation flux vanishes, it is the sum of the two 
momenta that remains constant. This sum we call the gravitational 
linear momentum and it is equal to a surface integral of an ex¬ 
pression consisting of linear derivatives of the y’s. 

. To give an example showing the difference between these three 
total linear momenta let us return once more to the example given 
at the end of the previous section; that of two particles and their 
total masses calculated in the Newtonian and post-Newtonian order. 
As before we have in the Newtonian approximation: 

P(iN) = P”®) = const; P° (F) = 0. (6.7.16) 

2 [2 2 

Then, in the post-Newtonian approximation, since the surface 
integral (6.7.10) over the infinitely great surface vanishes, 

.P(iN) = Total Kinetic Energy — 2 Potential Energy. 

-Pfo) = Total Kinetic Energy + Potential Energy =.= const. 

Ffa = 3 Potential Energy. 

We see from the expression for that in this approximation 

the surface integral oyer “the normal component of the Poynting 
vector” vanishes in this case for gravitational radiation. This is so. 
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because 

Pfo+Pfo = —~ J A m n m dS = 0, (6.7.17) 

that is, because of the Newtonian equations of motion the sum of 
kinetic and potential energy is constant. 

8. ON THE INVARIANCE PROPERTIES OF 

We ask: what are the transformation properties of Pfo under 
small transformations of the type 

= af + «*»? (6.8.1) 

We shall consistently neglect all* products of the «’s as well the 
products of the y’s and o’s. We recall that, according to (6.7.13), 
(6.5.15) and (6.5.16): 

Pfo - ~ f &**»*.** = ~ J(-y°V+ /" isMS, (6.8.2a) 

2 2 

r i’ 

(6.8.2b) 

We start by investigating the changes caused in P®g) by trans¬ 
formation (6.8.1). Thus wo have to use the formula: 

fix 

det * (6- 8 -3) 

Into this equation we have to introduce: 

( 6 - 8 - 4fi ) 

dx 

det —v = 1- a°, 
d» 1 


(6.8.4b) 
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From the last three equations it follows (since gf* v = rj^+y**) 
that 


y* 00 = y 00 +a 0 ,o-«V, 

y *om = y m _ a O |m+aTO|oj 

y*™ = y mn ~a m , n -a n lm + S mn a%+ d mn a° l0 . 

Introducing these into (6.5.16) and (6.5.16) we find: 


(6.8.5a) 

(6.8.5b) 

(6.8.5c) 




\fi 


*•"%--r-V- (*V-«"V,., 

_<r* «•„),.+ 


|m’ 


(6.8.6a) 

(6.8.6b) 


The bracketed expressions on the right-hand side of (6.8.6) 
are skew-symmetric in the indices m, s; therefore, because of the 
lemma proved in Section 2, we have: 


_ Jurats, ( 6 . 8 . 7 ) 

s 2 : 

or 

which shows that the components of the gravitational linear mo¬ 
mentum are invariant under small transformations of the coordinate 
system. This-would not be so if the products of the y’s and a 's were 
taken into account. We can see this immediately in the special case 
in which the gravitational flux vanishes. There, with respect to 
a Lorentz transformation the A 0a and T°° behave like tensor com¬ 
ponents, and therefore Pfo like a constant vector. Thus with respect 
to a Lorentz transformation the Pfo are not invariants. 


9. GRAVITATIONAL RADIATION AND THE CHOICE OF 
A COORDINATE SYSTEM 


We arrive now at the principal question: can we find a reasonable 
coordinate system for which the gravitational radiation vanishes? 
in other words: is the gravitational radiation as defined here some¬ 
thing that can be annihilated by a proper choice of a coordinate 
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system, or is it something that has an absolute meaning? This 

means: can -we find a coordinate system for which Pfo = 0 and 
P“(Gt) = constant, that is for which the law of conservation of the 
sum of the two linear momenta — the inertial and the field — is 
valid? 

Before we answer this question let us recall once more the 
definition of It was given in (6.7.13): 

Pb) = 8 ~ / (6.9.1) 

But because of (6.7.8) we can also write: 

P\g) - - J A^n^dS. (6.9.2) 

We can, therefore, calculate the change in linear momenta in two 
ways: first, through the linear expressions in the y’s, that is through 
(6.9.1), second through the non-linear expressions in the y’s, that 
is through (6.9.2). 

Now about the choice of our coordinate system: all our 
assumptions will concern only tho behaviour of our coordinate 
system and the y’s for r-»-oo, r being the “distance”'from a fixed 
point, that is r® = x k x k for very large r’s. We assume that for 
r -*■ oo our metric becomes Galilean at least like logr/r", with a > 0, 
that is: 

l/*l< . (6.9.3) 

T 

Furthermore we assume for r ->• oo tho following coordinate 
conditions: 


/ 1 \ 


= 0 {y+Zf’ > 

(6.9.4a) 

r'v-oUA 

(6.9.4b 
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We should like to draw attention to the condition (6.9.4b) which 
is different from de Donder’s condition y"*° !a = 0 used so extensively 
by Tock. The results presented here are essentially due to the choice 
of the coordinate condition (6.9.4b) which takes for granted the 
ideas on which the “new approximation method” is being built: 
that space and time have to be treated differently; this difference 
manifests itself in the choice of our coordinate conditions. Indeed, 
it was this coordinate condition that was chosen in earlier papers 
by Einstein, Infeld and Hoffmann. The chief difference between 
this work and that of other authors, notably Trautman, also 
consists in the different coordinate condition used here. 

How for this special coordinate condition, the take a very 

simple form: » 

= -iy°V> (6.9.5a) 

*’'“%= ~K/V+/%), (6.9.5b) 

as follows from (6.8.2). Also the field equations take on a simpler 
form in this coordinate system. At this moment wo are only 
interested in their (0, 0) and (0, Tc) components, and only for r -> oo, 
where we can forget the energy-momentum tensor: 

-£°° ,< V+A 0# = 0, (6.9.6a) 

* 0 %,+A°* = 0, (6.9.6b) 

or because of (6.9.5) we have in our coordinate system for r -* oo) 
simply: 

*y°V = -A #0 , (6.9.7a) 

= —A 0 *. (6.9.7b, 

Let us concentrate for a while on the second of these equations. 
As we know, A is a product of the y’s and their derivatives. If the 
y’s are of the order 1/r, the derivatives are of the order 1/r 2 or higher 
in 1/r. Therefore A 0 * must be at least of second order in 1/r It 
would seem, therefore, that A m can make a contributionto. the 
surface integral (6.9.2), and therefore that JPfo can be a function 
of time m this coordinate system. However this is not so. 
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We develop the essential expressions in A m , which may give 
a contribution to the surface integral (6.9.2), into a power series 
in x a jr, assuming, of course, that this is possible. That is: 


A m 






+• • • 


(6.9.8) 


where the coefficients a m etc. are functions of time only. Since, 
however, the left-hand side of (6.9.7b) is a pure Laplaeian it follows 
that the corresponding y m must be of the form: 

/yj|® 

y m = A 0m logr+A 0m 8 — +... (6.9.9) 

Here the A m etc. are such functions of the time that the equations 
(6.9.7b) are satisfied. But that is in contradiction to our assumption 
that the y m , like all the other y’s, become Galilean at infinity. 
Therefore A m can have no expressions of the order r ~ 2 ; it must 

start with a higher order in 1/r; therefore P°g) must equal zero, 
therefore P“©) must be constant! 

This result is bound very definitely to the choice of our coor¬ 
dinate system. It forces the appearance of a Laplaeian, instead of 
a d’Alembertian on the left-hand side of (6.9.7b). Indeed in the 
case of a d’Alembertian we could have had expressions for y om o. 
the order a[r and a/r a which are consistent with A m of the order r"*f 
Turning to the definition of Pfoas the integral of the linear 
terms of the y’s we see that 

P\o) = -g~/ ~ T7- J /V o^d8 - 0. (6.9.16) 

n £ 71 £ 

Let us assume that 

4Jf 1 

y 0# =-h terms of the higher order in • . 

T T 

Then, by taking the surface integral, we find: 

1 rim 


(6.9.11) 
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Therefore M must be constant and is equal to the total gravitational 
mass. 

We shall now turn towards the We shall prove that we 
may further specialize our coordinate system at infinity so that 
they, too, will be constant. That means: the coordinate transform¬ 
ation that we shall introduce now will not change the essential 
features of the previous coordinate system. By this we mean that 
the y*’s will vanish at infinity like log rfr a [a >0) and that y* 0 ^ 
and y will go to zero at least like l/r 2+ ° with a > 0. 

For the sake of simplicity we shall make a further assumption: 
that for r -»■ oo, y 0a goes to zero like 1 /r. We could have stuck to our 
assumption, that y 0a behaves like log rfr a with a > 0 for r ->■ oo, 
but our argument would then be slightly more complicated. 

We introduce the coordinate transformation: 

x* k = -x k + (6.9.13) 

where the C k are functions of t only. We shall, however, in trans¬ 
forming the y’s write out all expressions, linear, or non-linear, up 
to the order of l/r 2+ci with a > 0. From (6.8.3) and (6.9.13) wo find: 


y* 0 "* = y m +CTy +<T^ +0(^), 


y*™ 1 = y™-(rl- 

r) 

-) +d mn C‘l- 

r) + 

\ 1 

r / |n ' 1 

? 1 \m \ ' l 

r '1* 


• v®** • v m • • 1 / 1 \ 

+ crL- + cr^—+G m c n ~-+o[-^). 

r r r a \r* + / 

Since in (6.9.14b) y 00 | m /r is, at least, of the order 1 jr t+a 
we see that for r -> oo: 


(6.9.14a) 

(6.9.14b) 


(6.9.14c) 
(with a > 0) 


y io+/ 



a > 0, 


(6.9.15) 



VI. 9 GRAVITATIONAL RADIATION AND COORDINATE SVSTEMJ 187 

and similarly 



since y 0m is of the order 1/r and its space derivatives of the order 
XJt . We know that the linear terms in (6.9.14) do not have any 
influence upon P\ 6) . This was shown explicitly in the previous 
section. Furthermore, no terms of the order l[r 2+a have any influence 
upon the surface integral. Finally since in (6.9.5b) y* ok is differen¬ 
tiated with respect to space, it is sufficient to take into account 
in y*° only the expressions of the order 1/r. We see, too, that in 

(6.9.14c) the expression CTCT/r 8 cannot make any contribution to 
the surface integral, because it is multiplied by n m = x m fr and the 
surface integral vanishes because of the non-symmetry. Thus for 
the sake of finding P*q } , we may assume: 

y* 00 = y 00 , (6.9.17a) 

y ,#m = y m , (6.9.17b) 

y =7^+0“^ +O n ^ r . (6.9.17c) 

We shall introduce now the following functions of time: 

1 r y 0n 

^ s J c«- 9 - i8 > 

I 

Therefore, because of (6.9.1) and (6.9.5b) wo have: 

2(G) = ^(G) - (A£O™-fri*C fc ) |0O . (6.9.19) 

Thus, generally, we can, without invalidating our original conditions, 
annihilate -£*(<?). To do this we must solve the differential equation: 

JPfo = (<V m +A° 3 C k ) m . 


(6.9.20) 
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III it Pjo) and are known functions of time. Prom them we find 
the O’s and introduce with their help the coordinate transformation 
(6.9.13). In the new “starred” coordinate system we have: 

Pfa = 0, P*S) = const. (6.9.21) 

Thus in such a coordinate system the sum of the inertial and field 
linear momenta remains constant! 

|ln the case in which for r -> oo we have y <m — A m (t) , 

we would have to replace (6.9.13) by 



As far as we can see, there is only one argument that could be 
levelled against this reasoning: what happens if vanishes? 

This question can easily be answered by a generalization of the 
argument just presented. 

Indeed let us introduce instead of the transformation 

* C l 

x* k = a ?+■>— (6.9.22) 

r 

a slightly more general transformation 

x* k =a/ c +a k (6.9.23) 

where a k has the form: 

* b k bW 

a =T + ~^~ + - (6 ' 9 ' 24) 

and the 6’s in (6.9.24) are functions of t only. In other words we 
assume: a k is of the order r~ l and does not change its order by a tirriA 
derivative but lowers it by a space derivative. Thus we have: 



(6.9.25a) 
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y*om = + fl * o+ a m ]0 y°°+ o , (6.9.25b) 

y*' nn = y mn -<^ [n -a n {m + <5 mn a s |8 +a m |0 /" + 

+«’ l |oy 0m + a n |o« m |o+o(-^). (6.9.25c) 

oaly products o£ the type y • a and a • a can make contributions 
to P(G) and then only if they are of the proper order, we may assume 
(in order to find Pf^), similarly as in (6.9.17): 

/ #0 = y 00 , (6.9.26a) 

y* m = y™ 1 , (6.9.26b) 

=y™+«\Y M +a\y m +a m {0 a n l9 . (6.9.26c) 

Therefore, because of (6.8.2b): 

P (0) = ^ f («V*“+ a m {0 y ok +a\ 0 a m lo ) [0 n m dS: (6.9.27) 

£ 

Thus we may annihilate P*^ by any transformation satisfying 

*&) = J hf (a * ay0m+ fl * °y° k+ a > m ,o),o n m dS. (6.9.28) 

£ 

Even if y m is of the order r~ 2 , such a’s will always exist. In such 
a special case we have: 

%*(&) = J ( <,OT |oa fc |o)|o w, m^^* (6.9.29) 

71 £ 

It is sufficient to choose the a’s of the form: 



then we get the equation for a k and f!i 

1%) = i(a fc /?),o. 


(6.9.30) 
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This transformation leaves the coordinate conditions essentially 
invariant as we see from (6.9.25); that is we have: 



(6.9.31a) 



(6.9.31b) 


[10. ON THE GENERALIZATION OP THE COORDINATE SYSTEM 
Until now we have assumed that for the metric field at infinity: 

1- |y°*| <|jf*(*)-^- , tr > 0; (6.10.1) 

2 - f“|. = o(^). /"V-o(^); (6.10.2) 

3. A om , as far as it gives a contribution to the surface integral, 
can be developed for r oo into a power series in x a jr, the coeffi¬ 
cients being functions of time. 

We proved that under these conditions: 

A. P(Q) is constant; as a matter of fact to prove this we did not 
need to assume 1. It was sufficient to assume that 

\y m \ < M m (t) ^ with a > 0. 

B. The P“g) are invariant under a transformation 

■ x* a = x a +a a (6.10.3) 

if for r -> oo we neglect all products of ya and a-a. 

We shall csdl 

x* a =x a +a a (6.10.3) 

a proper transformation if, for r oo the following conditions are 
satisfied: <*“ is of the order 1/r; the space derivative lowers the 
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ON THE GENERALIZATION OF THE COORDINATE SYSTEM 

order by one, whereas the time derivative leaves the order unchanged; 
or: a a can be developed into a power series in x*fr with coefficients* 
depending on time only. Furthermore we assume that for a proper 
* ra 2 +ff 0r ^ aa ^ 0u K* m0,afi u) are to be calculated up to the order 
1 l r + a > 0; that is if necessary the non-linear terms in the 

transformation of y ttfi should be taken into acco unt 

With such a definition of the proper transformation we have 
proved: 

0. P®#) and the coordinate conditions (6.10.2) are invariant 
under a proper space transformation, that is if: 


x*° = x\ 

(6.10.4a) 

x* k — x k +a k . 

(6.10.4b) 

Also we proved: 

D. We can always choose such o’s in the proper transformation 
(6.10.4) that 

o’ 

II 

(6.10.5) 

Now we ask: what happens to and the coordinate conditions 

(6.10,2) if we nse not the proper transformation (6.10.4), but the full 
proper transformation (6.10.3)? To answer this question it is sufficient 
to consider only the proper time transformation, that is: 

Wo have: 

(6.10.6) 

y *oo = / o +a o |0 + a° |# /o+ O |-~-J, a > 0, 

(6.10.7) 


(6.10.8) 


(6.10.9) 


But the lineal' terms do not make any contribution to P^. The 
non-linear terms are of the order l/r 1+a and if differentiated once 
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more with, respect to the space coordinates will give no contribution 
to the surface integral. Therefore: 

E. Pgjj is invariant with respect to a proper time - space trans¬ 
formation. However, we can easily see that such a transformation 
does not leave the coordinate conditions invariant. Indeed we have: 

Y*\ = a0 |00‘^ a °i«0 a °l0+ ( ffl0 |o/°)|oH- 0 I'^S+o j> (6.10.10a) 

/“V = (-^). (6.10.10b) 

Here the a°’s are arbitrary, as long as the transformation is 
a proper one. From the above equations we conclude 

/“Vo ~ / V = O . (6.10.11) 

This means: we can prove that P®^ is constant under somewhat 
less stringent conditions than those enumerated at the beginning 
of this section. We can assume only the three coordinate conditions 
(6.10.11) instead of the four (6.10.2). The coordinate condition 
(6.10.11) if explicitly stated at the beginning of our argument would 
seem artificial. Here, however, we were led to them in a natural 
way by investigating the transformation properties of P®^. How¬ 
ever, one can give a simpler and more direct proof of the following 
theorem: 

F. We assume for r -> oo: 

1. \y m \ < ; (6.10.12) 

| 

2- rV-rV = o(^). «>0; (6.10.13) 

3. A m can be developed into a power series in x*/r, as far as 
the expressions significant for the surface integral are concerned. 
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Then, if these conditions are fulfilled, we have: 

P( 0 ) = const. 

The proof follows along similar lines to that given in Section 9 
and we shall indicate it here only briefly. We have, according to 
(6.7.13) and (6.7.5): 


= h --if 



(6.10.14) 


According to (6.9.6b) and (6.5.16) we can write the (0, fc) field*equa¬ 
tions in the form: 


i(y°V—y°“iaft+ y 8fc i« 0 ) + A ° k = o. (6.10.15) 

But, because of coordinate condition (6.10.13) we have, identically 
with (6.9.7b): 

= ~ A ° k + 0 (^r) • (6.10.16) 

From now on the argument is identical with that based on 

equation (6.9.7b) and leads to the conclusion that JP®^ = 0, 
therefore Pfa must be constant. We can see this in a still simpler, 
almost trivial way. Since y°* as a solution of (6.10.16) can be 
developed into a power series in x*\r, we conclude that 

“>°> (6.10.17) 

therefore our condition (6.10.13) can be written: 

y m Vo-/V = o(^r). (6.10.18) 


Motion and relativity 


IS 
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But a comparison with (6.8.2a) shows: 

y m ”,no-/°|om = -2^°’°^; (6.10.19) 

therefore obviously P®^ must vanish. 

The coordinate conditions (6.10.13) are, like P^, invariant 
with respect to any proper space-time transformation. From equa¬ 
tion (6.10.14) we also see that P®^ can be presented as a space 
integral. Obviously it does not depend on the coordinate system 
inside the region bounded, by the infinite sphere, as long as it leads 
to a proper transformation on the sphere. It will also be invariant 
with respect to a pure space transformation which for r oo has 
the form x* k = a?+a k (x), where a k is of the order zero in r, as long 
as its derivatives lower its order in r by one. Therefore it would 
seem that P®^, defined by (6.8.2a) and equal to the total gravita¬ 
tional mass, deserves to be regarded as the energy of the system. 

. It can be shown, however, that P®^ will not be constant in -all 
coordinate systems. Indeed, let us consider a simple transformation 
(as an example) for f -» oo: 

x'° = a?®-f slog-^-. (6.10.20) 

Jz 

Here e is constant and R is the very great constant radius over 
which we integrate. A straightforward application of the formula 
(6.8.3) gives: 

V" = >."+-£(6.10.21» 

8 

7 >m = y° m +e—+8y nu - J , (6.10.21b) 

y™ = y mn . (6.10.21c) 

How to apply them to our coordinate condition, we must remember 
that 

d d a d 
dx' p dx fi ea daf ' 


( 6 . 10 . 22 ) 
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Therefore 

y'V^V-VV^+O^) (6.10.23a) 

/”*,» - +o(^r)- (6.10.23b) 

The additional expressions are of the order 1 [r 1+a and they 
may change the surface integral for a = 1. Therefore the “gravita¬ 
tional radiation”, or rather the flow of its energy, can he created 
or annihilated by a choice of a coordinate system. However, as we 
have shown, there exist reasonable coordinate systems in which 
the “gravitational radiation” always vanishes. 

11. RADIATION AND THE APPROXIMATION METHOD 

Up to now we have shown that by a choice of the coordinate 
system we can annihilate the radiation. But we can also .show that 
it be created by a proper choice of the coordinate system and 
by a proper choice of the solution of the field equations. 

The method used in this chapter, in all sections with the excep¬ 
tion of the first, had little to do with the approximation procedure. 
We wish to come back to the ideas of the first section, applying 
them to General Relativity Theory. In the first section, the radia¬ 
tion depended on the equations of motion. Here we shall show 
how it depends on the particular solution we choose. 

Again we shall fix our coordinate system by adding four 
conditions, which need to be fulfilled only at infinity. These condi¬ 
tions will obviously not be of the type 

^10 = 0, /"",* = 0, (6.11.1) 

since they lead to no radiation. Instead we shall use de Donder’s 
condition 




(6.11.2 
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and investigate its influence upon radiation. How does this coor¬ 
dinate condition change the field equations? If we apply it to 
(6.2.3) and (6.2.2) we obtain very simple field equations: 

-iri af y M, iafl +A' tt, + 8n7 ,a = 0. (6.11.3) 


In the linear case it is the condition which leads immediately to 
the wave equations: 

□ y'"' = 0. (6.11.4) 

This is also the reason tor the popularity of this coordinate condition. 

How let us try to solve the field equations at infinity with the 
help of our approximation method. That is, we start with the 
solution for big r’s of 

Ay<>o = 0, (6.11.5) 

2 


and according to what was said in Chapter II: 



2 


N A 



( 6 . 11 . 6 ) 


For r-*- oo we can develop this expression, similarly as we did in 
Section 1: 


y 00 =—) +2B 9p l-) +... (6.11.7) 

s r \rl i« \r/i«p 

Here 

X s = M-'Yii*. (6.11.8) 

A 2 A 

n sp =££(1*t*-id tp £ r f r ). (6.n.9) 


The genesis of the last expression is simple; it springs from the 
following equality 
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radiation and the approximation method 

where p is an arbitrary function of time. Therefore it can always 
be so chosen that 


= 0 - ( 6 . 11 . 11 ) 

Now we turn towards the field equation for y 0 ” 1 . We have: 

. 3 

y m ise = 0 . ( 6 . 11 . 12 ) 

3 

Because of our coordinate conditions we then have: 

r- 4jfi ’r -^(vl+^L+dv),- (6 - aj3> 

Here p is an arbitrary function of the time and / sm an arbitrary 
skew-symmetric function of t. If we wish y m to vanish at infinity, 
we must assume M = const, because of the coordinate condition 
v\ = 0 . 

Now we can calculate / m \ Again, because of the coordinate 
condition, we have, if y mn is to become Euclidean for r -> oo: 

4 

X* = 0, = 0, (6.11.14) 

and 

y™ = (2i) mB + pd™ 1 ) -+... (6.11.15) 

4 4 r 

Now if we wish to take the radiation into account we must 
properly choose 

y 00 , y°", y mn , (6.11.16) 

$ [4 5 

As was shown in Section 1 of this chapter, jumping by two steps 
would give only “standing waves”, and therefore no solution 
corresponding to the retarded potential which alone gives hope 
for incorporating the gravitational radiation into our approximation 
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method. Therefore we should have to start with 


y 00 = -4Jf, y m = ±MX m , y™ = -(2(6.11.17) 

8 4 5 

Of these, only the “radiation term” in y m “ is different from zero. 

s 

Let ns, therefore, assume its existence and see .whether it leads to 
a finite energy radiation. Thns assoming, according to (6.11.4), 

X k = Qj jwn _ 0> (6.11.18) 

we have: 


r ._£ + w ( l ) w+ ... 1 


y™ = (21)^+a 77 ” 1 /?) - 2D mn — 

r 


To find the energy radiation we shall now use formula (6.7.13): 


= f K n0 ^^d8 = - f A om n m dS (6.11.20) 
2 2 

which in our coordinate system, because of (6.6.15), gives: 

-8 = / (y°V + / m |oo)«* = jA^n^dS. (6.11.21) 

2 2 

We may expect the radiation terms to be different from zero for 

the even orders of A m . The first • order in which this can happen 

is for A om . Indeed, without going into detail, we see that such 
8 

a contribution may come from the expression in A m which is: 

( 6 . 11 . 22 ) 

5 12 
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Such an expression, because D as = 0, will be proportional to 


pM. (6.11.23) 

6 2 

Therefore, it is generally possible to obtain radiation of the eighth 
order, but only if for r oo we have solutions of the kind (6.11.19) 

with p =£0. Otherwise, if we take a different world, with p = 0 
6 • 

we do not have any radiation of the eighth order. Thus the radiation 
can be annihilated or created by the proper choice of solutions. 
The simplest solution with p = 0 leads to no radiation. If we take 
P =£ 0 then we have radiation! From the left-hand side of equation 
(6.11.21), that is from the linear part of the surface integrals, we see 
that they must be equal to the right-hand side already calculated. 
This means that we must assume 

, (6.11.24) 

6 \r/|»n 


so that the integral of the left-hand side will be exactly equal to 
the contribution (6.11.23). If /? = 0, then of course <7 = 0 too. 

. Now let us assume that we have chosen p =-0 in (6.11.19). 
Then to tind the radiation we must proceed one step further, 
calculating A 0 ” 1 of the tenth order. For example, there will be ex¬ 
pressions of the type 

/"W (6.11.25) 

7 3 

Now 

(6.11.26) 


That is: y mn will have an expression of this type generated in the 


next approximation by D n . On the other hand, y has an expression 

S 3 

of the type 
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Therefore the product (6.11.25) will he of the order 1/r 2 and will 
give a contribution to the surface integral. It will be of the right 
order and of the form: 

i) mn ^jJ) mn . (6.11.27) 

at 

Generally we shall have many expressions of a similar form; this is 

A m ~f(lP). (6.11.28) 

10 10 

It would seem, therefore, that a radiation must appear in the tenth 
order. However, this is not so. At each approximation step we may 
add an arbitrary function of time. Therefore, let us say that in 
the 7th order we may add to y mn an arbitrary function of time of 
the form 

p(t) d mn . (6.11.29) 

7 

As we have shown before, such an added expression will give 
a contribution of the form 

pM. (6.11.30) 

8 1 

Therefore, again, by a proper choice of an additive harmonic function 
of the seventh order we can annihilate the radiation created by an 
additive harmonic function of the fifth order! We may add that 
such a harmonic function does not upset the accepted coordinate 
conditions. 

What is the moral of this chapter? The results are indeed meagre 
and mostly of a negative character. They show that it is hardly 
possible to connect .any physical meaning with the flux of energy 
and momentum tensor defined with the help of the pseudo-energy- 
momentum tensor. Indeed, the radiation can be annihilated by 
a proper choice of the coordinate system. On the other hand, if 
we use a coordinate system in which the flux of energy may exist, 
then it can be made whatever we like by the addition of proper 
harmonic functions starting with y mn as a function of time. 
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In the linear theory we were faced with the choice between 
the retarded and advanced potential. Here in the theory of grav¬ 
itation the choice is not so simple. Using the approximation 
procedure, we are faced with the choice between single and double 
jumps. We can speak only about radiation in the case of single 
jumps. However, its existence or non-existence or its value will 
depend upon the choice of arbitrary harmonic functions. 

Yet the full equation (6.5.1) makes sense, though not if we 
treat each of the sides as the definition of the energy-momentum 
flux. It makes sense since it is a consequence of the equations of 
motion. Indeed, the surface integrals taken around the singularities 
give us the equations of motion. Taken over the surface surrounding 
all singularities, they give us conservation laws which are a conse¬ 
quence of the equations of motion. 

Einstein often remarked: “We do not have any satisfactory 
classical theory of radiation. Bitz understood this fact. He was 
an intelligent man ...”i This remark seems especially apt if applied 
to gravitational radiation. 
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1. THE 8 FUNCTION 

la this work we use mostly Dirac’s <3 fuactions of a special kind; 
therefore we shall discuss their structure in detail. 

The ordinary d functions may he introduced into the theory 
in various ways. We shall mention three of them: 

(A) : the axiomatic method, 

(B) : the method of Fourier transforms, 

(C) : the realistic method, treating the 8 functions as a limit¬ 

ing case of spread- out sources. 

We shall begin with the one-dimensional 8 function. The 
generalization to a d funotion of more dimensions is almost trivial. 

(A) : In the axiomatic method we stipulate for <3(a>) the following 
properties: 

(AJ: one can formally differentiate d{x) as many times as one 
wishes, 

(A a ): for x =£ 0 we have d(x) — 0; for x = 0, we have d(0) = oo, 

(A„): for an arbitrary neighbourhood V(x 0 ) of x 0 and for an 
arbitrary function f(x) continuous at x 0 we have: 

OO 

f dxd(x—x 0 )f(^x) = f dxd(x—x 0 )f(x) =f(x 0 ). (A.1.1) 

V(io) -oo 

(B) : In the method of Fourier transforms, the symbol 8{x'—x") 
is defined by the following equation: 

8(x'-x") = £ 
x 


(A.1.2) 
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Hero the y> x (as) are the eigenfunctions of an Hermitian operator A 
labelled by the eigenvalues A. 

We assume that differentiation and integration are commu¬ 
tative with summation over the spectrum. Prom this and from the 
properties of a complete set of eigenfunctions one “proves” the 
validity of (A x —A 8 ). Of special importance is the definition of 6{x) 
through the eigenfunctions of the operator id Idee, that is {2n)~ 112 e~ ikx 
which according to (A.1.2) gives: 

00 

8(x'-x") = — f dhe ik(x '- x "K (A.1.3) 

2tv J 

—oo 

Such an “ordinary” Fourier representation of the 3 function plays 
an especially important role in calculations referring to Green’s 
functions. 

(0): We shall discuss the third, realistic method. In it we 
consider d(*) as a limit as s -*■ 0, of certain ordinary functions 
<5(e, x), with the following properties: 

(0 X ): <S(e,as) have all derivatives as functions of x and e, 

(0 a ): lim 8[e, x) = 0 for x ^ 0, 

*—►0 

(C 3 ): For an arbitrary neighbourhood 7(« 0 ) of a? 0 and for an 
arbitrary continuous function /(as) in F(as 0 ), there exists: 

lim J dx8(e, x—x 0 )f(x) = f{x 0 ). (A.1.4) 

*-*° V(ao) 

Indeed, by going to the limit e -»■ 0 and interchanging the limit¬ 
ing process with the integration in (C 3 ) we Justify the axiomatic 
properties (Aj.-A 8 ). 

(Sometimes we have to strengthen our assumptions by stipulat¬ 
ing how strongly d(e, x) goes to zero with increasing x.) 

The realistic method boils down to the following procedure: 
all calculations have to be performed, not on d(x) but on <5(e, x). 
The transition to the limit s -»• 0 has to be made in the last result. 
Only if such a result does not depend on the special foTm of <3(e, as), 
but is determined by its satisfying (0 X -0 8 ), can w r e expect the same 
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result -when operating, not -with <3(8, x), but with a d{x) satisfying 
(Aj^Ag). In other words, the realistic method can be regarded — 
at least by the physicists — as a justification for the axiomatic 
method. This justification lies essentially in the fact that the class 
of functions d(e, x) satisfying (0 l -0 8 ) is not empty. 

One of the possible examples of functions belonging to such 
a class (perhaps the simplest one) is: 

<5 («, x) = —L- 6~ x2 ^. (A.1.5) 

v2tib 

None of these three methods can be regarded as satisfactory 
from the mathematical point of view. However, generally speaking 
they are sufficient for the application of the <3 functions in theoret¬ 
ical physios. Even this statement is not always true. They are, of 
course, sufficient if in the special calculations only products of 
<3 functions with continuous functions appear. However, in general, 
and in particular here, in the problem of motion in G-. B. T., we 
shall have to deal with products of <5’s with functions singular at 
precisely the point at which the <3 functions are infinite. Hence 
we are faced with the problem of interpreting the expression: 

/ P > °* 

Such an integral is usually considered as divergent. It is this fact 
which produces many difficulties. 

What are the logical foundations of the statement that the last 
expression is divergent for a certain p > 0? Usually one appeals 
to (A s ). From this it appears to follow that the use of <3 substitutes 0 
for the argument x in l/|aj| p . Hence, we conclude that the integral 
(A.1.6) must be infinite. 

This result will seem less well-founded if we consider that A s 
is valid only for functions continuous at least at x = x 0 . Thus the 
last expression is neither divergent nor convergent. It is simply 
meaningless. 
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Indeed we can show that -we can consistently introduce modified 

Dirac 8 functions, -which we shall call <3 and <5, such that for any 
chosen p > 0 we hare either 


00 


/ 

—00 


dx 


h*) 

i*r 


= 0 , 


(A.1.7) 


or 

A 

00 A 

J = * s " 

—00 

ft 

where <o p is any arbitrarily assigned value. Obviously 8 as well as 8 
will depend on the particular p chosen. 

Thus we shall deline 8 by the following set of axioms: 

* * 

(A,) 8 has all derivatives for a? ^ 0. 

ft ft ft ft 

(A a ) 8(a>) = 0 if * # 0. 8(0) is.undefined. 8(x) under the integral 
should be treated as an ordinary function. 

(A,) For every continuous function f(a>) and for an arbitrary 
neighbourhood V(x Q ) of x 0 we have 


i* * 

J dx8(x — x 0 )f(x) =/(*o). 
F(®0) 

(A*) For a certain p we have 


(A.1.8) 


J dx 

m 


m 

i»i p 


Wj, 


(A.1.9) 


where co p is a previously assigned value. ^ s 

Let us start with the realistic introduction of the 5 and <3 which, 
as has been said before, can be regarded as a justification for their 
axiomatic introduction. TVe start from a realistic representation 
of the ordinary Dirac functions 8(e,x) — that is, such that satisfy 
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A 

(Cj —C*). Now we define d[s,x) and 8(e,x) in the following simple 
manner: 

(l 

d{e,x) = a(e) \x\ p x)), (A.I..L0) 


8(s, x) = d(e, x)[l + o) p \x\ p ). (A.t.1.1) 

First, let us show that 8(s,x) satisfies the axioms (An-A*) in 
the limit e -*■ 0. This is trivial for A x and A a . In order to satisfy (A s ) 
we have to choose a in (A.1.10) so that 

oo 

j d(e,x)dx = 1 (A.1.L2) 


is satisfied. This is certainly always possible and a turns out to bo 
~s~ p . (For example for <S(e, x) = (2jrs 2 ) -1/2 e _a!2/s *“ we find 


a 


o-» 


^U+wr(^±i)| 


-X 

.) We can also easily see that 


lim f 8(e,x)jxl n dx = 0 for n> 0. (A. 1.13) 

fi-M) —oo 

From the last two equations and from the fact that S(e,x) vanishes 

quickly enough for x ^ 0, A a follows, precisely as for ordinary 
Dirac functions. 

A 

Now, A t follows from: 


f he, a?) 

I i«r 


dx = a(e) d(e, x)x 


00 

-00 


= 0. 


(A.1.14) 


Therefore the axioms (J.!—A 4 ) are justified. Similarly we can find 
for <3 (a?) that 



(A.X.15) 
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However, we can. now go further. We can disregard any realistic 
representation, and thus write formally the following symbolic 
equations: 


da 

k«n) = a\x\ v — [xb{x% 
d(x) = 0(a>)(l + a> P |af) 


(A.1.16) 


where the S(x) are ordinary Dirac functions and the a’s infinite 
constants so chosen that 

J l dec = 1 = J 8dx. 

F(0) F(0) 

The generalization to more dimensions is simple enough. Let 
us assume a two-dimensional manifold (®, y) and a function f[x, y)j 
so that /(*, jO ->• 0 as x, y 0. We wish to find such 3(®, y), that 

f 5(a?, y)dxd/y = 1 (A.1.17) 

F(0) 

and 

f S(a 5 , y)\f\~ l dxdy = 0. (A.1.18) 

F(0) 

The answer is symbolically written: 

3(a?,y) = (A.1.19) 

assuming, of course, that |/| is such, that the integral (A.1.17) does 
not become identically zero for an arbitrary a. 

From the point of view of our applications the most important 
case is the three-dimensional one in which the <3’s, the fs and 

therefore also the d and 3 depend on r only, where r 2 = cc -\ry . 
Such central symmetrical cases can be reduced to one with only one 

dimension. 
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For an ordinary three-dimensional Dirac <3(x) function we hare: 

j'dxd(x) = 1. 

We can satisfy the last equation, as well as other conditions which 
should he satisfied by the 8(x) function, putting 


d(x) 



where 


/ 


d[r)dr = 


(A.1.20) 


However, for 8(r) and 8{r), besides the normalization conditions 
identical with those above, we also wish the following equations to 
be satisfied: 

A 

f = 0 ail< * f = (A.1.21) 

o o 

Thus if we define: 


«(*)■ 


2ji r* 


d(r), <5(x)-> 


2jw‘ 


Hr) 


(A.1.22) 


we see that 


/ 


dx 


8(x) 


= 0 , 


/ 


dx 


«(*) 


W — M v 


(A.1.23) 


These 3 (or 3) functions behave as kernels of integral operations 
with respect to continuous functions exactly like ordinary Dirac 
<5 functions. They allow us, however, to connect definite meanings 

with integrals of products of these <5 or <5 functions with those that 
become infinite for x -»• 0. From our point of view, this is especially 
important for the <5 functions. We shall assume that they are so 
chosen that the integrals over these 5 functions times functions 
that become infinite at x = 0, vanish in each particular case. We 
shall call such properly chosen 8 functions the “good” 8 functions. 
The application of these functions is equivalent to the regularization 
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procedure, which, in this way is incorporated into the mathematical 
tool which we arc using. 

Sometimes we look for good three-dimensional <S functions 
which satisfy our condition (A.3,.23), not only for a certain p hut 
for a finite sequence of the p’s, say p =1,2,3,..., Te. We shall 
show that we can just as well construct a realistic model of such 
a function. Lot us start with a model of an ordinary Dirac ^function 
satisfying the following condition: 

<5(e, x) = 8~*A |~"j (A.1.24) 

where A(z) is such that 

OO 

D P = in f dzz z ~ p A(.z) m j dzA[\z\)\z\~ v (A.1.25a) 

0 oo 

with 

Do — 1, D p t* 0, (p -1,2,...,*). (A.1.25h) 

Furthermore, wo assume that A[z), defined in the interval 
(0, oo), has all its derivatives and for «->■ oo they reduce to zero 
at least exponentially. 

For practically every model of an ordinary Dirac <3 function, 
we can find a model of a A function satisfying (A.1.25). If the chosen 
function A (z) does not have this property, it can he made to have 
it hy multiplying it hy z* -2 and renormalizing. For example, if 

6(e, x) mm ( 27 E)“ 3 / 2 e“ s « _1/a| * |! '“*, 


or 

A{z) « (2«)- 3/i o _1/J * 3 
does not havo tills property, then 


A (*) = -.. 

2n-2 {k+1)li r 
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does. In this case the iys are: 

Now that we have the A(z) satisfying (A.1.25), and therefore 
the model of an ordinary Dirac <5 function, we canf ind a simple 
example of the good <5(e,x) function. This is: 

8(e,x) = T e 8(s,x), (A.1.26) 

where T t is the simple operator 

*•“!(£)*'• (A - 1 - 271 

Indeed, we then, have: 


/ ^ vfc OO 

J! 8(s, x) Ixl-^dx =• — I— j e k ~ p j A(z)e 2 ~ p de 

oo 0 



1 for p — 0 , 

0 for p = 1 , 2, ..., Jc. 


(A.1.28) 


It is slightly more complicated to define d(e, x). As before, we first 
introduce a 8 function which satisfies (A.1.25). Then we define: 



(A.1.29) 

(A.I.30) 


First we can. see that T, is a generalization of T t and goes over into 
T a for co p = 0, p = 1, 2,..., Jc. 

We shall now calculate the integral 

- OO 

fdxd(e, x)|x|'* > = T,e~ p 4nf dzt?~ p A{z) = D P T,e~ p . 
oo . o 


(A..1.31) 
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Wo know that 






« for 



s > p, 


for 


Therefore 



1. 

M 






» < P- 
:(A.1.32) 

(A.1.33) 


Wo also know that 


(Ap +-* » I 0 !or 

uW ( (Aj--- j>)! 

Therefore, finally 


* < V , 

for 8 = p. 


(A.1.34) 


j dxd(e, x)|x| 
00 


\p/ \p/ A! 

1. for p — 0, 

<»„ for p 1., 2, k. 


(A.1.35) 


There in one more formal problem that we must discuss here 

and which in connected with the une of 8 (or 5) functions. We shall 
start our discussion by an example. Let us assume that in some 

calculations, operating with the (one dimensional) 8 functions we 
came across an expression 


/ 


8(a> ’)8(a") 

]P~ ■'»"[' 


fki'dx". 


(A.1.36) 


We ask: what is the meaning of such an expression l! In its place, 
lot us write: 


liin f 

PO 


\.v' xi'\ . 


dto'dx". 


(A.1.37) 
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But the expression in brackets is well-defined and, according to 
axiom (A s ), is equal to 


lim 


/ 

—oo 


<3(e, <o") 
W'\ 



k°n 

I®"! 


dx". 


(A.1.38) 


This expression, however, vanishes for a properly chosen 6. We 
also see that the result does not depend on whether we spread out 
<3(®') or <3(a>"). This is so because the function |®'— x"\~ l is symmetrical 
in x\x". 

Obviously we can generalize this theorem for an integral over 
a product of an arbitrary number of one-dimensional <3’s. The 
integral is defined if one of the <5’s is the <3 function and the rest 
of them spread out 5(e, x) functions with the same e. Thus we> have 
the following theorem: 

...,®»)$(®i). ..d{a> n )iko l ...dx n = f{ 0,..., 0) (A.1.39) 

OO 

if /(0,..., 0) is finite, and vanishes if /(0,..., 0) = oo and if the <3 
can be and are so chosen that fd(x k )f(0 ,..., x k , ...0)dx k vanishes 

OO 

for every Tc. Further generalizations of this procedure for a, say, 
twice-taken three-dimensional manifold are trivial. The rule formu¬ 
lated here will allow us to avoid all indefinite expressions appearing 
in our calculations. 

“Mass points”, “charge points”, are, of course, fictional concepts 
introduced into the theory in order to describe simply the compli¬ 
cated reality. It is just these concepts that are adequately described 
by Dirac’s ^function. 

On the other hand, if we have a continuous distribution of 
density, say, concentrated around a point a? 0 , then the space 
integrals of the products of g times |x—x 0 j -p will generally exist 
and have definite values. This is obviously true it g(x 0 ) ^ 0 and 
p < 3. The values of such integrals, corresponding to our co^’s 
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characterize to a certain extent the continuous structure of the 
density function. 

' Thus our <3 reduces such a continuous distribution to a point 
in such a way as to leave in our calculation a skeleton of the 
internal structure of the “mass points” or “charge points” in the 
form of Op . However, this can always he wiped out hy the renor¬ 
malization process. Therefore, it is more convenient always to use 

the d instead of the <3, or instead of the <3’s to avoid both the infin¬ 
ities and the renormalization procedure. 

To conclude: in this book we use this convenient tool of reason¬ 
ing, that is the <3, corresponding to the choice (o p = 0. IE nothing, 
is said to the contrary, our <3’s are the good <3’s, denoted in this 
section by <3. 


2. THE FIELD VALUES ON THE WOELD-LINES 

The solutions of field equations which we meet in this book 
are functions of world points x a = (a>°, cf ); but they depend 
functionally on motions of singularities. 

Let us assume that we have N singularities and that their mo- 

A A A 

tions are given parametrically by'£“ ?= i a (A). Here the index A 
labels the singularity (and so A = 1,2, ...N). 'The summation 

A 

over A is to be performed only if explicitly stated. By A we denote 
the invariant parameter of the A’th world-line. Therefore, if we 
denote a component of the field by /, we have: 

/-/(«“)[!, f,.•■»'?], ( A - 2 - 1 ) 

where the square bracket is reserved for functional dependence. 
Let us now consider the equation: 

- f # (A) 
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and assume that it determines i uniquely as a function of as 0 : 

2 — j^o). Then the equations of the A’th world-line can be ex¬ 
pressed in terms of the parameter <r°: 

f a = o = l, 2, 3; A= 1,2,..., .AT. 

These functions, if analytic on the x° axis, are fully determined 
by their values and those of all their derivatives at one arbitrary 
point. Therefore, the functional dependence can bo changed into 
ordinary dependence on an infinite set of arguments, that is, on all 

A A 

.the derivatives at a certain moment x°. The can bo so chosen as 
to he the same for all particles and the same as the oo° appearing 
in Therefore, we can write the / of (A.2.1) in the form: 

/ = /(»<>, f(x°), f°(A„ •••)• (A.2.2) 

Because of the approximation method used, in tho actual 
applications we shall have only time derivatives of the first and 
second order in expressions of the above type. 

The function / may (and usually will) have a singularity on 
the JL’th world-line. Tot, in spite of this, we can connect a certain 
definite meaning with the value of / along tho A’th world-lino. 
To do this we shall use our good <3’s. We define f along the world-line 

A 

A, denoted by f, by means of the equation: 

f = jdxd(x-i(x°))f{a> 0 , x)[i, |,..., £]. (A.2.3) 

Without much loss of generality, we now assume that thore 
is only one world-line, and we drop the writing of “1” above the 
We ask: what is the meaning of defining / along tho world-line, 
or, as we shall say for short, of the “tweedling” process? 

We start with a simple example: 

/ = + «0+«.(a> 8 - ft + for(® 4 - ft (a r - ft -h. • (A.2.4) 
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where the a’s are some functions of x>° only. Then using our definition 
(A.2.3), we have: 


1 — f\s — fpr — (A.2.5) 

Thus the process of tweedling consists in two things: first: in 
ignoring the singular part of /; second: in introducing the expres¬ 
sions £*( 8 °) instead of x k into the regular part of f. In this way the 
tweedled expression becomes a function of only. 

Let us now assume that / is a sum of expressions of the follow¬ 
ing type: 


1 (a? 1 - lYC® 2 - £ 2 ) 9 (a> 3 - £ 3 )* 

|x-fr ix-tr 5 * 5 


(A.2.6) 


Then we see that 

1 °: <p — 0 for n > o; 

2 °: q> = 0 if any of p, q, s is odd. This is true for arbitrary n, 
therefore also for n — 0 , because d'(x) is spherically symmetric 
in x; 

3°: If n — 0 and p, q,s are even, <p can easily be calculated 
by the use of a polar coordinate system. 

In making use of the tweedling operation prescribed by (A.2.3) 
we must be careful to distinguish between 

7\» and —- = (A.2.7) 

which we shall also occasionally denote by 95 ,*. 

These are 

<P\s — f 9V ^ ( x ~ #) d x (A.2.8) 



(A.2.9) 
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Therefore tp^ and 99 ^ are equal only if 


d(p 

w 



0 , 


(A.2.10) 


that is if tp does not depend explicitly on f®. Otherwise, we have, 
in general: 


- dtp dtp —' 

9 ’ 1 * = ~FF ~ HF 95,8 


(A. 2 . 11 ) 


One more formula which will play an important role later and 
which follows from the definition (A.2.3) is: 

- dcp d r -- '-’ 

<Pi0 = l^ = d^J = Vio+^i^io (A.2.12) 

and because £°| 0 = *° 10 = 1 , we have: 

<P[0 = ^la^io = 9 J |0+9’|*^|0- (A.2.13) 

Let us now assume the existence of two field quantities ip and tp, 
both singular at the world-line £. We ask: is it true that 

tpf = 9991 ? (A.2.14) 

This is certainly so if neither of the fields is singular. But if at 
least one of them is singular this does not need to be true. Indeed, 
the singular part of, say, ip, multiplied by the regular part of tp, 

v—' m* ~ 

may give a contribution to ftp not present in ftp, in which only the 
regular parts of both functions appear. But there is one special 
case, important in later applications, in which (A.2.14) is satisfied. 
This is the case in which the singular parts of f and tp involve only 
odd powers of |x— #|, e. g.: 


-(38+1) , I a -3 

Ix-ll 2 ®* 1 |x—f | 8 

~(2P+1) ■ | & -3 

i*-#r +i ix-#i s 


+ 

+ 


«~i 

l*-#l 

b-i 

I*—II 


+ regular part 
-(- regular part. 


(A.2.15) 
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We see, in this case, that the only possible contributions coming 
from the products of singular and regular expressions are of the 
type: 


|*-f | 2n+1 


where s+p+q = 2n+l, 


since the regular parts of y and <p contain only terms like the numer¬ 
ator of the last expression. 

But the integral of such an expression multiplied by <5(x) is zero, 
because not all of the s, p, q are even and because <5(x) is spherically 
symmetric. 

We shall see that all the expressions with which we shall deal 
in practice are of this type; that for them 


ft 


yy = yy. 


(A.2.16) 


This equality, to which we shall refer briefly as the tweedling 
of products, is assumed throughout this book. Later we shall see, 
in our concrete calculations, that this equality is indeed valid. 


3. THE COVARIANT CHARACTER OF THE <5’s. 

TENSORS ON WORLD-LINES 

We ask: what are the transformation properties of an ordinary 
Dira- ^function in four dimensions, defined by a generalization 
of the axioms J. s ) of Section 1 to this number of dimensions? 
The answer is simple enough. Obviously we must have 

J Ax (so) — 1. (A.3.1) 

This equality must have a covariant character. Since an integral 
is invariant only when the integrand is a scalar density, we deduce 
that (5( 4 ) is a scalar density, so that \t)( — g ) -1/2 is a scalar. (It would 
have been more consistent to use a different type for d. We shall, 
however, refrain from doing so, keeping in mind that <5 (4) is a scalar 
density.) 
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We assume now that |°(A) is a world-line, such that through 
f® = x° = f°(yl) the parameter X can he expressed as a function 
of »°. Consider now a tensor Tp[\\\^[X) defined only on such 
a world-line. This means: under an arbitrary coordinate transfor¬ 
mation of — x a (jx') these T’s transform according to the formula: 




duf 1 


dx'°p dx^ 


dx v P dx' fl 


ox** \| 
dx' p v 




(a.3.2) 


We may apply the rules of tensor algebra, but not those of 
tensor analysis, to such tensors defined along the world-line. To be 
able to apply the latter, we must have tensor fields. We can, at 
least symbolically, change a tensor defined along a world-line into 
a tensor field, or rather into a tensor density field. This can be done 
in the following way: 

00 

J dxs^x-sw) zfcjgw. (A. 3 . 3 ) 
—00 

Of course such a tensor density vanishes off the world-line and 
is singular on the world-line. Yet formally it can be treated as 
a tensor density field. Its transformation properties as a tensor 
density field can easily be verified from (A.3.2), from the fact 
that d (4) is a scalar density, and finally from the fact that (for 
continuous functions): 

= *<4)(»-*)/(*)• 

As a special application of this idea, we now introduce for T"*.(l) 
a vector defined only along the curve, namely dt/dX. We can change 
this vector field symbolically into a vector density field: 

00 

u“ = J dX (5( 4) (x -1) ^. (A.3.4) 

—00 

We shall be especially interested in the zero component of this 
field: 

00 

«°(») = / di°d w (x-i) = «5(,)(x-fO»°)) = d(x-§). (A.3.5) 
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This equation may be regarded as the definition of (5 (3) as far as 
its transformation properties are concerned; that is: it is the zero 
component of a vector density field. 

We have introduced the <3 (4) , which we shall meet only rarely 
in practical calculations, solely in order to find the transformation 
properties of <5 <3) . Since the <3 (8) should be one of our good <3’s we 
must redefine <3 (4) , so that the good 8 appears in (A.3.5). This is 
easy enough; the <3(4) must satisfy the axioms (J. 1 -A 8 ) for a four¬ 
dimensional manifold; besides, it must satisfy (A.3.6) where the 8 
is a good 8. Thus (A.3.5) defines both the transformation properties 
of d (4 ) and the structure of 3 (4) . 

Although <3 (3) has rather complicated transformation properties, 
we do not need to be disturbed by this since it will appear (in 
practice) always in conjunction with dx. Indeed, having a function 
/(a 0 , oo k ), multiplying it by 8 and integrating over a neighbourhood 
7(f) of the point f on the world-line, at a certain moment ® 0 , we find, 
because of (A.3.5): 

Jf8(x-§)dx=J J fdS°8 w (w-£)dx =/(f°, r*). (A.3.6) 

m r(f) -oo 

Putting f — 1, we find: 

Ja(x-f)dx = l (A.3.7) 

which means that ddx can be treated as an invariant. 

Since we have the connection (A.3.5) between d ^ and d, we 
can write down in a different way the equations (A.3.3): 

—00 

OO 

J «%(•-««•» 

— 00 


(A.3.8) 
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The last equation, shows that the singularity of the symbolic 
tensor density field £ is a singularity of the <5 type along the world¬ 
line. One can deduce from it the opposite formula expressing 
a tensor along a world-line in terms of the tensor density field: 

J< A - 3 - 9 > 

a 

i 

where a is a hypersurface on which the point = £“(4) lies. 

These symbolic fields, defined by (A.3.3), will play an essential 
role in our analysis. Such fields conserve formally the usual prop¬ 
erties of the ordinary tensor densities, with respect to parallel 
displacement and covariant differentiation. 
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